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Abstract 


Research into micro energy harvesting systems has been driven by the need 
for self-sustaining, stable power sources for interconnected wireless sensors. 
Waste heat, particularly at temperatures below 200 °C, presents a promising 
but challenging energy source to recover using current conversion technolo- 
gy. The progress of high-performance thermomagnetic generators (TMGs) 
has been impeded by the lack of advancements in thermomagnetic materials. 
This study builds upon previous research on small-scale TMGs. The main 
objectives are; 


e Development of a Lumped element model (LEM) simulation of the 
TMG to analyze and optimize the performance. 


e Utilizing LEM and experiments to understand the effect of various de- 
sign parameters. 


e Upscaling of the TMG’s active material volume to increase the abso- 
lute power output of a single device. 


e Upscaling of TMG by operating multiple TMGs in parallel to increase 
the lateral size. 


e Extending the operation range of heat source temperatures to near 
room temperatures without losing resonant self-actuation. 


Experimental measurements and LEM simulations explore TMGs based on 
various materials such as the Ni-Mn-Ga Heusler alloy film, Gadolinium and 
La-Fe-Si-H. The effect of various design parameters on the performance of 
the TMG is investigated. LEM simulations describe the coupled dynamic 
properties of TMGs utilizing magnetic shape memory alloy (MSMA) films. 
The TMG employ self-actuation by leveraging temperature-dependent 
changes in magnetization and rapid heat transfer through thermomagnetic 
thin films. 


Abstract 


Detailed LEM simulations examine the temperature change, magnetic field 
change, and resulting magnetization of the TM films over time and position. 
Optimal conditions for resonant self-actuation are achieved through careful 
device parameter design, leading to continuous, undamped oscillation of the 
cantilever. The thesis discusses various design parameters affecting resonant 
self-actuation for the case of Ni-Mn-Ga thin films, emphasizing the im- 
portance of fine-tuning each parameter for maximum output power. The 
impact of factors like magnet, coil turns, cantilever stiffness, load resistance 
(Rx), Curie temperature (Tc), heat transfer coefficient (hr), and thermal 
resistance (Rv) are investigated to understand their influence on the TMG 
performance. LEM simulations reveal critical values of hr and Rp» that enable 
stable energy generation with significant stroke and frequency, resulting in 
substantial electrical output. 


Upscaling the TMG with Ni-Mn-Ga thin film demonstrates opposing effects 
on power output and footprint, with an improved electrical power per foot- 
print achieved by increasing film thickness from 5 to 40 um. At a tempera- 
ture change of only 3 °C and a frequency of 146 Hz, values of 50 uW/cm? are 
achieved. The parallel architectures are essential for generating sufficient 
energy for direct application. Thermal cross coupling affects the dynamic 
performance and power output of TMGs operating in parallel. Notably, 
thermal effects become noticeable at low device separations and high heat 
source temperatures, but no magnetic or mechanical interaction is observed 
between the parallel-operating TMGs. 


Using gadolinium as the active TM film, near-room temperature operation is 
possible at a low heat source temperature (Tsource) of 40 °C. The TMG is able 
to generate 1.3 uW of power at a frequency of 54 Hz at this Tsource, equivalent 
to a power output per footprint of 10 uW/cm?. At a Tsource of 65 °C, this value 
sharply rises to 24 1. W/cm? with an oscillation frequency of 117 Hz. Even 
though Tamb of 11 °C is required for this optimum performance, the device 
can sustain resonant self-actuation up to an ambient temperature (Tamb) of 
19°C while still producing 8.7 uW/cm? of power for a Tsource of 50 °C. 
Moreover, the sharp boundaries of operation temperatures in terms of Tsource 
and Tamb are investigated and presented. The TMG using hydrogenated La- 
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Abstract 


Fe-Si alloys as active TM film can generate 38 uW/cm? from a Tsource of 
90 °C while operating at a frequency of 137 Hz. TMGs are shown to be able 
to function across a broad range of Tsource temperatures while producing 
significant power output per device footprint. Thus, making them extremely 
competitive not just with other TMGs but also with the best commercially 
available thermoelectric generators. 
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Kurzfassung 


Die Forschung im Bereich der Mikro-Energiegewinnung wurde durch den 
Bedarf an autarken sowie stabilen Energiequellen für vernetzte und drahtlose 
Sensoren vorangetrieben. Abwärme, insbesondere bei Temperaturen unter 
200 °C, stellt eine vielversprechende, aber mit den derzeitigen Umwand- 
lungstechnologien schwer zu gewinnende Energiequelle dar. Der Fortschritt 
von thermomagnetischen Generatoren (TMGs) mit hoher Leistung wurde 
durch den Mangel an Weiterentwicklungen von thermomagnetischen Materi- 
alien behindert. Diese Arbeit stützt sich auf frühere Forschungsarbeiten zu 
TMGs im kleinen Maßstab. Die Hauptziele sind: 


e Entwicklung eines LEM-Modells (Lumped Element Model) zur Simu- 
lation des TMG, um die Leistung zu analysieren und zu optimieren. 


e Nutzung von LEM und Experimenten, um die Auswirkungen ver- 
schiedener Designparameter zu verstehen. 


e Die Hochskalierung des Volumens des aktiven Materials eines TMG, 
um die absolute Ausgangsleistung eines einzelnen Generators zu er- 
höhen. 


e Die Hochskalierung des TMG durch Parallelbetrieb mehrerer TMGs 
zur Vergrößerung der lateralen Größe. 


e Erweiterung des Betriebsbereichs der Wärmequelle auf Temperaturen 
nahe der Raumtemperatur, ohne die resonante Selbstaktivierung zu 
verlieren. 


Zunächst werden mittels experimenteller Messungen und LEM-Simulationen 
TMGs, die auf verschiedenen Materialien wie dem Ni-Mn-Ga Heusler- 
Legierungsfilm, Gadolinium und La-Fe-Si-H basieren, grundsätzlich er- 
forscht. Die Auswirkung verschiedener Designparameter auf die Leistung des 
TMGs wird untersucht. Dabei beschreiben LEM-Simulationen die gekoppel- 
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ten dynamischen Eigenschaften von TMGs, die Filme aus magnetischen 
Formgedächtnislegierungen (MSMA) verwenden. Die TMG nutzen Selbstak- 
tivierung, indem ein temperaturabhängige Magnetisierungsänderungen und 
einen schnellen Wärmetransfer durch thermomagnetische Dünnschichten 
ausgenutzt wird. 


Detaillierte LEM-Simulationen zeigen die Temperaturänderung, die Magnet- 
feldänderung und die daraus resultierende Magnetisierung der TM-Filme 
über Zeit und Position. Optimale Bedingungen für eine resonante Selbstakti- 
vierung werden durch sorgfältiges Design der TMG-Parameter erreicht, was 
zu einer kontinuierlichen, ungedämpften Oszillation des TMG-Ausleger 
führt. In dieser Arbeit werden verschiedene Design-Parameter erörtert, die 
sich auf die resonante Selbstaktivierung im Falle von Ni-Mn-Ga- 
Dünnschichten auswirken, wobei die Bedeutung der Feinabstimmung jedes 
Parameters für eine maximale Ausgangsleistung hervorgehoben wird. Die 
Auswirkungen von Faktoren wie Magnet, Spulenwindungen, Auslegerstei- 
figkeit, Lastwiderstand (RL), Curie-Temperatur (Tc), Wärmeübergangskoeffi- 
zient (hr) und Wärmewiderstand (Rb) werden untersucht, um ihren Einfluss 
auf die TMG-Leistung zu verstehen. LEM-Simulationen zeigen kritische 
Werte für hr und Ro, die eine stabile Energieerzeugung mit signifikantem Hub 
und Frequenz ermöglichen, was zu einer deutlichen Steigerung der elektri- 
schen Leistung führt. 


Die Hochskalierung des TMG mit Ni-Mn-Ga-Dünnschicht zeigt gegensätzli- 
che Auswirkungen auf die Leistungsabgabe und die Grundfläche, wobei eine 
verbesserte elektrische Leistung pro Grundfläche durch eine Erhöhung der 
Schichtdicke von 5 auf 40 um erreicht wird. Bei einer Temperaturänderung 
von nur 3 °C und einer Frequenz von 146 Hz werden Werte von 50 uW/cm? 
erreicht. Die parallelen Architekturen sind entscheidend für die Erzeugung 
ausreichender Energie für die direkte Anwendung. Die thermische Kreuz- 
kopplung beeinträchtigt die dynamische Leistung und die Leistungsabgabe 
von parallel betriebenen TMGs. Thermische Effekte machen sich vor allem 
bei geringen Abständen zwischen den Bauelementen und hohen Temperatu- 
ren der Wärmequelle bemerkbar, wobei jedoch keine magnetischen oder 
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mechanischen Wechselwirkungen zwischen den parallel arbeitenden TMGs 
beobachtet werden. 


Bei Verwendung von Gadolinium als aktiver TM-Schicht ist ein Betrieb bei 
einer niedrigen Wärmequellentemperatur (Tsource) von 40 °C möglich. Der 
TMG kann bei dieser Tsource eine Leistung von 1,3 uW bei einer Frequenz 
von 54 Hz erzeugen, was einer Ausgangsleistung von 10 uW/cm? pro Fläche 
entspricht. Bei einer Tsource von 65 °C steigt dieser Wert sprunghaft auf 24 
uW/cm? bei einer Frequenz von 117 Hz an. Obwohl für eine optimale Leis- 
tung eine Tamb von 11 °C erforderlich ist, kann das Bauelement die resonante 
Selbstaktivierung bis zu einer Umgebungstemperatur (Tamb) von 19 °C 
aufrechterhalten und dabei immer noch 8,7 uW/cm? Leistung bei einer Tsource 
von 50 °C erzeugen. Außerdem werden die scharfen Grenzen der Betriebs- 
temperaturen in Bezug auf Tsource und Tamb untersucht und vorgestellt. Ein 
TMG, bei den hydrierten La-Fe-Si-Legierungen als aktiven TM-Film ver- 
wendet, kann 38 uW/cm? aus einer Tsouree von 90°C erzeugen, während es mit 
einer Frequenz von 137 Hz arbeitet. 


xi 


Abbreviations and parameters 


LEM 
TMG 

TM 

Tamb 
Tsource 
SMA 

Te 

c axis 
Ni-Mn-Ga 
Gd 
La-Fe-Si 
La-Fe-Si-H 


La-Fe-Mn-Si- 
H 


TM film 
SmCo 
NdFeB 


TM cycle 


Lumped element modeling 
Thermomagnetic generators 
Thermomagnetic 

Ambient temperature 

Source temperature 

Shape memory alloy 

Curie temperature 

Magnetic easy axis 

Nickel manganese gallium alloy 
Gadolinium 

Lanthanum iron silicon alloy 


Hydrogenated lanthanum iron silicon alloy 


Hydrogenated lanthanum iron manganese silicon al- 


loys 
Thermomagnetic material 
Samarium Cobalt magnets 
Neodymium Iron Boron magnets 


Thermomagnetic cycle 
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Abbreviations and parameters 


PDE - Partial differential equation 

hf - Contact heat transfer coefficient 
Tfilm - Temperature of film 

Rb - Bonding layer thermal resistance 
RL - Load resistance 

AT - Temperature change in TM film 
Tcant - Cantilever temperature 
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1 Introduction 


With the increase in IoT (Internet of Things) devices and connected technol- 
ogy, there is a growing need for energy in remote and hard-to-access places. 
Health monitoring of buildings by wireless sensor networks embedded inside 
constructions, wearable devices that are embedded with sensors to monitor 
human health, implantable devices, and smart home sensors are a few exam- 
ples of devices that have such an energy demand.'? These devices usually 
require electrical power in the range of tens of uW up to 100 mW.? Such high 
demands for remote access to energy are currently met by batteries but these 
pose a great disadvantage as it has limited storage capacity which means they 
have to be either recharged or replaced at regular intervals.* Thus, new 
innovative solutions to generate energy in situ from environmental sources 
are required if the world is moving towards a connected world which increas- 
es energy demand and at the same time needs to keep the goal of carbon 
neutralization in the near future.” 


1.1 Energy harvesting 


Energy harvesting is a term used to define the collection of ambient energy in 
small amounts (in range of uWs to mWs) to power a remote wireless device 
where using batteries or wires is impractical. The principle is similar to large- 
scale renewable energy generation but at miniature scale.° Energy harvesters 
generate very low amounts of power at miniature scale making it attractive 
for inaccessible devices and networks with a large number of nodes.’ The 
block diagram shown in Figure 1 illustrates a complete energy harvesting 
circuit. When an energy harvester or power generator is mentioned in this 
thesis, an energy converter is referred. There are several forms of ambient 
energy sources like mechanical energy in the form of vibrations, solar energy, 
wind energy, radio frequencies, and finally thermal energy.* 
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The load system refers to a microsystem that needsvery little amount of 
energy to operate like a wireless temperature sensor or a strain gauge. Most 
of these devices are on standby and only have an activity mode for a few 
milliseconds at regular intervals where it powers up, takes a reading, sends it 
wirelessly to a hub, and goes back to standby mode with very low power 
consumption.'® 


ey => (corre Energy => Input = Energy => (noun) Output = (ce) Load 
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Figure 1 Energy harvesting circuit adapted from Bader et al.” 


Kinetic power generation 


Kinetic energy harvesters consist of devices that convert kinetic energy like 
mechanical vibrations or movements into electrical energy.®!! Vibrational 
devices convert ambient vibrations into electrical energy using electromag- 
netic induction” or by straining piezoelectric crystals.!? There are piezoelec- 
tric energy harvesters that directly convert a mechanical strain into electrical 
energy.!*!5 Typical vibrational energy harvester contains a cantilever beam 
with an attached tip mass, which vibrates in response to ambient vibrations. 
The device only performs with high efficiency when the device resonance 
matches the ambient vibrations.'° But there are also other forms of energy 
harvesters that harvest vibrational ambient energy.!” 


1.1.1 Solar power generation 


Solar energy conversion by utilizing photoelectric conversion is one of the 
most common ways of energy harvesting.'®:!° Photovoltaic power generation 
is scalable to very large scales and there are several large solar power plants 
currently operating without output range of Megawatts around the world.” 
One of the advantages of a solar energy harvester is that it generates a DC 
output that can directly be utilized by microsystems. The disadvantage is that 
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sunlight is only available during the day and inaccessible areas might not 
have adequate access to sunlight. 


1.1.2 RF power generation 


This is an interestingly popular research topic that has gained attention after 
the introduction of RFID technology (Radio Frequency Identification). But 
due to impedance mismatch and losses (e.g., conduction, thermal, reflective, 
conversion etc.) that affect the performance of an RF energy harvesting 
device, it is not widely adopted.?! An RF harvesting device converts electro- 
magnetic energy from the ambient into usable electrical energy in the form of 
DC voltage. One of the advantages of RF energy is that this is not limited by 
location or time of the day and is available outdoors and indoors throughout 
the day in most locations around the world.” 


1.2 Thermal power generation 


Thermal energy is abundantly available as ambient waste energy. Most of 
this energy is available in the low temperature regime of below 200 °C which 
falls below the range of normal power plants and conventional energy gen- 
eration cycles.7*> This is of particular interest as a large proportion of 
energy consumed in the industrial sector and urban areas is rejected as waste 
heat. There are several forms of thermal energy harvesting principles, a few 
relevant principles will be introduced. 


Thermoacoustic engines utilize sound waves to pump heat or use heat to 
pump sound waves. Typically, a gas-filled acoustic resonator is positioned 
appropriately with thermoacoustic elements.” Thermogalvanic energy 
harvesting operates based on the dependence of electrode potential on tem- 
perature. An electrochemical cell is used to charge and discharge at different 
temperatures to create a cycle that will convert heat into electricity.” Howev- 
er, as the temperature gradient becomes small, the conversion of energy 
becomes inefficient.” The current state-of-the-art technology for energy 
harvesting at a miniature scale is thermoelectric generation using the Seebeck 
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effect.” >! Thermoelectric devices can generate an electric potential from a 
constant temperature gradient. These devices suffer from scaling restrictions 
and low efficiency at low temperature differences.’ Disadvantage of thermo- 
electric devices is the requirement for a large heat sink which usually is 
several times larger than the device itself to maintain a reasonable tempera- 
ture gradient required for energy generation. There are also other forms of 
thermal energy harvesters like pyroelectric energy harvesting which utilizes 
the pyroelectric effect.’ Disadvantage is that it requires a temperature 
change (AT) in the material to generate electricity, unlike TEG which only 
requires a constant gradient. 


Scaling requirements 


The energy harvesting devices typically need to be scaled according to the 
power requirement of the microsystem. Hence scalability is one of the major 
criteria that need to be fulfilled when considering an energy harvesting 
device.°’” Thermoelectric devices are not scalable to smaller dimensions due 
to the inability to maintain a temperature gradient without the use of a larger 
heat sink.*? Thus, a sustainable energy harvesting device should possess an 
active heat dissipation mechanism that can maintain highly efficient power 
conversion even when scaled to miniature dimensions. 


1.3 Thin film microsystems 


Miniaturization is considered a fundamental element of technological ad- 
vancement, since more devices are being made in miniature scales to increase 
portability, decrease cost, and reduce weight.** The development of MEMS 
(micro-electro-mechanical systems) has pushed the technology to downsize 
rapidly. MEMS devices tend to be highly durable, portable, and only require 
a very small amount of energy to operate.*? Thin film actuation plays an 
important role in microsystems developments due to the limited number of 
actuation mechanisms available at this dimensional scale. The major actua- 
tion mechanisms at the microscale include electrostatic and thermal micro- 


actuators. SMA actuators,*%*! microscanners” and magnetic thin film actua- 
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tors? based on thermomagnetic (TM) actuation utilize thin films to achieve 
small-scale actuation. Thin film devices have the advantage of a very large 
surface-to-volume ratio enabling a high frequency of operation due to the fast 
heat transfer. Thin films are usually fabricated by either vapor deposition, 
cold rolling, or melt spinning. More details about thin film fabrication are 
detailed in section 2.7. Thin film fabrication is also quite straightforward, and 
most structures can be realized using laser machining or UV lithography 
followed by wet etching.***5 


1.3.1 Smart materials 


Thin film actuators and energy harvesters are fabricated out of smart materi- 
als. These are categorized as materials that exhibit multifunctional properties. 
Such materials show a significant change in their properties in response to an 
external stimulus in the form of a thermal, magnetic, electric, or mechanical 
stress field. The majority of smart materials are ferroic or multiferroic mate- 
rials, which exhibit several ferroic properties like ferroelectricity, ferromag- 
netism or ferroelasticity. Such materials show a coupling between the prima- 
ry order parameters and the multiferroic property, which allows the material 
to perform smart actuation.*® 


Figure 2 shows the coupling between thermal, magnetic, and mechanical 
properties of a multiferroic material. In this case, by changing the magnetic 
field, the magnetization of the material could be changed but at the same 
time, a strain can also be achieved in the material. Heusler alloys are a prime 
example of smart materials.“ These are ordered intermetallics that exhibit 
multiferroic transformations as electronic spin interactions strongly depend 
on interatomic distances and local symmetry. Some of the examples of 
Heusler alloys are Ni-Mn-Z and Ni-Co-Mn-Z (Z=Ga, In, Sn, etc.) that exhibit 
giant effects including magnetic-field induced strains,** magnetoresistances,”” 
magnetocaloric effects*’ and magnetic field induced shape memory effects.°! 
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Magnetoelasticity 
Elastomagnetics 


Figure 2 Mechanical, magnetic, and thermal coupling in multiferroic materials.” 


1.4 Aim ofthis work 


The overall goal of this work is to develop novel miniature thermomagnetic 
generators (TMGs) operating in resonant self-actuation mode that are scala- 
ble in terms of size and power output to suit the application. The are several 
objectives listed as follows: 


e Using low-grade waste heat below 200 °C for power generation 


e To develop a LEM model to represent the multiphysical characteris- 
tics of the TMG operation, so as to optimize and maintain resonant 
self-actuation, and large power output when the various device pa- 
rameters are changed. 


e To use LEM to analyze the thermal processes in the TMG. 


e To utilize experimental measurements and LEM simulations to devel- 
op in depth understanding of effect of various TMG design parame- 
ters. 


1.4 Aim of this work 


To upscale the TMG’s active material volume to increase the power 
output of a single device and to understand the effects of scaling on 
performance. 


To upscale the TMG in lateral size by operating TMGs in parallel to 
increase the overall output from the array and to understand the cou- 
pling effects. 


Maintain resonant self-actuation that results in large power output 
over a wide range of operating sources and ambient temperatures. 


To develop TMG’s operating near room temperature while still main- 
taining resonant self-actuation even though the temperature change 
achieved in the active material volume would be small. 


1.4.1 Outline of this document 


Section 2 describes the relevant background information regarding TM 
energy generation, fabrication methods, and measurement techniques. 


Section 3 presents different thermomagnetic generators that have been 
demonstrated in literature over the years. The previous generation of TMGs 
from the Smart Materials and Devices group in KIT is also presented, based 
on which this work has been completed. 


Section 4 introduces the Lumped element model utilized to analyze the TMG 
performance and the thermal processes within the TMG. The operation 
principle of the TMG and results from the TMG using Ni-Mn-Ga is present- 
ed. 


Section 5 discusses various parameters that have an impact on the perfor- 
mance of the TMG and ability to sustain resonant self-actuation. LEM 
simulation results for parameter study are presented. 


Section 6 describes the upscaling of the TMG when the active material 
volume of the device is increased and when the lateral size is increased by 
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multiple devices operating at close distances. The increase in power output 
and the coupling effects of multiple device operation are presented. 


Section 7 describes TMGs operating near room temperature with Gadolinium 
film and La-Fe-Si film as active TM film materials. The performance of the 
TMG using Gd film and La-Fe-Si-H films are investigated. The boundary 
conditions observed are then analyzed. 


Section 8 discusses the effects of different parameters and the performance of 
different TMGs. 


Section 9 summarizes the presented results of the thesis. 


2 Basic principles 


The most important principles like ferromagnetism, shape memory effect, 
power conversion and heat transfer, and LEM modelling information that are 
relevant to this thesis will be described in this section. 


2.1 Magnetism 


There are several types of magnetism discovered and explained in physics 
like ferromagnetism, paramagnetism, diamagnetism, and antiferromagnetism. 
But ferromagnetism is the strongest type of magnetism of the group. When 
exposed to a magnetic field, ferromagnetic materials become highly magnet- 
ized. The magnetic moments of the atoms within ferromagnetic materials 
align in the same direction, producing a net magnetic field that is significant- 
ly larger than the applied magnetic field. Iron, cobalt, gadolinium, and nickel 
are some of the materials that are ferromagnetic in an elemental state. Ferro- 
magnetism also depends on the lattice structure. A combination of non- 
magnetic materials could be ferromagnetic like Heusler alloys if they contain 
at least one ferromagnetic element, and a combination of magnetic materials 
could be non-ferromagnetic like stainless steel.°>> 


2.1.1 Ferromagnetism 


An electron carries a charge but also has a magnetic dipole moment. A 
magnetic dipole moment is a measure of the strength and orientation of a 
magnetic dipole. A magnetic dipole is a fundamental entity that consists of a 
current loop, or a pair of opposite magnetic poles separated by a small dis- 
tance. When a magnetic dipole is placed in an external magnetic field, it 
experiences a torque that tends to align the dipole with the field. The electron 
acts like a small magnet producing its own magnetic field because of the 
quantum mechanical spin. This can either be pointing up or down depending 
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on its state and is the source of ferromagnetism. The contribution from the 
orbital angular momentum of the electron when revolving around the nucleus 
also adds to the ferromagnetism. 


For a magnetic dipole consisting of a pair of opposite magnetic poles separat- 
ed by a distance d, the magnitude of the magnetic dipole moment is given 
by:°° 


m = IA.ê 


Eq (1) 


where m is the magnetic dipole moment vector, I is the current flowing 
through the loop, and A is the area of the loop and ê is the unit vector corre- 
sponding to the direction. 


|m| = IAsin(0) = ud 


Eq(2) 


where u is the magnetic moment per unit volume and 9 is the angle between 
the magnetic dipole moment vector and the direction of the external magnetic 
field, while d is the distance between magnetic poles. 


In a macroscopic view, a high field would be generated when magnetic 
dipoles of electrons in a piece of material are aligned in the same direction. 
But the electrons are always present in pairs of opposite spin, hence an atom 
with a filled shell would mean a total dipole moment of zero. But the material 
in which atoms have partially filled shells would have a net dipole moment. 
By Hund’s rules,* the first few electrons in a shell tend to have the same spin 
and hence the net dipole moment tends to be large. Exchange interaction, also 
known as exchange coupling, is a type of interaction between magnetic 
moments in a material. It results from the transfer of electrons between 
nearby atoms or ions, which can cause the alignment or anti-alignment of 
nearby magnetic moments. One common way to represent exchange interac- 
tion in a spin model is through the Heisenberg Hamiltonian equation, which 
is given by:” 
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Figure 3 Schematic diagram of magnetic dipole moment arrangement in various type of 
magnetic materials. Adapted from Hui et al.*’ 


H = -X Jy Si - Sj- uB Y Si 
Eq (3) 


where Si and Sj are the spin vectors at sites i and j, respectively, % < i, 
represents the sum over all nearest-neighbor pairs of spins, and J is the 
exchange interaction energy. The exchange interaction energy J in this 
Hamiltonian equation reflects the intensity of the contact between adjacent 
spins. Si. Sj shows whether the interaction is ferromagnetic (positive number) 
or antiferromagnetic (negative value). The term "magnetic anisotropy" 
describes how a material's magnetic characteristics vary depending on its 
direction. It specifically refers to a material's propensity to exhibit various 
magnetic characteristics (such as magnetization or susceptibility) depending 
on the direction or orientation of the measurement. Although, exchange 
interaction keeps the magnetic dipoles aligned but not necessarily in any 
particular direction due to magnetic anisotropy.” Figure 3 shows the mag- 
netic dipole moment orientation in various materials. Ferromagnetic materi- 
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als are usually divided into small regions or areas called magnetic domains as 
the exchange interaction is a short-range force. In other words, it is a small 
region within a material where the magnetic field created by the electrons in 
the atoms all point in the same direction. While, within these domains, the 
material tends to have spins aligned in one direction, these domains being not 
all aligned together would cause them to cancel each other and the effective 
magnetization would be either low or zero since the bulk material usually is 
found in its lowest energy state. These domains are not fixed in any way and 
can be altered by using an external magnetic field. When exposed to high 
external fields, the magnetic domains tend to reorient themselves to keep 
more of the dipoles in the same direction as the external field. This is called 
the magnetization of the material. Total Magnetization M of the material is 
given by Eq ( 4). 


M= Nyoı < M > 


Eq (4) 


Usually, the material tends to return to its unmagnetized state once the 
magnetic field is removed. Materials that retain a large portion of magnetiza- 
tion upon removal of the external magnetic field are hard magnetic materials 
and are used to make permanent magnets which are used in induction motors 
and other magnetic applications. The remaining magnetic field in the material 
is known as the remanent magnetic field. 


2.1.2 Curie temperature 


First-order changes are distinguished by an abrupt shift between the two 
phases, during which the substance experiences a discontinuous change in its 
characteristics across the phase border. These transformations generally entail 
a shift in the material's internal energy and are followed by heat release or 
absorption. The melting of a solid to create a liquid is an established example 
of a first-order transformation. The solid experiences a first-order phase 
transition and melts into a liquid at the melting point, resulting in a signifi- 
cant change in the material's physical characteristics such as density, viscosi- 
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ty, and thermal conductivity. In contrast to first-order transition, second-order 
transitions are distinguished by a gradual shift in the properties of the sub- 
stance rather than a sharp transition between discrete stages. These changes 
frequently involve a shift in the material's internal energy but do not usually 
involve the release or absorption of heat.°'' When the temperature of a 
material increases, random thermal motion, and entropy of material causes 
reorientation of dipoles. When the temperature increases beyond a specific 
point (specific to the material), called Curie temperature (Te), the material 
undergoes a second-order transformation losing its ability to maintain spon- 
taneous magnetization. The material will still paramagnetically respond to 
external fields but will not attract to it like in the case of a ferromagnetic 
state. Te is a critical point at which the magnetic susceptibility is infinite in 
theory. Being a second-order transformation, the material transformation at 
Te is not associated with any hysteresis or latent heat. 


2.2 Shape memory effect 


2.2.1 The one-way shape memory effect 


The one-way shape memory effect is the ability of a material to restore its 
original shape or configuration above a characteristic temperature. Shape 
memory alloys are a group of material alloys that exhibit shape memory 
effect. This ability has made SMAs interesting for a wide range of applica- 
tions. The shape memory effect is shown in the stress-strain curve of a typical 
SMA material in Figure 4. Let us assume a SMA material that is currently at 
a sufficiently low temperature denoted by point 1 in Figure 4. Martensitic 
phase transformation is a first-order solid-state phase transformation in which 
a parent phase changes crystal structure into a new phase. External variables 
such as mechanical deformation, thermal shifts, or chemical preparation can 
all cause this transition. The martensitic transition usually includes a rapid 
and irreversible shift in the material's crystal structure, resulting in a new 
metastable phase. The new low temperature phase, known as martensite, has 
a twisted and/or twinned crystal structure and can have distinct mechanical, 
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electrical, and thermal characteristics. The martensite has a low symmetry 
lattice structure which allows for twin boundaries to be present, due to 
various orientations of the variants. Twin boundaries are interfaces or planes 
within a crystal where the orientations of two areas of the crystal lattice are 
different. These areas, known as twin domains, are mirror copies of each 
other across the twin boundary. The twin border is a narrow planar defect that 
divides the twin domains. These twin boundaries can move easily and cause 
reorientation of martensite. When this material is loaded, the initial response 
from the material is elastic and the material could be assumed to be in a phase 
of twinned martensite, where all several martensite variants exist together. 
Once the strain and stress reach point 2 in Figure 4, the material exhibits 
reordering of martensite variants and the response is similar to plastic defor- 
mation. The reordering is completed once the material reaches point 3 in 
Figure 4, after which the material behaves elastically when a further load is 
applied up to point 4. 


cooling 


Lo) esting 


"RE 


Figure 4 Shape memory effect, the stress-strain-temperature diagram of the one-way effect of 


shape memory alloys (SMAs).‘° 
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The material has now been de-twinned completely by reorientation of various 
variants along the direction of the dominant variant. From point 4 to point 5 
in Figure 4, the material is unloaded, and the resulting behavior exhibited by 
the material is elastic. But once it is fully unloaded a residual deformation 
remains in the material. This residual deformation could be reversed by 
heating the material above a specific temperature, denoted by point 6. After 
point 


6 is reached the material experiences a phase change into de-twinned austen- 
ite. The lattice structure of austenite has high symmetry and thus only a 
single shape is possible at this high-temperature state. Further increase in 
temperature causes the material to completely transform into austenite at 
which point the material has also recovered its original shape completely. 
Cooling the material does not have any effect on the material shape and the 
material transforms into twinned martensite. This completes a closed loop 
and is also known as the one-way effect. 


2.2.2 Superelasticity 


The free energy equilibrium of martensite and austenite affect the transfor- 
mation temperature of SMA. Superelasticity occurs when there is a change in 
external energy of the material by a mechanical load resulting in a shift of 
transformation temperatures. This effect is shown in Figure 5. When stress is 
applied to SMA in its austenite state, the material behaves elastically up to 
critical stress. 


At critical stress, the material starts transforming to martensite, displaying a 
large change in strain and if further stress is applied, then the material trans- 
forms completely to martensite releasing heat in the process and elastic 
behavior reappear after transformation. When the load is released, the materi- 
al transforms back to austenite by absorbing energy, and the original shape is 
recovered. The large stress-strain plateau is shown in Figure 5. The heat 
release and heat absorption of material on loading and unloading of the 
material form the basis of elastocaloric effects. 
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Figure 5 SMA superelastic effect stress-strain characteristics.°* 


2.3 Magnetic shape memory effect 


In the conventional shape memory effect explained in section 2.2, the materi- 
al has high-temperature austenite and a low-temperature martensite phase. 
The switching between these two phases could either be driven by a change 
in temperature or by mechanical deformation. In Magnetic shape memory 
alloys, a magnetic field could be applied to initiate a field induced reorienta- 
tion. When an external magnetic field is applied to the MSMA material, the 
magnetic domains reorient themselves in a particular direction (field induced 
reorientation), resulting in a shift in the shape or mechanical characteristics of 
the material. A typical schematic of a single crystal Ni-Mn-Ga is given in 
Figure 6. It can be seen that on the application of load, causing deformation 
in the material, the crystallographic variants are oriented along the direction 
of load with the magnetic easy axis (c axis) parallel to the force direction. 
Given its crystallographic structure or anisotropy, the magnetic easy axis is 
the direction in which the magnetic moment of the magnetic material chooses 
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to position itself. It is the path along which the magnetization of the material 
needs the least amount of energy for alignment. The material remains in this 
state and the magnetic-field-induced-strain could be recovered by applying a 
magnetic field perpendicular to the c axis, which causes the magnetic do- 
mains to realign and results in a shift in the material's shape due to the mag- 
netostriction effect. This reorients the c axis in the direction of the magnetic 
field. These transformations or reorientation of variants occurs in the material 
by a movement of twin boundaries. In contrast the magnetic properties of the 
material changes as the magnetic moments of the material are reoriented 
when it is subjected to a mechanical stress or strain. This inverse magneto- 
striction effect is called Villari effect and allows to control the magnetic 
properties of the material by using mechanical deformation. 


Deformation Magnetitic Field 


\7 
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Figure 6 Schematic diagram of the tetragonal structure of Ni-Mn-Ga single crystal showing 
reorientation.* 
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2.3.1 Ferromagnetic shape memory alloys 


The discovery of large magnetically induced strains in Ni-Mn-Ga alloys by 
Ullakko et al.** showed that magnetostriction as a result of magnetic field 
induced reorientation takes place in such alloys up to a few percent have 
sparked high attention to such alloys. These belong to the family of Heusler 
alloys that exhibit multiferroic phase transformations.**°!6 Such material 
is ferromagnetic in both martensitic austenitic phase and undergoes a second- 
order transformation at Tc, after which the magnetization of the material 
abruptly drops and becomes paramagnetic. Figure 7a shows the Ferromag- 
netic-paramagnetic phase transformation in Ni-Mn-Ga. The initial small 
change in magnetization denotes the first-order transformation in the materi- 
al. The Differential scanning calorimetry (DSC) measurement curve in 
Figure 7b shows the martensitic start and finish temperatures, as well as the 
austenitic start and finish temperatures. 
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Figure 7 Material characterization of Nis3.s-Mn23.8-Ga22.7 alloy, (a) Magnetization measurement 
with the second-order transformation, (b) Differential scanning calorimentry 
(DSC) measurement showing the first-order transformation. Adapted from Jo- 
seph et al.” 


2.3.2 Metamagnetic shape memory alloys 


Metamagnetic shape memory alloys are a subset of magnetic shape memory 
alloys that are not ferromagnetic shape memory alloys. Ni-Mn-X (X=In, Sb, 
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Sn) is a representative material of the category.°® These materials tend to 
have para or antiferromagnetic martensitic phase and the austenitic phase is 
ferromagnetic below the Tc. Such materials have a large rapid change in 
magnetization at the phase transformation, this can be observed in Figure 8b. 
But in this case, the material has latent heat and hysteresis associated with the 
transformation, being a first-order transformation. 
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Figure 8 (a) DSC measurement and (b) Magnetization measurement data for a sample of 
Ni4sCosMnaoSnıo.°” Reprinted with permission from AIP Publishing. 


The quantity of heat energy absorbed or released during the phase transfor- 
mation of a material without producing a change in its temperature is referred 
to as latent heat. While the phenomenon of hysteresis describes how a mate- 
rial's response to a specific stimulus, such as temperature, pressure, or mag- 
netic field, relies on its history of exposure to that stimulus. The hysteresis 
can be tuned by changing the lattice parameters of martensite and austenite 
unit cells.’ Figure 7a shows the DSC measurement which shows the latent 
heat associated with the transformation and Figure 8b shows the magnetiza- 
tion data plotted against temperature. The austenite phase experiences a 
reduction in magnetization with an increase in temperature at Tc. 
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2.4 Energy conversion to electrical energy 


2.4.1 Electromagnetic induction 


The physical interaction that is observed in electrically charged particles 
moving through a magnetic field is commonly called electromagnetic force 
or Lorentz force. This has magnetic and electric components and plays an 
important role in determining the force on a moving point charge in an 
electromagnetic field. The magnitude of the Lorentz force on a charged 
particle is provided by the equation:*? 


F = q(v XB) 


Eq(5) 


Where F represents the force, q represents the particle's charge, v represents 
its velocity, and B represents the magnetic field. The right-hand formula 
determines the force orientation, which is perpendicular to both the velocity 
vector and the magnetic field vector. The process of producing an electric 
current in a circuit by varying the magnetic field through the circuit is known 
as electromagnetic induction. Faraday's rule of induction states that the 
electromotive force (EMF) generated in a circuit is proportional to the rate of 
change of the magnetic flux through the circuit. Electromagnetic induction is 
the fundamental operating principle of transformers, inductors, electrical 
motors, generators, and solenoids. It is also known as Faraday’s Induction.” 
Faraday’s law of induction explains how an emf is generated in a conductive 
loop when the magnetic flux through the surface enclosed by the loop is 
varied in time. 


EMF = —d®/dt 
Eq (6) 


Where EMF is the electromotive force, ® is the magnetic flux, and dt is the 
time interval over which the magnetic flux changes. It also works in reverse 
and magnetic flux can be generated by a conductive loop which is the work- 
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ing principle of solenoids. The emf generated will try to resist the varying 
magnetic field and hence this principle is used for damping applications. ” 


2.4.2 Pyroelectricity 


Polarization is the separation of positive and negative charges in a material, 
which produces an electric dipole moment. When the material is polarized, it 
acquires a net electric charge, with one end of the material positive and the 
other end negative. Pyroelectric effect is observed in polarized materials and 
can be explained as the change in polarization due to the temperature 
change.”! This change in polarization results from a reorganization of electric 
charges in the material, produced by a temperature-induced shift in crystal 
symmetry. This shift in polarization can be detected as a shift in the electric 
charge on the material's surface. There are two main effects: 


e Direct effect- the temperature induced shift in charge distribution. 


e Indirect effect — The change in temperature causes temperature in- 
duced strain in the material and this results in strain induced polari- 
zation of the material. 


Materials that exhibit the pyroelectric effect are called pyroelectric materials 
and the change in electric charge per unit area of surface for a temperature 
change (AT) of 1°C is known as the pyroelectric coefficient. Polarization 
change of the material is given by: 


AP = p.AT 


Eq(7) 


Where AP denotes polarization and AT denotes temperature, while p is the 
pyroelectric coefficient. Eq ( 7 ) explains the linear connection between 
polarization and temperature in a pyroelectric substance and can be used to 
compute polarization at a given temperature or to calculate temperature 
change based on observed polarization change. 
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Figure 9 Schematic of pyroelectric effect. Adapted from Whatmore et al. 7 


The change in surface charge is collected by using an electrode coating on the 


surface and this develops a potential whenever the material experiences a AT. 
72 


AQ = pxAx AT 


Eq(8) 


Where, Q is the electric charge generated, p is the pyroelectric coefficient, 
which is a material-specific constant that relates the change in electric polari- 
zation to the change in temperature, AT is the change in temperature and A is 
the surface area of the material. According to Eq ( 8 ), the amount of the 
electric charge produced by the pyroelectric effect is proportional to the 
pyroelectric coefficient and the temperature change. Figure 9 explains the 
simple working principle of a pyroelectric effect. 


2.5 Heat transfer 


Heat transfer is the term used to broadly describe the various processes by 
which thermal energy is transferred either within the material or from one 
material to another. There are three types of heat transfer mechanisms or 
processes, conduction, convection, and radiation. Heat transfer by conduction 
is the result of mutual interaction of vibrating atoms without changing their 
mean position. This type of heat transfer mechanism is only observed in 
solids. Contact heat transfer is a special case of heat transfer by conduction in 
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which one material is in contact with another (without bonding and some- 
times intermittent). Contact heat transfer would be further explained in detail 
in Section 2.5.6. 


Heat transfer by convection involves the movement of atoms or molecules 
(also aggregates of molecules) carrying the heat from one position to another. 
Such heat transfer is observed in liquids and gases. Convective heat transfer 
can also occur at the boundary between solids and liquids/gases. Radiative 
heat transfer is a heat transfer mechanism through space as photons of elec- 
tromagnetic radiations of wavelengths of more than 10000 A. Unlike con- 
duction and convection which require a medium for transfer of heat from one 
position to another, the radiative heat transfer can occur without any material 
medium present and no physical contact is required between the two materi- 
als (or objects) between which the heat transfer is occurring.” 


2.5.1 Conductive heat transfer 


Conductive heat transfer, also known as thermal conduction, is the movement 
of thermal energy between objects in direct contact with each other or be- 
tween various temperature zones within a single object. Heat is transmitted 
through conduction in a solid substance by the impact of vibrating atoms and 
free electrons with adjacent atoms and electrons. Thermal energy is moved 
from the hotter to the colder regions of the item until the temperature is 
uniform throughout. Fourier's law governs the rate of heat transfer by con- 
duction, which says that the rate of heat transfer is proportional to the tem- 
perature difference in the direction of heat movement and thermal conductivi- 
ty of the substance, while inversely proportional to the distance. 
Mathematically, Fourier's law can be written as:”° 


q = —kA(dT /dx) 


Eq(9) 
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Where, q is the heat flow, k is the thermal conductivity of the material, A is 
the cross-sectional area through which heat flows, dT /dx is the temperature 
gradient in the direction of heat flow. 


2.5.2 Convective heat transfer 


Convection is the heat transmission between a solid surface and a moving 
fluid (liquid or vapor). It happens when the temperature of the solid surface 
and the fluid differs, causing the fluid to flow and transfer heat away from the 
surface. Convection heat transfer entails both the transfer of energy through 
heat transmission between the solid surface and the neighboring fluid layer 
and the transfer of energy through fluid motion. The rate of heat transmission 
by convection is determined by fluid characteristics such as velocity, temper- 
ature, and viscosity, as well as solid surface geometry and flow circumstanc- 
es. Convective heat transfer Q is given by:”° 


Q=h*A * (T,; — T;) 
Eq(10) 


Where, Q is the heat flow, h is the convective heat transfer coefficient, A is 
the surface area of the solid in contact with the fluid, T, is the temperature of 
the solid surface, Tris the temperature of the fluid. 


2.5.3 Thermal capacity 


Thermal capacity of a material is the amount of heat that needs to be trans- 
ferred to the material in order to increase its temperature by 1 °C. The specif- 
ic heat capacity is the amount of heat that is required to increase the tempera- 
ture of a unit mass of the material by 1 °C. Thermal capacity per unit volume 
can be calculated by multiplying specific thermal capacity by the material 
density. Thermal heat capacity of a sample of material is denoted by C, and 
is given by:’° 
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A: AQ 
P = AT 


Eq (11) 


Where the amount of heat that must be added to the material sample (of mass 
M) in order to increase its temperature by AT is denoted by AQ. Based on the 
material sample’s initial temperature (T) and the pressure (P) that is being 
applied to it, the value of C, often varies significantly. 


2.5.4 Thermal conductivity 


The thermal conductivity of a material defines the ability of a material to 
conduct heat. It is the ratio of heat flow through the unit area of a material to 
the temperature difference across its two surfaces. This relation is given by 
the equation Eq ( /2 ) which is derived from Fourier’s Law given in Eq (9). 
It is a material property, and the negative sign indicates that the flow of heat 
is from the high-temperature to the low-temperature side.’ 


Eq (12) 


Where Q is the heat flow through the material, « is the thermal conductivity 
of the material, and VT is the temperature gradient between the two surfaces 
of the material. Under steady-state conditions, the value of the temperature 
gradient remains constant, and this results in a constant heat flow through the 
material. But in the case of changing temperature gradient in a general 
transient state, Eq ( 73 ) describes the heat transfer through the material. 7° 


ƏT ƏT T 
OT = u 
E 32) = om 


Eq(13) 


25 


2 Basic principles 


The transient temperature distribution within a material over time when heat 
conduction is not in a steady state is given by Eq ( /3 ), and this depends on 
the density of the material (p), the thermal capacity of the material (Cp), the 
spatial gradient and the thermal conductivity of the material (V(«VT)). The 
heat flow in the material and VT experienced by the material will follow Eq ( 
13 ) whenever the material is exposed to a disturbance in the steady state 
condition either by intermittent contact with a heat source or sink, or by 
internal heat generation or a change in boundary conditions until a steady 
state temperature distribution can be reached again. 


2.5.5 Thermal resistance 


Thermal resistance describes how effectively the material or component 
opposes against the transfer of heat through it. Materials with very high 
thermal resistance per unit length are called insulators and materials with low 
thermal resistance per unit length are called thermal conductors. Thermal 
resistance can be calculated by Eq ( 14). ” 


l 
KA 
Eq (14) 


Where | is the thickness of the material,« is the thermal conductivity of the 
material and A is the area of the material side (cross-sectional area). This 
shows that thermal resistance is inversely proportional to thermal conductivi- 
ty which is a material property. Since the thermal resistance also has dimen- 
sional dependence like cross-sectional area and thickness, it is possible to 
tune the thermal resistance by tuning the dimensions. 


2.5.6 Contact resistance 


When two solid surfaces are in contact, the actual contact occurs only in 
discrete spots of micro areas and the rest of the space is either filled with 
vacuum or with some other media (air, water, oil etc.) due to the roughness 
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that is present in the solid surfaces. This media would have a thermal conduc- 
tivity that may be largely different from the thermal conductivity of the solids 
in contact and hence would influence the heat transfer through such a contact. 
The limitation in heat transfer due to the influence of the interface material 
(media or vacuum) is called thermal contact resistance. In classical heat 
transfer, conduction, convection, and radiation are the three mechanisms by 
which heat is transferred. In this case, it could be similar but, in most cases, 
only conduction needs to be considered since the convection and radiation at 
this interface is negligible compared to conduction.”® 


F 


Figure 10 Contact heat transfer between the solid surfaces in contact. Adapted from Wang et 
al.”8 
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Figure 10 shows the micro area contact between surfaces. The total heat 
transfer will be the sum of all the individual heat transfers through the micro- 
contact areas. The rate of contact heat transfer between two surfaces can be 
estimated using the following equation: 


Q=hxAxJAT 
Eq (15) 


Where, Q is the heat flow between the two surfaces, h is the overall heat 
transfer coefficient, which includes the contributions of both materials in 
contact, A is the actual contact area between the two surfaces, AT is the 
temperature difference between the two surfaces. The overall heat transfer 
coefficient (h) depends on the thermal conductivity and roughness of each 
material surface in contact, as well as the nature of the contact between the 
two surfaces (such as whether there is air or other fluids present in between). 


2.6 Lumped element modeling 


In electrical and electronic systems, resistance is considered to be a single 
physical component with linear relation V = IR, but this involves interactions 
of atoms and electrons with enormous complexity in geometry, random 
thermal motion, and collisions. This type of simplification by the concept of 
the “law of large numbers” in probability allows us to make a simplified 
model of a complex system by breaking it down into discrete blocks with a 
few inputs and outputs. Each block would describe the behavior of a discrete 
element at an average or aggregate level by employing algebraic, numerical 
of conceptual abstractions. If the limits of those abstractions and approxima- 
tions are understood and employed properly, these blocks can be cascaded 
together by using their inputs and outputs. A perfect lumped element model 
would be the simplest model that can be used to design or analyze a complex 
system but not simpler as this would cause errors to get cascaded and approx- 
imations would start to deviate the solution too much from the actual solu- 
tion. 
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2.6.1 Mechanical systems 


A mechanical system with elastic behavior, several damping factors, and 
attached mass can be simplified into a spring mass damper system. For 
example, a cantilever beam that bends upwards would have its tip moving in 
two dimensions. But this could be simplified into a one-dimensional spring 
damper system. With the addition of tip mass and cantilever mass into the 
equation, this would represent the oscillatory behavior of the system. Such 
systems could be used to describe complex dynamic systems by simplifying 
them into manageable spring mass damper systems. 


2.6.2 Thermal equivalent circuits 


When heat flow is considered to be analogous to electric current, the voltage 
will be analogous to temperature and the thermal resistance and thermal mass 
will be analogous to resistance and capacitance. Thus, a complex heat trans- 
fer problem can be broken down into blocks of electrically equivalent com- 
ponents and a circuit could be designed that would represent the heat flow 
and AT in a device. Such thermal equivalent circuits are used commonly to 
solve complex heat transfer problems. In the course of this thesis, MATLAB 
SIMULINK (SIMCSAPE) is used to design and evaluate thermal circuits to 
describe the device’s heat transfer characteristics. Commonly used approach- 
es include Foster and Cauer thermal networks.” 


2.7 Fabrication methods 
Several fabrication methods are used either directly at IMT or at partner 


universities to fabricate the films, cantilever, pick-up coil, and substrates of 
the device. The following section outlines the important fabrication methods. 
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2.7.1 Magnetron sputtering 


Magnetron sputtering is a physical vapor deposition method that falls under 
the category of vacuum coating techniques. This method is commonly used 
to deposit metals, alloys, and other compounds onto a wide variety of materi- 
als and up to a thickness of a millimeter. This method, although developed as 
a coating technique is used to deposit thin films on substrates that can be later 
released to have free-standing thin films. Magnetron sputtering is capable of 
producing thin coatings ranging in thickness from a few nanometers to 
several micrometers. However, parameters such as deposition rate, uniformi- 
ty of deposition, and target substance can restrict the thickness range.*° The 
high deposition rates, ease of sputtering any metal, high adhesion of films, 
and high purity film deposition make it very attractive over other deposition 
methods to fabricate thin films of materials. 


Figure 11 shows the schematic of the magnetron sputtering machine. When 
power is supplied to the magnetron, the target is supplied with a negative 
voltage (typically 300 V or higher). This negative voltage pulls positive ions 
to the target surface, generating a lot of kinetic energy at the same time. 
When a positive ion collides with atoms at the surface of a solid, an energy 
transfer occurs. Primary recoil atoms can be formed if the energy delivered to 
a lattice site is larger than the binding energy. These atoms can then collide 
with other atoms and disperse their energy via collision cascades. If the 
energy transferred in a direction normal to the surface is three times the 
surface binding energy, sputtering occurs (approximately equal to the heat of 
sublimation). One of the probable outcomes of ion bombardment of a 
surface is the sputtering of a target atom. Aside from sputtering, the emission 
of secondary electrons from the target surface is another significant activity. 
The glow discharge can be prolonged thanks to these secondary electrons. 
The target materials for the sputtering process are nearly limitless, ranging 
from pure metals that can be sprayed with DC power to semiconductors and 
isolators that require either RF power or pulsed DC. Deposition can be done 
with single or multi-elemental targets in either nonreactive (inert gas only) or 
reactive (inert and reactive gas) discharges.” 
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Figure 11 Schematic of Magnetron Sputtering. Adapted from Hassan et al.*! 


2.7.2 Melt spinning 


The melt spinning method is used to rapidly solidify metallic alloys by 
melting the metal and then fast cooling it by exposing it to a revolving wheel 
or cylinder chilled by a stream of gas or liquid. The molten metal is swiftly 
quenched when it comes into contact with the spinning wheel, resulting in a 
fine-grained solidified substance with unique characteristics. The method is 
frequently employed in the manufacturing of amorphous metallic alloys, also 
known as metallic glasses, which have mechanical, physical, and chemical 
characteristics that differ from those of typical crystalline metals. Unfortu- 
nately, only the side of the ribbon that touches the wheel directly is complete- 
ly amorphous. Amorphous metallic alloys may be manufactured in a number 
of forms, such as ribbons, wires, and powders, and are utilized in a variety of 
applications. The melt spinning process offers several advantages, including 
the capacity to generate huge amounts of material in a short period of time 
and materials with unique qualities that are not achievable with typical 
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casting and processing procedures. However, the process is difficult, necessi- 
tating precise control of the temperature, cooling rate, and other process 
factors in order to generate materials with the appropriate qualities. Figure 
12 illustrates the schematic of melt spinning. 
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Figure 12 Schematic of fabrication by melt spinning. 


2.7.3 3D printing 


Additive manufacturing (also known as 3D printing) is a parallel design and 
production process that incorporates materials, structures, and functions. 3D 
printing works by placing consecutive layers of material to create a three- 
dimensional structure. Fused deposition modeling, inkjet printing, selective 
laser sintering, electron beam melting, stereolithography, and digital light 
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processing are some of the 3D printing processes available. Typically, poly- 
mers are used as printing material either as filaments or in liquid resin form 
that is cured using UV light. However, selective laser sintering and electron 
beam melting are commonly used to print metal parts. 


Recently magnetic shape memory alloys have been 3D printed. Induction 
melting of high purity Ni, Mn, and Ga elements in an argon environment 
produces Ni-Mn-Ga polycrystalline ingots. A Retsch planetary mill is used to 
crush and ball grind Ni-Mn-Ga ingots for 3 hours. Sieving is used to obtain a 
powder with a particle size of less than the desired limit. These powders are 
used to 3D print Ni-Mn-Ga magnetic shape memory alloys with regulated 
densities ranging from 45 to 99 percent using binder jet printing. The 3D 
printed Ni-Mn-Ga alloy is then subsequently sintered to avoid oxidation, 
composition changes, and gradients. Further information regarding the 3D 
printing of magnetic shape memory films and their properties can be found in 
Mostafaei et al.*? 


2.7.4 Heat treatment 


Heat treatment is an effective method used to achieve the desired mechanical 
properties for metallic alloys. The phase transformations occurring within the 
material can be fine-tuned by controlling the temperature and time of heat 
treatment. Both the temperature and time need to be precisely controlled to 
achieve repeatable results for the material properties. Heat treatment can 
produce a variety of alterations in the microstructure of the SMA film, 
including recrystallization, grain growth, and phase transition. In the instance 
of a nickel-titanium (Ni-Ti) SMA, the heat treatment procedure can result in 
the creation of a crystalline phase known as R-phase, which improves the 
shape memory effect. The heat treatment procedure can also enhance the 
SMA film's mechanical characteristics, such as rigidity and ductility. The 
chemical composition of the magnetron sputter deposited films is determined 
by the sputtering power, and the thickness of the film is determined by the 
duration of sputtering. The material properties of the deposited films can be 
adjusted by heat treatment post deposition. As an example, in Ohtsuka et al. 
84 the martensite phase was changed from cubic to tetragonal after heat 
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treatment which is confirmed using XRD measurements. The deposited films 
had columnar grain structure as deposited, which increased in its width upon 
heat treatment. The heat treatment would slightly change the composition if 
any of the contents are evaporating at high temperatures.** In some cases heat 
treatment process results in an increase in magnetization of the Ni-Mn-Ga 
films and sharper transition with higher magnetization change per unit 
change in temperature. 85 


2.7.5 Optical lithography 


Optical lithography is a method used to pattern and copy a desired image 
onto a substrate, such as a silicon wafer. A chemical compound called photo- 
resist, which is light-sensitive, is coated on a surface or substrate. Usually, 
photoresist is coated on the substrate surface using spin coating process. 
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Figure 13 Schematic of fabrication using a lithography process. 


The photoresist will be exposed either by using a mask in case of UV lithog- 
raphy or directly writing on the resist by using a UV laser in case of mul- 
tiphoton lithography. Both of the methods will result in a 2D image of the 
structure being transferred to the resist layer. This layer is then developed 
using a suitable developing solution removing either the exposed areas in 
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case of positive resist or the unexposed areas in case of negative resist. The 
selection and application of a positive or negative resist depends on by the 
intended design and the processing specifications. Positive resists are usually 
used for simpler patterns with bigger features, whereas negative resists are 
used for more complicated patterns with smaller features. Figure 13 shows 
the process flow diagram for optical lithography by which a pattern can be 
impressed on the photoresist. 


2.7.5.1 UV lithography 


UV lithography uses electromagnetic waves of ultraviolet wavelengths to 
produce a pattern on the spin-coated photoresist. A photomask is used to 
generate the pattern on the photoresist by blocking some areas and exposing 
other areas selectively as shown in Figure 14. Light diffraction restricts 
optical system resolution, which limits the smallest feature size that can be 
resolved. The Abbe diffraction limit determines the diffraction limit, which is 
affected by the wavelength of light used for the procedure.®° 


UV lithography is a common procedure used in the semiconductor industry 
to create integrated circuits (ICs). It's also utilized to make microelectrome- 
chanical systems (MEMS), microfluidic devices, and other micro- and nano- 
scale structures. It is a popular choice for many applications due to its high 
resolution, rapid throughput, and outstanding repeatability. 
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Figure 14 Schematic using UV lithography. 
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2.7.5.2 Multiphoton lithography 
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Figure 15 Schematic of multiphoton lithography or direct laser writing. 


Figure 15 illustrates the schematic of the direct laser lithography. Multipho- 
ton lithography also known as direct laser writing is a process of exposing 
negative-tone or positive-tone photoresists with laser light of a certain wave- 
length. Structuring is performed similarly to conventional photolithography 
methods. However, the avoidance of reticles or photomask is the key distinc- 
tion. In place of that, two-photon absorption is used to cause a significant 
alteration in the solubility of the resist for the right developers. This tech- 
nique makes use of a multi-photon absorption procedure in a substance that is 
transparent at the laser's pattern-creating wavelength. A chemical change 
(often polymerization) takes place at the laser's focus point and may be 
managed to produce any desired three-dimensional periodic or non-periodic 
pattern by scanning and modulating the laser appropriately.°” This technique 
has been used to quickly prototype intricate structures. 
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2.7.6 Wetetching 


Wet etching is a method for removing materials off substrates utilizing liquid 
etchants or chemicals. Photoresist masks on the wafer specify the precise 
patterns that are achieved by lithography. Liquid chemicals are used to etch 
away areas that are not shielded by this mask. In cases where the etching 
times are longer, a hard mask with chromium and gold deposition on top of 
the material is used. Multiple chemical reactions that take place during a wet 
etching process use up the initial reactants and generate new ones. Three 
fundamental steps may be used to summarize the wet etch process. (1) The 
liquid etchant is diffused into the structure to be eliminated. (2) The interac- 
tion between the substance being etched away and the liquid etchant. Typi- 
cally, a reduction-oxidation (redox) process takes place. In this process, the 
substance is first oxidized, and subsequently, the oxidized material is dis- 
solved. (3) Diffusion of reaction byproducts away from the surface that 
underwent the reaction.®® 


2.7.6.1 Isotropic and anisotropic etching 


A material will either be removed by isotropic etching (uniformly in all 
directions) or by anisotropic etching when it is attacked by a liquid or vapor 
etchant (uniformity in the vertical direction). 


Photoresist 
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Figure 16 Schematic of various wet etching processes (a) anisotropic etching, (b) partially 
anisotropic, (c) isotropic etching. 
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The schematic illustration of various etching processes is given in Figure 16 
Wet etching often removes material at rates that are quicker than many dry 
etching methods, and it is simple to alter these rates by adjusting the tempera- 
ture or the concentration of active species. 


2.7.7 Laser cutting 


By use of so-called "ablation processes," the energy of a concentrated laser 
beam may be utilized to remove material. A high-power laser's output is 
typically directed through optics in a laser cutter. To focus the laser beam on 
the material, laser optics and CNC (computer numerical control) are utilized. 
Using a motion control system, a commercial laser used for material cutting 
can follow a CNC or G-code that specifies the pattern to be etched onto the 
substrate. When a concentrated laser beam is pointed at a substance, it either 
melts, burns, vaporizes or is blasted away by a jet of gas, leaving an edge 
with a high-quality surface finish. 


2.7.7.1 Ultra-short pulse laser ablation 


Because of the speed at which the energy is supplied by ultrashort laser 
pulses, ablation occurs before heat conduction to the bulk material. In order 
to treat materials without thermally impacting the surrounding bulk material, 
ultrashort laser pulses are utilized. With significantly less heat-related nega- 
tive thermal effect, this may be utilized to fabricate structures more accurate- 
ly. Treating materials with ultrashort pulses has the additional benefit of 
processing almost any substance. This results from the extraordinarily high 
energy densities that permit multi-photon absorption. It is also possible to 
process materials that are transparent to the laser's wavelength as explained in 
the section 2.7.5.2. 


2.7.8 Semi-automatic coil winder 
The pick-up coil for TMGs is made in this work utilizing an in-house manu- 


factured semi-automatic coil winder. By spinning the core and counting and 
showing the number of turns in the pick-up coil, the semi-automatic winder 
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enables for reliable and reproducible pick-up coil construction. This guaran- 
tees that the pick-up coil's internal resistance remains within the tolerance 
limit, resulting in more reliable and precise readings. The gadget is outfitted 
with a programmed Arduino Nano microprocessor that offers fine control 
over the winding process. In addition, the gadget has a motor with an inte- 
grated magnetic encoder, allowing for precise and exact measurements of the 
number of turns in the pick-up coil. Figure 17 depicts the setup as 3D render- 
ing from CAD software. Furthermore, the combination of the Arduino Nano 
and motor with magnetic encoder allows for further expansion of the winder 
to be automatic completely getting rid of human intervention requirements. 


Figure 17 3D model of the Semi-Automatic coil winder built in-house for fabricating pick-up 
coils. 


2.8 Thermomagnetic materials 
Magnetic shape memory alloys like Heusler alloys and ferromagnetic materi- 


als like gadolintum (Gd) and lanthanum iron silicon (La-Fe-Si) alloys that 
undergo ferromagnetic ordering at a particular temperature can be commonly 
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categorized as thermomagnetic materials (TM films). These materials exhibit 
a change in their magnetization in response to changes in magnetic field and 
temperature. Both ofthese can be used for achieving self-actuation in TMGs. 
Figure 18 a and b, show the magnetization change in the material Ni-Mn-Ga 
and La-Fe-Si-Mn-H with respect to change in both magnetic field and tem- 
perature. 
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Figure 18 (a) Magnetization vs temperature graph for a Ni-Mn-Ga alloy. (b) Magnetization vs 
temperature graph for a La-Fe-Si-Mn-H alloy (c) Magnetization vs magnetic 
field for Ni-Mn-Ga alloy at room temperature (d) Magnetization vs magnetic 
field for La-Fe-Si-Mn-H at room temperature. 


In some alloy compositions, the material might still have remanent magneti- 
zation even after the transformation of the material at T. as shown in Figure 
18b. The dependence of magnetization value on the magnetic field is also 
dependent on the alloy composition. Figure 18c shows that the given Ni-Mn- 


40 


2.9 Operation principle of TMG self-actuation 


Ga alloy composition has a very low saturation magnetic field at room 
temperature, thus showing less dependence of the magnetization value on the 
magnetic field. The given composition of La-Fe-Si-Mn-H alloy has a higher 
saturation field as shown in Figure 18d, as the value of magnetization is 
increasing with the increase in the value of the magnetic field. It is also 
important to note that the ferromagnetic to paramagnetic transformation in 
TM films is a continuous second-order transformation. But there are also 
materials like metamagnetic shape memory alloys (MSMAs) which present 
paramagnetic or antiferromagnetic martensite to ferromagnetic austenite 
transformation through first-order transformation. In first-order transfor- 
mations, the material would abruptly transform into a phase often accompa- 
nied by Hysteresis. 


2.9 Operation principle of TMG self-actuation 


Figure 19 illustrates the design layout and operating principle of the TMG. A 
cantilever with the free end mounted on a ceramic substrate. Cantilevers 
often have a greater natural frequency than bridge beams, allowing them to 
operate at higher frequencies and power densities.° Cantilevers can be 
created utilizing well-established and reasonably affordable microfabrication 
processes. A cantilever's resonance frequency may be changed by modifying 
its size, making it easier to optimize performance for a specific application. 
Brass is chosen for the cantilever material due to its mechanical and thermal 
properties. The mechanical and magnetic parts allow for each to be tuned 
separately as shown by Gueltig et al.. At the free end of the cantilever, a 
thin film of thermomagnetic material (TM film) is placed by using a bonding 
layer. A pick-up coil is placed underneath the beam which is used to convert 
mechanical motion into electrical energy using electromagnetic induction. 
The whole cantilever assembly is placed beneath a tilted permanent magnet 
to achieve maximum contact once the cantilever is bent. The permanent 
magnet is made of Samarium Cobalt (SmCo) which can tolerate high operat- 
ing temperatures up to 350 °C and serves as the heat source for the system 
with a temperature higher than the Te of the TM film. When the temperature 
of the TM film is below its Tc, the material is ferromagnetic and hence has 
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large magnetization when placed under the magnetic field. The magnetic 
field gradient and magnetization results in large magnetic attraction force 
between the TM film and the magnet, which pulls the free end of the cantile- 
ver tip to bend up to the magnet surface. 


a) b) 
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Figure 19 Schematic of the thermomagnetic generator self-actuation (a) heating at contact with 
the permanent magnet as the heat source, (b) cooling during the self-actuation. 
Reprinted from Joseph et al.”! 


This increases the temperature of the TM film above its T. and the material 
becomes paramagnetic causing a sudden drop in magnetization of the materi- 
al, resulting in an abrupt change in magnetic force. When the restoring force 
is larger than the magnetic force, the cantilever bends back to its original 
position with an overshoot due to the inertia of the tip mass. Once the TM 
film loses contact with the magnet, it starts to cool down and this results in 
increasing magnetization and hence magnetic force. When the magnetic force 
is larger than the restoring force, the cantilever is bent back to the magnet 
surface. This cycle results in a repeating self-actuated movement of the 
cantilever tip towards the magnet surface and away from it. 


2.9.1 Thermomagnetic cycle for self-actuation 


The thermomagnetic cycle (TM cycle) involved in the operation of the TMG 
under self-actuation, could be broken down into four stages. Figure 20 
illustrates all four stages of the TM cycle followed by the TMG under self- 
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actuation. In an ideal TM cycle, Stage I is an increase in the magnetic field 
which increases the magnetization of the TM film, this is achieved when the 
cantilever tip moves towards the magnet. Stage II denotes the decrease in 
magnetization in the TM film, and this is the result of increasing temperature 
at contact with the magnet surface. Stage III describes the decrease in the 
magnetic field, this is achieved when the cantilever tip moves away from the 
magnet and the TM film then experiences a decrease in the magnetic field 
and therefore the magnetization also reduces. Stage IV indicates the TM film 
cooling down, causing an increase in magnetization within the material. The 
actual TM cycle under which the presented TMG operates does not have 
defined Stages III-IV-I as the TM film starts to cool down as soon as it leaves 
the surface of the magnet. Hence there is no isothermal actuation in Stages I 
and III. The actual cycle is represented by the shaded area inside. This results 
in a more gradual increase of magnetization over the cooling cycle which is 
also combined with the change in magnetization due to the reduction and 
increase in magnetic field during the cooling cycle. 
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Figure 20 TM cycle showing stages I-IV in a TMG. Adapted from Joseph et al.” 
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2.9.2 Resonant self-actuation 


Self-actuation is the key working principle utilized in this work. There is heat 
intake at contact and heat dissipation during the oscillation away from the 
magnet surface. This results in tip displacement as shown in Figure 21a. The 
short duration in which the displacement 
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Figure 21 Self-actuation without resonance, (a) displacement, (b) current, (c) power, (d) plotting 
displacement vs current to check cyclic behavior. 


is constant without oscillations denotes the TM film in contact with the 
magnet for heat intake. The cooling of TM film is achieved by the damped 
oscillations that eventually cool the TM film by conduction and convection, 
enough for it to be attracted back to the magnet surface. Figure 21b shows 
the current generated in the pick-up coil for cantilever tip displacement in 
Figure 21a. It can be observed that the current value depends on the stroke 
and since the stroke is decreasing in each oscillation until it goes back to the 
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magnet, the current decreases accordingly. Figure 21c shows the power 
generated for the oscillations and that there is still potential to improve the 
power output by increasing stroke of these additional oscillations or by 
increasing the frequency of operation. 
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Figure 22 (a) displacement measures (b) current measurements, (c) instantaneous power 
derived from current, (d) current vs displacement to check cyclic behavior. 


Since the power is largest when the stroke is largest, it would be ideal to have 
the cantilever making contact with the surface of the magnet in each oscilla- 
tion and not after having any additional damped oscillations. This can be 
achieved by optimizing device design, which includes: 


e fabrication technology (bonding layers and material selection) 
e electrical load resistance 


e optimization of cantilever geometry affecting footprint 
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e heat intake during contact between the TM film and heat source (per- 
manent magnet) 


e heat dissipation by conduction and convection 


Both heat intake and heat dissipation are dependent on the dynamics of 
cantilever motion. The heat taken in during the contact should be just enough 
to cause a drop in magnetic force, which results in restoring forces taking 
over. This would enable fast operation at high frequency by taking only a 
small amount of heat in each cycle. Heat dissipation also needs to be opti- 
mized by blocking most of the heat at contact so that the temperature increas- 
es quickly and then dissipate the heat into the cantilever over the cooling 
cycle. This can be achieved by using a non-conductive bonding layer with 
large thermal resistance between the TM film and cantilever. Once the heat 
intake and heat dissipation are brought in balance with the oscillation of the 
TMG, current and power output becomes periodic with identical cycles.?®”! 
Resonant self-actuation was first presented in Gueltig et al.” Figure 22 
shows the displacement, current and instantaneous power output of the TMG 
in resonant self-actuation. The major challenge is to have TMGs with heat 
transfer optimized to have heat intake in a short contact time with the magnet 
and to have heat dissipation completed over the rest of the cycle to have 
regular and periodic AT. If AT is not regular the device will experience heat 
buildup and would either have an additional oscillation or a complete stop. 
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Initial concepts were first proposed in the nineteenth century and numerous 
theoretical and experimental investigations on the use of ferromagnetic 
materials, such as cobalt (Co), iron (Fe), and gadolinium (Gd), have been 
conducted since then. The interest in TM generation has recently regained 
momentum.”>34 Although the theoretical relative efficiencies of ferromag- 
netic materials operating at or near their Tc can reach up to 55 % of Carnot 
efficiency, the actual values achieved in experiments are much lower.” 
Several concepts have been introduced over the years and all of them could 
be categorized into two different types, the first one being direct conversion 
in which the change in magnetization is directly converted into electrical 
energy. The second type of device uses an intermediate stage of mechanical 
self-actuation achieved by periodic heating of TM material and then using the 
mechanical actuation for energy conversion. Direct conversion is usually 
called active TM conversion and indirect conversion is called passive TM 
conversion.” 


3.1 Early devices 


Early TMGs had iron as the ferromagnetic material which needs 700 °C for 
demagnetization and had very complicated designs. This limitation stagnated 
further development as such high temperatures had unfavorable effects on the 
permanent magnet used in these devices.” 


3.1.1 Pyromagnetic motor 


Thomas Edison proposed the idea of a pyromagnetic motor in 1888, which 
produced heat by burning coal or wood.”® Design of the pyromagnetic motor 
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by Thomas Edison is shown in Figure 23. The device is built with an intersti- 
tial armature with interstices or tubes extending longitudinally through it 
using a thin iron sheet. The armature is positioned between the magnetic 
poles of permanent or electromagnets and is installed vertically using a shaft. 
Underneath the armature is a furnace with two outputs covering the lower 
ends the armature tubes on the opposing sides. These two groups of tubes are 
subjected to hot air blowing from the furnace outlets, while the remaining 
group of tubes situated in between the hot tubes is subjected to cold air. 
There are cold tubes in the center and hot tubes on the opposing sides. The 
configuration did not restrict the armature's ability to rotate. The magnetiza- 
tion of the armature varied across various portions because of the temperature 
variation, producing an unbalanced force. As long as the armature tubes are 
periodically heated and cooled, this resulted in a continuous rotation of the 
armature. 


Figure 23 Pyromagnetic motor from US patent 380100 by Thomas Edison” 


3.1.2 Thermomagnetic motor 


Nikola Tesla patented a thermomagnetic motor under US Patent 396121 in 
the year 1889.?’ Figure 24 shows some of the early concepts of Nikola Tesla. 
By using the combined effects of heat, magnetism, and stiffness of a spring, 
the presented devices are able to generate mechanical energy through a 
reciprocating or rotational action. In Tesla's thermomagnetic motor, an 
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armature is in contact with an electromagnet or permanent magnet that 
produced a magnetic field. A hinged arm with the armature attached to it 
could freely rotate. The armature is separated from the permanent/electro- 
magnet when the temperature of the armature is increased over its Tc as a 
result of the heating burner. After a while, the armature is cooled below the 
Te, at which point the magnetic force increased more than the spring force 
and returned the armature to its initial position. As long as the heating and 
cooling cycles persisted, this process is repeated. However, because of their 
poor performance, these TMGs were never deployed in either domestic or 
industrial applications. 


Figure 24 Designs of the thermomagnetic motor by Nikola Tesla in US patent 396121” 


3.1.3 Pyromagneto-Electric Generator 


A patent for a pyromagnetic electric generator was first introduced by Nikola 
Tesla in 1890.” A device that directly produces electric current without 
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mechanical motion is described in this patent. An armature core made of a 
few iron tubes is used to bridge a permanent magnet, as seen in Figure 25. 
Around the armature core are two conductors that would serve as the source 
of the electric current. A furnace placed beneath the center of the armature 
core will heat it, while steam originating from a boiler can cool it. Alternately 
cooling and heating the armature turned on and off the magnetic circuit, 
changing the magnetic field, which eventually led to an electrical signal in 
the leads. 


Figure 25 Pyromagneto-Electric Generator patented by Nikola Tesla in US patent 428057°® 


3.1.4 Pyromagnetic Generator 


Another type of pyromagnetic generator with bundles of thin iron tubes 
arranged circumferentially inside two iron rings was invented by Thomas 
Edison in 1892.” Figure 26 provides an illustration of this gadget. These 
iron tube configurations are magnetized using strong permanent magnets or 
electromagnets so that the rings formed the two poles of a long magnet. Each 
bundle of tubes had a wrap around it. Under the lower iron ring is a furnace, 
and some of the tubes' lower ends are shielded with shields in the shape of 
half disks. The open tubes heated up after receiving the hot furnace air. 
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However, because the barrier shielded the tubes from the furnace's heat, they 
are able to cool. The tubes are always thermally unstable because the shield is 
gently rotated around its axis by a motor. Iron tubes are gradually cooled and 
heated, changing their magnetization in the process, which produced current 
in the coil turns. This device along with the previous device designs ex- 
plained until now is not used in any form in any commercial or domestic 
applications due to the very low efficiency, high temperature requirement, 
and lack of permanent magnet materials that had high Te to sustain perfor- 
mance at such high temperatures. 


Figure 26 Pyromagnetic generator design patented by Thomas Edison in US patent 476983 
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3.2 Recent thermomagnetic energy generators 


Due to the extremely low efficiencies that the available materials could 
achieve, research on TMGs has advanced very gradually over the past centu- 
ry. Brillouin and Iskenderian reexamined the concept of TM power genera- 
tion for the first time in 1948 and reviewed the design of a static TMG for the 
direct conversion of heat into energy.!° Gadolinium is an intriguing material 
and is used as an active material for subsequent research.'°'! Gadolinium, 
which was then more widely available, is primarily considered for experi- 
mental work since it has a Te close to room temperature.!” Kirol and Mills 
put out various theories based on the fluidic heating and cooling of TM 
film.!® Solomon carried out more computations and optimizations in 1988 
and created a general model." Devices described from section 3.2.1 to 
section 3.2.3 are examples of direct conversion or active TMGs. Devices 
presented from section 3.2.4 to section 3.2.7 are examples of passive TMGs. 


3.2.1 Direct energy conversion using first-order 
transformation 


Srivastava et. al. presents a TMG that directly converts the change in temper- 
ature to electricity using the change in magnetization of a shunt material. 
Figure 27 shows the image of the prototype demonstrator and the schematic 
of the device. A pick-up coil surrounds the sample, which is secured (to 
restrict movement) close to a pole of a permanent bar magnet. An external 
heat gun provides the heat. The fundamental concept is that through a trans- 
formation during heating, the specimen's average magnetization M increases, 
producing a significant change in magnetization for a unit temperature 
change (AM/AT). The proximity of the permanent magnet prevents the 
demagnetization of the material. A field is generated during transformation 
and the pick-up coil is designed with the goal of providing the maximum 
component of the field parallel to the wire, which will cause a current to flow 
in accordance with Faraday's law. On the reverse phase change after cooling, 
a potential difference of opposite polarity is produced across the pick-up coil. 
Although the promised high-power density competes with cutting-edge 
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thermoelectric energy harvesters, the experimental efficiencies are quite low. 
In reference,!”’ a theoretical model of the efficiency is provided. The hystere- 
sis of the phase transformation must be near zero for energy harvesting at low 
temperature differences. The demonstrator generates around 50 uW of power 
and this corresponds to a power output per footprint value of 28 W/cm? in 
ideal conditions where the temperature is externally oscillating between 
approximately 25 °C and 40 °C. 


Figure 27 Adapted image from,'® (a) photograph of the demonstrator, (b) schematic of the 
demonstrator. Reprinted with permission from John Wiley and Sons. 


NiasCosMnaoSnio, a quaternary Heusler alloy, is presented by Srivastava et 
al.,'%® demonstrating highly sharp transitions and hysteresis of only 6 °C. 
Along with narrow hysteresis, good martensite, austenite phase compatibility 
and good first-order phase transition reversibility ensure low functional 
fatigue. 


3.2.2 Static thermomagnetic generator 


The device shown in Figure 28 from Christiaanse et al.'” consists of two 
active regenerators made from (Mn,Fe)2(P,As) materials suspended in a 
magnetic circuit, and a permanent magnet acting as a field source. The 
regenerators are made up of 48 disks of (Mn,Fe)2(P,As) material with four 
different Curie temperatures, and microchannels are laser cut into each plate. 
By stacking these plates in series, two regenerators with a gradient tempera- 
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ture span are built. As a result, magnetic stray fields are prevented since the 
permanent magnet's generated magnetic flux is always directed through one 
of the TM disk. To determine the highest power output at the optimal fre- 
quency, the material is cycled at various frequencies. 


The device uses the regenerative Ericsson cycle, which consists of two 
isothermal curves and two isofield curves, to achieve higher efficiencies. The 
researchers tested the device by cycling the materials at different speeds to 
find the maximum power output at the optimum frequency. They found that a 
more tapered design and different materials would be critical to a better 
design. The device uses the regenerative Ericsson cycle, which consists of 
two isothermal curves and two isofield curves, to achieve higher efficiencies. 
The researchers tested the device by cycling the materials at different speeds 
to find the maximum power output at the optimum frequency. They found 
that a more tapered design and different materials would be critical to a better 
design. The tapered design ensures a greater external field and a larger 
external field change across the material. 


a) Paramagnetic MnFe(P,As) Plates b) 


Ferromagnetic MnFe(P,As) Plates 


Figure 28 Adapted image from Christiaanse et al.'” 


proposed device in Christiaanse et al..'°° Paramagnetic state magnetocaloric 
material in red, and Ferromagnetic state magnetocaloric material in blue. b) 
Shows the MnFe(P,As) plates with microchannels. Reprinted with permission 
from Springer Nature. 


a) Three-dimensional representation of the 
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Additionally, the materials for devices that use water as a heat-transfer fluid 
should have a lower latent heat while maintaining their magnetic properties to 
increase cycling speed and improve power output. This will increase the 
speed of the process cycle and boost the device's power output. At a frequen- 
cy of 0.17 Hz, the power production is at its highest. The heat transfer to the 
regenerator is decreased when frequency is increased further, preventing the 
material from cycling through the magnetic phase transition, which results in 
the power reduction. The fundamental cause of this relatively poor cycling 
behavior is the laser treatment's degradation of the magnetic characteristics 
while the microchannels were cut. In essence, the walls of the water channels 
are transformed by heat from the cuts into a material that is not magnetic, but 
simply has thermal resistance and heat capacity. 


3.2.3 TMG with pretzel-like magnetic topology 


The magnetization should, in an ideal scenario, abruptly change AM around 
Te of the TM film to consider the TM material as a thermal switch. Waske et. 
al. has shown a pretzel-like topology for the magnetic circuit to avoid stray 
fields and to double the flux change in the circuit.!!° This is achieved by the 
sign reversal of the magnetic flux within the connecting yokes when the 
temperature of the TM materials is changed alternatively. This design is an 
extension of the thermomagnetic generator presented by Nikola Tesla.’ 


Figure 29 Thermomagnetic generator with a pretzel-like magnetic topology introduced in 
Waske et al.,''° (b) image of the TMG prototype, (b) different topology of 
magnetic circuits discussed. Reprinted with permission from Springer Nature. 
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The magnetic circuit depicted in Figure 29 includes the TM material as a 
component. Permanent magnets are the source of the magnetic flux in this 
circuit, and two soft magnetic yokes are used to direct the flux. The TM 
materials of this magnetic circuit act as a thermal switch for the magnetic 
flux. The developed TMG with the new topology allows the device to 
achieve a high power output of 1.25 mW from a very small AT of 30 °C 
resulting in a power per footprint of 32 uW/cm?. However, the pump and 
other fittings to circulate the fluids for heat transfer are not taken into ac- 
count. By making some design changes to improve the cycling frequency of 
the device, it is claimed by authors to be possible to increase the power 
output and efficiency of conversion. 


3.2.4 Miniature ferromagnetic thermal energy 
harvester 


A miniature energy-harvesting prototype developed by UCLA researchers is 
based on the change of Gadolinium’s magnetization with temperature 
change.'!! As seen in Figure 30, a piezoelectric spring mechanism carries the 
ferromagnetic material between a hot and cold reservoir. The findings include 
a description of the mechanical energy present in the system and an analytical 
method for calculating the potential electrical output. 


(b) Leave Spring Ambient/Cold Side 
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Soft Ferromagnet 


Hot Side t Hard Magnet 
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Figure 30 TMG using ferromagnetic Gadolinium as active material presented in, (a) actuation 
device, (b) cross-sectional view during cooling, (c) cross-sectional view during 
heating.''! Reprinted with permission from AIP publishing. 
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High-efficiency piezoelectric energy conversion is proposed by the authors 
capable of converting up to 3.6 mW/cm? from a temperature difference of 
50 °C. In the following research by Bulgrin et. al.!!?, the system design is 
theoretically examined in terms of power output and efficiency, and it is 
concluded that with optimizations, an average power output of more than 
250 uW is feasible, translating to a power per footprint of up to 62 uW/cm?. 
However, the experimental verification is still pending for the claimed large 
power per footprint. 


3.2.5 La-Fe-Si based thermomagnetic generator 


In 2019, Ahmim et. al. presented a TMG based on La-Fe-Si-H alloy.!!3 After 
two years, in 2021 an updated version of the device is presented.'!* The TMG 
shown in Figure 31 is constructed using an magneto caloric material (MCM) 
plate, a field source, a thermal energy source, a heat sink, and a disc spring. 
Two piezoelectric buzzers are used to create the disc spring. Through a 
miniature CuBe ring, they are inversely bonded to one another. NdFeB 
parallelepiped magnets are used to construct the Halbach structure, which 
serves as the field source. 


Magnets 


2 Piezoelectric 
<< MCM 


buzzers 


Heat sink 


Figure 31 Schematic of the TMG with La-Fe-Si as active material and piezo spring as energy 
transducing material.''* Reprinted with permission from Elsevier. 
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In comparison to the field produced by a single magnet, this combination is 
claimed to provide a larger force with a steeper change across small dis- 
placements. While the cold reservoir is constructed with a water-cooled heat 
exchanger, the temperature of the heat source is maintained through electrical 
resistance heating. The MCM has direct thermal contact with the heat sink 
through a hole in the piezoelectric buzzer. The presented device is capable of 
generating a maximum output of 4.2 uW from a small temperature difference 
of 35 °C. The resulting power per footprint is around 0.12 uW/cm?. The low 
value is due to the large size of the device even though the power density 
with respect to active material volume is around 240 uW/cm’. 


3.2.6 Hybrid thermomagnetic generator 


A hybrid device that combines TM and triboelectric effects to produce 
electricity is demonstrated by Rodrigues et al.!!5 In the presence of slight 
temperature gradients close to room temperature below 30 °C, the TM effect 
enables the periodic and sustained movement of the ferromagnetic material 
with the second-order transition. Using a triboelectric nanogenerator 
(TENG), this mechanical actuation is subsequently transformed into electrical 
energy. A TENG is generally made up of two materials having opposing 
electron affinities that are brought into contact with one another. Electrons 
move from one substance to the other when the materials are rubbed or 
forced together, resulting in charge separation and the production of an 
electrical potential. This energy may be captured and utilized to power small 
electrical gadgets or charge batteries.'!® The presented hybrid TM and triboe- 
lectric nanogenerator is depicted in Figure 32. This concept's application is 
shown in both the cold (15 to 37 °C) and hot (60 to 90 °C) sides of a wide 
range of operating temperatures. The electrical output produced by the hybrid 
TENG is further claimed to have 35 times more than that produced by a 
typical TM generator driven by magnetic induction. The device generates a 
large power output of 481 uW but with a power per footprint of 
5.47 uW/cm?. Owing to the large area of the device, the power per footprint 
is still low compared to other TMGs with magnetic induction.” 
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Figure 32 Hybrid thermomagnetic and triboelectric nanogenerator presented in Ahmim et al.!!> 


Reprinted with permission from John Wiley and Sons. Legend: TENG- triboe- 
lectric nanogenerator, PP tape- polypropylene tape, PTFE- polytetrafluoroeth- 
ylene. 


3.2.7 Pyroelectric conversion using thermomagnetic 
actuation 


A novel method of energy generation that uses a hybrid method of both TM 
actuation and pyroelectric energy transduction is presented in Lallart et al.!!7 
The authors of this study describe the development of a heat engine for 
small-scale energy harvesting in which a TM film transfers its heat to a 
pyroelectric element. It converts thermal energy into electrical energy, 
resulting in a hybrid method of direct energy conversion path as opposed to 
conventional electromechanical heat engines as shown in Figure 33. Addi- 
tionally, faster cooling of the TM film, which accounts for a higher thermal 
cycling frequency, is made possible by thermally separating the pyroelectric 
element from it. The power per footprint generated by the device is 
0.02 uW/cm?. 
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Figure 33 The schematic of the TMG with a pyroelectric element for energy conversion directly 
from heat to electricity.''’ Reprinted with permission from Applied Energy. 


These results show that hybridizing does yield some advantages but until 
now none of the presented results have generated more power than the TMG 
with magnetic induction. Several other devices have also been published in 
recent times which utilize the TM energy generation principle using various 
designs.''*!3 These devices have shown that TM energy generation is 
gaining huge research interest due to its high-power density and high effi- 
ciency in low temperature differences. 


3.3 Preceding work on TM energy generation 


This thesis work is built upon previous generation TMGs that are developed 
at KIT, 42439092,124-126 This section explains the previous work at KIT. 


3.3.1 Thermomagnetic actuations 


TM conversion principles can be employed for thermal actuation in addition 
to generating power from heat. Kohl et. al provides an overview of miniature 
magnetic SMA-based actuators.”” A micro scanner is presented that utilized 
magnetic shape memory alloy film in the form of a cantilever to achieve 
actuation based on the electrical input.**!?7 A mirror is mounted to the front 
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end of a double beam cantilever that is formed from free-standing Ni-Mn-Ga 
films and anchored to a substrate at the other end as illustrated in Figure 34. 
The double-beam cantilever may be heated by recurrent electric heating 
pulses since it is electrically connected. A magnet is positioned above the 
cantilever, deflecting it in its initial state. The temperature of the cantilever 
beam is controlled by the electrical pulses (heating) and convection is used to 
cool it. Thus, periodically heating the beam allows the microscanner to 
achieve oscillations with frequencies up to 200 Hz with extremely large 
deflections which allow this system to be adapted to portable mobile optical 
sensor systems. This actuation principle laid the basis for the subsequent 
development of energy harvesters with an oscillating cantilever design. 


Aa 


| Fsma 


Figure 34 Microscanner using Ni-Mn-Ga film as cantilever material and a micromirror attached 
to the tip.” Reprinted with permission from IOS press. 


3.3.2 TMG with Ni-Co-Mn-In 


TMGs using MSMA (Niso.4Co3.7Mn32.8In13.1) are presented in Gueltig et al.'74 
Figure 35 shows the schematic of the device. This study demonstrated that 
MSMA films provide many pathways for energy conversion and intrinsic 
actuation. It is necessary to minimize heat transfer durations while also 
optimizing the effective change in magnetization (AM) in a constrained 
temperature change (AT) in order to make the best use of these multifunc- 
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tional features for thermal energy harvesting. It has been shown that the 
Faraday effect's direct conversion of thermally induced AM to electricity is 
less effective for MSMA film devices than its indirect method, which in- 
volves resonant stimulation of large oscillation amplitudes at high frequency. 
An example MSMA film cantilever device is shown and put to the test for 
various energy conversion configurations. The device produced significant 
power densities when it is self-actuating with mechanical frequency up- 
conversion. The average power density is 0.5 mW/cm? with respect to the 
TM film, which is approaching the range of thermoelectrics, for a AT of 
10 °C. Initially, the MSMA film is heated using a laser to induce AM in the 
MSMA material, this caused an electrical pulse in the pick-up coil that is 
attached to the cantilever tip as shown in Figure 36a. Then in the second 
method, the cantilever is moved using an external actuator to come in contact 
with the heat source to induce AM as depicted in Figure 36b. The magnet is 
attached with the cantilever tip in both cases. 
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Figure 35 Schematic of energy harvester based on MSMA film Nigs-Mn4o-Inis that is presented 
in Gueltig et al.'** Reprinted with permission from John Wiley and Sons. 


Finally, the cantilever is separated from the magnet and the self-actuation 
capability of the MSMA film is used to break contact with the heat source 
and to move towards the magnet. There are two mechanisms contributing to 
the energy conversion. The self-actuating MSMA device resulted in an 
induced change in magnetization (AM) due to thermal energy and also 
additionally due to the changing magnetic field, the overall change in mag- 
netization is higher, which resulted in a higher electrical signal in the pick-up 
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coil. Additionally, the pick-up coil is also experiencing an additional change 
in magnetic field gradient due to the actuation. Thus, this study is the founda- 
tion of energy generation using the self-actuation of magnetic shape memory 
films with a cantilever-based design. MSMA films like Ni4s-Mnao-Inis under- 
go the first-order transformation from a paramagnetic martensite phase to a 
ferromagnetic austenitic phase. This phase transformation is abrupt and large 
but due to the presence of hysteresis, the optimization of the thermal cycle is 
limited. 
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Figure 36 Schematic of the TMG, (a) heating MSMA with laser, (b) external actuation of the 
MSMA film.'™ Reprinted with permission from John Wiley and Sons. 


3.3.3 TMG with N-Mn-Ga 


The second generation of thermomagnetic generators developed at KIT is 
presented in Figure 37.” In this generation, the drawbacks of optimization of 
the thermal cycle are addressed by using a magnetic shape memory alloy 
with second-order transformation (ferromagnetic shape memory alloys). 
Heusler alloy film of composition Nisi.4 -Mnzs.3- Gazo3 and a thickness of 
5 um is used in this device. Optimizing the heat intake and heat dissipation 
from the film by changing the polyamide cantilever to brass to increase the 
thermal dissipation of the device, resonant self-actuation is achieved for the 
first time. This resulted in a very large power output of 2.4 uW at 84 Hz for a 
heat source temperature (Tsource) of 160 °C. 


63 


3 State-of-the-art thermomagnetic devices 


a) b) 


Heatable 


magnet (Tu) Ni-Mn-Ga film 
Temperature sensor 


Ni-Mn-Ga 
film 


N 
a) om Pick-up coil 


Computer/ Laser 
Labview sensor 
Figure 37 Schematic of second-generation TMGs from KIT, using FSMA film Ni-Mn-Ga,” (a) 


measurement setup with the schematic of the device, (b) photograph of the de- 
vice. Reprinted with permission from John Wiley and Sons. 


Cantilever 


Experimental evidence has shown that the TMG is for the first time at eye 
level with the state-of-the-art thermoelectric generators at miniature length 
scales (mm and cm-range). It became possible to gather local thermal energy 
with a small temperature differential thanks to the high power density that is 
observed. A very high power density of 118 mW/cm? is achieved with 
respect to the volume of the active material. Even when the footprint is taken 
into account, the device still has power per footprint of 14.2 uW/cm?. Several 
other iterations are done prior to reaching these milestones, which are pre- 
sented in detail in.”° 


64 


4 Multi-physical analysis of TMG 


The multi-physical LEM simulation model of the TMG is created using the 
MATLAB SIMCAPE toolbox within MATLAB Simulink. Simulink is a 
graphical programming environment with a unidirectional linear flow of data 
through blocks. Although for most simulation models this is acceptable and 
sometimes preferred, in the case of a TMG there are several different parame- 
ters affecting the operation of the TMG and each of them defines the current 
state of the device at every time step. There are several different data loops 
within the model which makes the Simulink toolbox in its native form in- 
compatible with creating the LEM model for the TMG. The SIMSCAPE on 
the other hand has bidirectional data transfer between blocks. SIMSCAPE 
blocks include a variety of partial differential equations (PDEs), linear 
equations, and expressions. The connections and dependencies between these 
equations are defined by the physical ports. Thus, each time step is solved 
numerically to converge. This allows for incorporating multiple data loops 
without following a linear pattern of data flow. Since the blocks are also 
solved numerically additionally to the blocks being connected analytically, it 
allows for modelling complex multiphysics systems in a dynamic environ- 
ment with less computational power. 


4.1 Operation principle 


The TMG device is shown schematically in Figure 38 and is made up of a 
cantilever beam that is affixed to a substrate and has an active TM film and a 
pick-up coil affixed to its freely movable end. The magnetic force that pulls 
the free end of the cantilever towards the magnetic surface is produced by a 
heated magnet, which is subsequently utilized to produce electrical energy. 


‘Section 4.1 has been published in similar form by the author in Joseph et al.” 
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Heat transfer takes place as the cantilever makes contact with the magnet, 
creating a ferromagnetic transition in the TM film. The TM film cools as the 
cantilever oscillates back, which restores the magnetic attraction force. 
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Figure 38 Schematic of the TMG, (a) Side view of the TMG, (b) front view of the TMG. 
Reprinted from Joseph et al. "° 


The result is an uninterrupted oscillation that is maintained by the cantilever 
front's inertia force due to the tip mass of the coil. As a result, the oscillatory 
action aids in cooling the TM film. When cooling and heating result in an 
adequate change in magnetization and associated magnetic actuation force 
within a cycle, resonant self-actuation begins to occur. This is recognized by 
periodic oscillation of the cantilever tip without showing any damped oscilla- 
tions and each cycle is a repetition of the previous one. In this instance, 
thermal energy is most effectively transformed into mechanical and magnetic 
energy. The pick-up coil then transforms this mechanical and magnetic 
energy into electrical energy. 


4.2 Lumped element model 


By combining the mechanical motion of the cantilever front, magnetization 
changes in the MSMA film, and heat flows during film heating and cooling, 
the LEM should represent the energy conversion processes in the TMG. 


Section 4.2 has been published in similar form by the author in Joseph et al." 
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By maintaining the harvester in resonant self-actuation mode, the model 
should enable optimization of the design parameters to achieve maximum 
power and efficiency. For each position and time step during simulation, the 
various physics components exchange data. The LEM parameters used for 
simulations are given in a table for each section of the results. The model is 
divided into 4 parts, mechanical, magnetic, thermal, and electrical sections. 
Each section deals with a part of the TMG cycle. The mechanical section 
deals with oscillations and damping using the magnetic field, magnetic field 
gradient, and magnetic force derived from the magnetic section. The thermal 
part simulates the heat transfer processes and finally, the electrical section 
derives the electrical output. 


4.2.1 Mechanical section 


The center of mass at the movable cantilever beam tip is moved in one 
dimension to approximate the mechanical deflection of the cantilever tip 
while 


Contact Hardstop Magnetic force 


Magnetic section 
Effective (Mag ) 


Electromagnetic 
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(Magnetic section) 
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Figure 39 Mechanical section with inputs of magnetic force and electromagnetic damping 


force and outputs of velocity and position of the beam front. Reprinted from 
Joseph et al.!*8 
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ignoring the lateral displacement brought on by bending. Gravitational force 
affects the movable mass while using the TMG in a vertical direction, how- 
ever, it is insignificant in comparison to magnetic and restoring forces. The 
magnetic force Fmag, inertia force, damping force, and the elastic force of the 
cantilever all contribute to the net force Fret acting on the cantilever front. For 
computing the effective force at the cantilever tip, the section uses the mag- 
netic force and electromagnetic damping force as input parameters. Thus, the 
damping forces brought about by structural damping, viscous air damping 
Fair,!?? and electromagnetic damping Fem!*° are taken into consideration. The 
center of mass’s position and velocity are taken as the output parameters by 
solving the equation of motion considering the influencing forces. 


2 dx 
Fret = M z + C +kx+ F nag = Fam = Fair 
Eq (16) 
The movable mass of the m is taken as:!*! 
33 
m= 140 X Meant + Mcoil + Mfilm 
Eq (17) 


The cantilever tip experiences an additional elastic impact force as it comes 
into contact with the magnet surface, which causes it to abruptly halt move- 
ment at the surface. Figure 39 shows a schematic representation of the 
mechanical section of the LEM. With a damping coefficient of c, the damper 
represents structural damping, and the spring specifies the stiffness of the 
beam cantilever, k. The contact hard stops with contact stiffness and the 
impact damping coefficient provide the impact force at contact. With the use 
of a firm reference, the beam is fastened at one end. The mass serves as a 
point load at the opposing free-moving end. 
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4.2.2 Magnetic section 


The magnetic section is depicted schematically in Figure 40. The tempera- 
ture of the TM film, the location of the cantilever tip, and the current flowing 
through the pick-up coil are input parameters. A lookup table that contains 
experimental data on the relationship between the temperature and magnetic 
field and the magnetization of the TM film enables the simulation to choose 
the magnetization value that is most suited given the inputs for temperature 
and magnetic field as there is no hysteresis in transformation. The magnetic 
force Fmag on the TM film is calculated with the inclusion of magnetization. 
The magnetic field is calculated along the symmetry axis z direction to 
determine Fmag, using an adapted version of Biot Savart law.!?*133 


Z (magnetic easy axis) 


Br ab 
B(z) = — jarctan (——==) 
T 2zV4z2 + a2 + b2 


e] 
— arctan | ————— 
2(c + z)y4(c +z)? + a? + b? 
Eq (18) 
Eq ( 18 ) lets us calculate the magnetic field at a distance z from the magnet 


surface. This can also be extended to three dimensions if required.'” Taking 
the derivative of the equation Eq ( 18 ) with respect to z gives us Eq ( 19 ), 
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which is the magnetic field gradient at any distance z. Here the assumption is 
made that the cantilever deflection is a 1- dimensional motion. 


Current (Electrical section) 


Magnetic field 


2-D Lookup table 


| Position 


(Mechanical section) 


E Temperature 
Thermal section) 


Magnetic force 
Magnetic field gradient 
Electromagnetic damping force 
Figure 40 Magnetic Section with the temperature of the TM film and position of the beam front 


as inputs. The magnetic force, magnetic field gradient, and electromagnetic 
damping force are the outputs of the section. Reprinted from Joseph et al." 


OB(z) 
0z 
2 ab z 2ab 
_ Br 222422 + a? + b? (422 +22 + b2)2 
= 2h2 
i P ee +b?) +1 
= ab _ 2ab 
2(c + z)VJ4(c +z)? + a2 + b? (4(c+z)2 +a2 + b2)2 


u a?b? 
4(c +z)? (4(c + 2z)? +a? + b?) v4 


Eq (19) 
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oB 
=V, 


F magMr 2 


mag 
Eq (20) 


Equations Eq ( 18 ) and Eq ( 19 ) are used to predict analytically the magnetic 
field Bz and field gradient (OB,/dz) as a function of the position of the 
MSMA film (distance from the magnet surface along the z-direction). The 
TM film volume is taken as Vmag. Using the generated current in the pick-up 
coil and magnetic field as inputs, the electromagnetic damping is calculated. 
Magnetic force Fmag, magnetic field gradient, and electromagnetic damping 
force are the output parameters. These forces depend on the temperature and 
position of the TM film. Eq ( 20 ) gives the final magnetic force acting on 
the magnetizable material. 


4.2.3 Thermal section 


The numerous heat transfer processes that take place in the device are sche- 
matically shown in Figure 41. Heat transfer between a heat source and TM 
film (process 1) depends on the contact heat transfer coefficient (hr), and the 
temperature difference between the Tsource and the temperature of TM film 
(Trim). A large contact area is an advantage of thin film geometry due to its 
large surface-area-to-volume ratio. The permanent magnet's polished surface 
results in a smooth surface finish, which is crucial for optimum thermal 
contact. Heat convection from the TM film to the surrounding air via the top 
and side surface is denoted by processes 2 and 3, respectively. Along with the 
temperature difference between the TM film and the surrounding air, the 
convective heat transfer coefficient, which is heavily influenced by the speed 
of the cantilever tip determines these processes. When the Tsource rises, the 
oscillation frequency increases, and the heating time required to achieve a 
significant AT decreases.”® Because the cantilever front moves more quickly 
with increasing frequency, there is a rise in convective heat transfer, which 
enables the device to adapt to changes in the temperature of the heat source. 
The thermal section of the LEM simulation models the heat flow and thermal 
interaction between different components of the TMG. The common ap- 
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proach for modeling heat flow is by utilizing either a Foster or Cauer model 
in which an electrical analogy is used with resistance and capacitors repre- 
senting thermal resistance and thermal capacitances of different compo- 
nents.”? 


(®©) Heat transfer from heat source to TM film 

© Convection from film to ambient air 
(enhanced by oscillation) 

© Lateral convection to ambient air (TM film and bonding layer) © © TM film 

© Heat transfer from film to bonding layer 

©) Heat transfer from bonding layer to cantilever = y Bonding layer 

© Heat conduction through cantilever 


© Convection from cantilever to ambient air mO! @= Cantilever 
(enhanced by oscillation) SULIT 
Heat transfer to ceramic substrate oO Pick-up coil 


Figure 41 Heat transfer processes in the TMG assembly. Reprinted from Joseph et al.?' 


Although both models can give an accurate heat flow, the Foster model does 
not represent the system on the component level, thus it is hard to model the 
heat flow path within the components and also to know the temperature of 
each component. Hence in the presented LEM model, a Cauer thermal model, 
which is a continuous fraction model, is used. This allows modeling of 
thermal behavior and performance on a component level within the TMG. 


The TMG’s thermal section is shown schematically in Figure 42 with differ- 
ent heat transfer processes marked with a box that corresponds to Figure 41. 
Eq ( 21 ), which gives the overall thermal impedance Zin of the device, 
represents the continuous fraction equation of the thermal network used in 
the work: Thereby ‘n’ is the total number of nodes, and Cmo) and Rund) are the 
thermal capacitance and resistance of a given node, respectively. A section of 
the device assembly is represented by each node in the model; for example, 
the first section would be the TM film, the second would be the non- 
conductive bonding layer, etc. 


72 


4.2 Lumped element model 


TM =TM film 
BL _=Bonding layer Tamb 


C_f =Cantilever front 


C_b =Cantilever body 
conv =Convective thermal resistance C_Dyon 
res =Conductive thermal resistance t 


cont =Contact thermal resistance 
Tambo =Ambient temperature of air 
Tsource =Temperature of heat source 
Tsu  =Temperature of substrate T. 

Sw _=Switch heat source & ambient sub 


Temperature 
Position 
(Mechanical section) 


Figure 42 LEM of thermal section takes in position as input and generates the temperature of 
different components (like film, cantilever, etc.) as the outputs. Reprinted from 
Joseph et al.”! 


Zin = 
Cna) + 
Rat, 


1 


th(2 
@)+R 7 


th(2)+-—— 
ChM tR 


Eq (21) 


Using this method, it is also possible to think of individual parts as a combi- 
nation of several thermal resistors and capacitors, allowing us to calculate the 
temperature gradient across each component of the device with a resolution 
equal to the number of thermal resistors and capacitors. However, only one 
thermal resistor and thermal capacitor are utilized for each part generally 
because the heat flow and average layer temperature are of greater im- 
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portance. Depending on its dimensions and thermal characteristics, each 
section's AT can be calculated and the heat flow direction can be identified.'** 


4.2.4 Electrical section 


The schematic of the electrical section is shown in Figure 43. It computes the 
induced voltage Vina using Faraday's law based on the inputs of the magnetic 
field gradient and pick-up coil velocity. 


Eq (22) 
Nooit and Acoil are used to denote the pick-up coil's area and the number of 


f a i ; . 
turns, respectively. x represents the film/pick-up coil velocity. 


Current 


£ Velocity 
(Mechanical section) 

Magnetic field gradient 

(Magnetic section) 


Load Coil 
resistance parameters 


Figure 43 Electrical section takes in the velocity of the beam front and magnetic field gradient 
as the input to the section and generated current is given out as output. Pick-up 
coil parameters and Rı do not change during the simulation of one configura- 
tion and hence are given as constants. Reprinted from Joseph et al." 
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The electrical section's output, the induced current is calculated from Vina, 
and the internal pick-up coil resistance as a function of load resistance (Rt). 
As a minor effect in the case of TM film, the change in flux density brought 
on by the change in magnetization of the TM film is disregarded. This is due 
to the low volume of the active material, which causes a change in magnetic 
flux density which is several orders less than the change due to the stroke of 
the beam. 


4.3 Experimental setup 


A modular measurement setup is assembled using position stages and adapt- 
ers. The magnet holder is designed to achieve 360° rotation while keeping an 
edge of the magnet at the same position. The position of the magnet in 
relation to the sample should only change in one direction when adjusting to 
achieve proper contact with the TM film. Achieving proper physical contact 
between the TM film and magnet surface is significant to having a large 
value of contact heat transfer coefficient (hr). The thin films (Ni-Mn-Ga and 
Gadolinium) are usually smooth as they are fabricated by magnetron sputter- 
ing, the magnet surface in contact is wet polished to a mirror finish to achieve 
maximum contact. In some fabrication processes, like the melt spinning used 
for La-Fe-Si-H samples, the surface of the TM film prepared is rough and 
hence the value of hr at contact is lower compared to the sputter deposited 
films. 


Figure 44 illustrates the electromechanical circuit diagram of the TMG. A 
resistive heater built into the magnet holder is used to heat the permanent 
magnet, and an RTD sensor (PT100) affixed close to the magnet's tip is used 
to monitor its temperature. The RTD sensor is connected to a temperature 
sensor module USB-TEMP from MCC, and temperature data retrieved using 
InstaCal software on the PC. For the purpose of measuring deflection, a 
Panasonic HL-G103SJ laser triangulation sensor is employed. The sensor is 
connected to a USB 6211 A/D card from National Instruments and is used in 
analog mode. The cantilever moves in two dimensions because it is bending 
upward, which causes a measuring error that is equivalent to a 5 percent error 
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in terms of the entire stroke. To measure the generator's output, a load resis- 
tor (RL) is attached to the TMG. To measure the current flowing through the 
Rx, a current pre-amplifier with the model number SR570 from Stanford 
Research Systems is connected in series. Under low noise conditions, a 
sensitivity of 100 uA/V is used, with a resolution of 60 pA. Additionally, an 
A/D card is attached to the output of the current pre-amplifier. The USB 6211 
A/D card is connected to a PC running LabVIEW, which is used for data 
acquisition at 5000 Hz and performs basic data processing, including the 
computation of instantaneous power from the measured current, which is 
exported as an excel file (xlsx format). Data is recorded for a duration of 
either 3s or 10s, and 10 sets of such measurements are taken for every point 
of operation of the TMG (for averaging). Different loads are connected to the 
TMG to investigate the optimum Rz for the maximum power output of the 


TMG. A python script is used to analyze the data. Average power is calculat- 
ed by integrating instantaneous power for the entire time and dividing it by 
the duration of data acquisition. By doing peak to peak analysis of deflection 
measurement data, the stroke is computed and averaged over the entire 
measurement. 


heated permanent 
-magnet 


„MSMA film stack 
“cantilever 
Im 


Y “~~ pick-up coil 


Figure 44 Electromechanical circuit diagram of a TMG. A spring-mass-damper system with 
stiffness K and damping constant D make up the mechanical component. The 
pick-up coil's internal resistance (Ri) and the Rı are both included in the elec- 
trical component. Im stands for the induced electrical current. Adapted from Jo- 
seph et al.” 
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The frequency of the device operation is determined by converting deflection 
measurements from the time domain into the frequency domain by doing an 
FFT analysis. Power, stroke, and frequency are calculated for ten measure- 
ments taken consecutively and then the mean is calculated to get the average 
power, average stroke, and average frequency for that measurement point. 
Figure 45 shows a schematic model of the experimental setup used to char- 
acterize samples. 
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Figure 45 3D model of the experimental setup used to measure the output of the thermomag- 
netic generators, allowing for several degrees of freedom. 


4.4 TMG with Ni-Mn-Ga 


4.4.1 Material properties 


Heusler alloys are ordered intermetallics that display multiferroic phase 


transformations because electronic spin interactions are highly influenced by 


variations in interatomic distances and local symmetry.*8°!,6- 135.136 


“Section 4.4 has been published in similar form by the author in Joseph et al.”*'7 
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The alloys Ni-Mn-Z and Ni-Co-Mn-Z (Z = Ga, In, Sn...) are notable exam- 
ples. These materials display strong effects such as magnetic field-induced 
strains,*%13” magnetoresistances,*” magnetocaloric effects,’ and magnetic- 
field-induced shape memory effects.*! In this study, the focus is on 
Nis3.sMn23.sGa22.7 films, 
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Figure 46 (a) DSC measurement of Ni-Mn-Ga alloy (b) Magnetization measurement for 
varying temperatures for Ni-Mn-Ga film. Reprinted from Joseph et al.?® 


which exhibit both a first-order martensitic transformation and a second-order 
ferromagn-etic transition. By using RF magnetron sputtering to create the Ni- 
Mn-Ga films, the Ni concentration and, consequently, the phase transition 
temperatures, are controlled.'**"° Due to limitations in sputtering duration 
and associated film homogeneity, the film thickness has usually a thickness 
below 10 um. The material properties are finalized by heat treatment of 10 
hours at 800 °C. Differential scanning calorimetry is used to determine the 
film’s phase change over temperatures (DSC). Figure 46a depicts the DSC 
measurement for the Ni-Mn-Ga film. With transformation temperatures 
between 57°C and 83 °C, the peaks represent the first-order martensitic 
transition. 


The fact that the transition temperatures are heavily influenced by chemical 
composition and the associated electron concentration per atomic unit (e/a- 
ratio) provides for a wide range of tailoring.'4! A temperature-dependent 
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magnetization measurement of a Ni-Mn-Ga film at 0.05 T ofthe low magnet- 
ic field is shown in Figure 46b. The material is in the martensitic and ferro- 
magnetic state at room temperature. The martensitic phase change is respon- 
sible for a small step-like characteristic that appears at around 70 °C. The Tc 
occurs at about 102 °C, the magnetization undergoes a rapid ferromagnetic 
transition with a substantial change in magnetization (AM) occurring within a 
small AT which is much less influenced by chemical composition.'”° Due to 
the second-order nature of this transition, it does not show hysteresis. To 
drive the TMG’s resonant self-actuation in a magnetic field gradient at a low 
magnetic field (<0.5 T), we use the abrupt change of AM/AT at Tc. The 
magneto-caloric effect is thought to be minimal in situations with a low 
magnetic field and separated first-order and second-order transitions. 


4.4.2 Fabrication 


In a high-purity argon environment, RF magnetron sputtering is used to 
create the Ni-Mn-Ga films. The Ni concentration and subsequently the phase 
transformation temperatures are controlled by adjusting the sputtering 
power.!?® To customize the temperature-dependent change in magnetization 
at the Tc, a heat treatment lasting 10 hours is conducted at 800 °C. By using 
an inductively coupled plasma technique, it is discovered that the crystallized 
freestanding films are composed of Nis3.sMn23.sGaz22.7. The sputtering time is 
regulated to produce a layer with a thickness between 5 and 10 um. Ohtsuka 
et al.'* provide information about the films within the experiment's structure. 
The pick-up coils are created by sandwiching a 250 um thick PMMA core 
with a 1 x 1 mm? surface area between two 25 um thick Kapton foils with a 2 
x 2mm?’ surface area. Using a semi-automatic setup (Semi-automatic coil 
winder) and insulated copper wire with a diameter of 15 um, the pick-up 
coils are fabricated. 


The device is fabricated in five steps as shown in Figure 47. Step I show the 
ceramic substrate which is prepared using laser cutting from an alumina sheet 
using a pulsed laser. The brass cantilever is fabricated using laser processing 
by an ultra-short pulse femtosecond laser from a brass sheet. In step II this 
cantilever is attached to the ceramic substrate. The pick-up coil fabricated by 
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the semi-automatic coil winder is attached to the underside of the cantilever 
tip in step III. Using a non-conductive bonding layer, the TM film is attached 
to the cantilever tip in step IV. 


Ill Iv A 
N Ry 


Figure 47 Fabrication steps for preparing a TMG prototype. 


Finally, in step V the connector bush pins are added and attached to the pick- 
up coil. The connector bush pin allows external connections to the pick-up 
coil for electrical measurements. The completed device assembly is placed 
under a permanent magnet for providing the magnetic field to magnetize the 
TM film and the magnetic field gradient for the magnetic force. Using a 
digital microscope, the magnet is placed in relation to the Ni-Mn-Ga film 
surface at the ideal distance and angle to guarantee maximum contact while 
minimizing the impact on the magnet surface during operation. The wires are 
only 15 um, which is thin enough to avoid any mechanical effects on the 
cantilever motion. Before conducting experiments, the magnet is heated until 
a stationary temperature is attained. Data are gathered by varying Ru of each 
reading while maintaining a constant temperature at source. For every set of 
Tsource, these procedures are repeated. 
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4.4.3 Device characterization 


Since several different TMGs are fabricated and performance is analyzed for 
all of them, it is necessary to have a standard performance characterization 
for the device so that different materials could be compared with each other. 
The size of the magnet is kept constant (3 mm x 3 mm x 8 mm) to avoid the 
effect of the magnet on the performance comparison between the devices. To 
keep consistency, the lateral dimensions of the TM film are also kept con- 
stant, the magnetic force is only tuned by changing the thickness of the film. 
The width of the cantilever is kept at 2 mm to avoid any overhanging which 
increases air damping. The stiffness of the cantilever is adjusted by changing 
the length of the beam, if the length of the beam exceeds 6 mm, the thickness 
of the cantilever is reduced. 
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Figure 48 Effect of Ri on the performance of the TMG (a) stroke of deflection, (b) frequency 
of oscillation, (c) power output from the pick-up coil. Part of data pre- 
published in Joseph et al.78°!!°8 
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If the length of the cantilever required for appropriate stiffness to balance 
mechanical and thermal cycle is less than 4.3, the thickness of the beam is 
increased, and the length is tuned accordingly. This is to avoid increasing the 
stroke more than 2.5 mm (in the case of long cantilevers) and to avoid a very 
large bending angle within the cantilever (for short cantilevers). A stroke of 
more than 2.5 mm would result in a very low magnetic field gradient. This is 
explained in detail in section 5.1. A large bending angle would result in 
higher stress within the cantilever beam causing fatigue and other undesirable 
effects. The dependence of performance of the TMG on Ru and Tsource is 
shown in Figure 48. Figure 48 shows the performance of a TMG with a 
magnetron sputtered Ni-Mn-Ga film as the TM film. The composition of the 
film is Nis3.:-Mn23.s-Ga22.7. 


The thickness of the film is 10 um with a lateral dimension of 2 mm x 2 mm. 
The cantilever is a single beam with a lateral dimension of 2 mm x 5.7 mm 
and a thickness of 20 um. The cantilever is attached to the substrate using a 
non-conductive bonding layer. The TMG is operated at a Tamb of 23 °C. The 
maximum performance is reached at a Tsource of 170 °C and an Rı value of 
400 Ohm. The power output increase as the Rı is increased from 200 Ohm 
and reached a peak value of 3.1 uW. The power output shows a reduction 
when the Rı is increased further. The stroke increases as the Tsource is in- 
creased. The stoke continuously shows an increase as the Rı is increased 
from 200 to 600 Ohm without presenting an optimum value of the Rt. Fre- 
quency is also affected by the change in Rı and shows an increase as the 
value of Rı is increased. Although both the stroke and frequency increase 
with an increase in Tsource, the transfer efficiency decreases as the Rı is 
increasing, as the value is moving further away from the internal resistance of 
the device. The performance of the device increases as the current through 
the pick-up coil is reduced when Rx increases, causing a reduction of elec- 
tromagnetic damping. Thus, these contradicting features cause the device to 
have an optimum Rı of 400 Ohm at which the maximum power output is 
observed. 
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4.4.4 Model validation 


For a Tsource of 170 °C and an Rı of 400 Ohm, time-resolved experimental and 
simulated properties of center-of-mass displacement and corresponding 
electrical current induced in the pick-up coil are examined for a TMG with 
Ni-Mn-Ga as TM film are shown in Figure 49. The simulation model is 
validated by the good agreement between simulation and experiment results. 
The results shown are obtained under stationary operating conditions, which 
are attained after a short period of time in continuous resonant oscillation and 
then last for hours with no discernible performance change. 
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Figure 49 Model validation for the case of the TMG with NiMnGa at a Tsource of 170 °C and an 
Rx of 400 Ohm (a) displacement measurement, (b) LEM simulation of dis- 
placement, (c) current measurement, (d) LEM simulation of current. Reprinted 


from Joseph et al.’ 
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Once the time resolved performance of the simulated the TMG is matched 
with experiments, it allows for the estimation of the contact heat transfer 
coefficient (hr) and the heat transfer coefficient of convection. These values 
are dependent on the surface conditions, impact force, velocity of the cantile- 
ver tip and the ambient conditions. These values could not be measured 
experimentally without affecting the performance ofthe TMG. 


4.5 Thermal cycle in TMG under resonant self- 
actuation 


The LEM is used to investigate the time-dependent variations in temperature, 
magnetic field, magnetic field gradient, and associated magnetization of TM 
film in a TMG with Ni-Mn-Ga as the active TM film during TM cycling. 


Cantilever motion can be divided into four phases, as shown in Figure 50a. 
Stage (1) corresponds to the instance when the Ni-Mn-Ga film comes into 
mechanical contact with the heat source, resulting in the transfer of heat into 
the film. The magnetization of the film is reduced due to the increase in 
temperature of the film, even though, the magnetic field experienced by the 
film is constant as can be inferred from Figure 50d and Figure 50c respec- 
tively. Stage (2) represents the instance when the Ni-Mn-Ga film and pick-up 
coil are accelerated to maximum speed upon deflection from the magnet, 
which corresponds to a maximum induced voltage. At the greatest distance 
from the magnet, stage (3) is attained which also corresponds to the lowest 
field and magnetization within the TM film. Stage (4) denotes a second 
maximum induced voltage that occurs in the opposite direction when the 
pick-up coil is moving at its fastest speed towards the magnet as it attracts the 
Ni-Mn-Ga film. The actual cycle is more gradual with increasing magnetiza- 
tion over the entire oscillation, unlike the case of an ideal cycle with a heat 
sink causing the magnetization to increase suddenly at stage (3). Figure 50b 
shows the time-resolved AT that occurs during resonant self-actuation in the 
Ni-Mn-Ga film. 


‘Section 4.5 has been published in similar form by the author in Joseph et al.'* 
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The film depicts a AT of around 10 °C in a period of time of about 8.7 ms 
during one oscillation cycle. The temperature of the film exhibits a sharp rise 
within approximately I ms while in mechanical contact (1). The subsequent 
temperature drop happens more gradually as it moves through stage (2) - (4). 
Rapid heat transfer during mechanical contact with the heat source makes the 
relatively substantial AT possible, whereas heat loss by heat conduction and 
convection is constrained. The temperature difference between the Ni-Mn-Ga 
film's lowest temperature and the cantilever below it is less than 0.5 °C. The 
cantilever temperature rises over time and saturates with very small tempera- 
ture cycling. The film cools by both conduction and convection, as the film 
approaches the temperature of the cantilever, the conduction heat transfer is 
reduced, but if the oscillation is still not complete the convection continues as 
the ambient temperature is at a much lower temperature. 
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Figure 50 (a) Device operating stages (b) temperature change within Ni-Mn-Ga film, (c) 
magnetic field change, (d) magnetization change in the Ni-Mn-Ga film. Re- 
printed from Joseph et al. "8 
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This causes the Ni-Mn-Ga in some cases to cool down to a temperature 
slightly below cantilever temperature. This is aided by the thermal resistance 
of the bonding layer that prevents a sudden heat flow from cantilever back 
into the film. 


4.5.1 Achieving resonant self-actuation by design 


To achieve resonant self-actuation, it is crucial to tune the TMG parameters 
like the volume of the active material, cantilever dimensions, bonding layer 
thermal resistance and coil parameters (size and number of turns) according 
to the TM film properties like the magnetization under the applied magnetic 
fields and the change in magnetization per change in temperature (AM/AT). 
From the magnetization measurements, the AM/AT values can be obtained by 
taking derivative of the magnetization based on temperature. Figure 5la 
shows the typical magnetization curve for Ni-Mn-Ga and Figure 51b shows 
the derived AM/AT. The peaks in the AM/AT values show the point of opera- 
tion for the device where the AT in the material should be concentrated to 
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Figure 51 Typical properties of Ni-Mn-Ga. (a) Magnetization of the Ni-Mn-Ga film, (b) AM/AT 
derived from magnetization to observe the change in magnetization with re- 
spect to temperature change. The data for this graph was extracted from Guel- 
tig et al.” 
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achieve the largest change in magnetization for a small AT. The key is to 
adjust the stiffness, bonding layer thermal resistance (Ro), and Tsource range 
accordingly to achieve the AT in the TM film around the peaks in AM/AT 
values. Achieving this ensures that the TMG is operating with maximum 
change in magnetization allowing the device to operate at high frequencies. 


The LEM model is supplied with the measured magnetization values as a 
lookup table. Then the parameter ranges are supplied to an optimization 
program that finds the operating range of temperatures under which the 
device is in resonant self-actuation and then further optimizes the perfor- 
mance by adjusting the parameters. The procedure for optimization is illus- 
trated using a schematic in Figure 52.Once the optimal parameters are found, 
like the R» and the stiffness of the cantilever beam, a device prototype is 
fabricated with the optimized parameters by keeping tolerance to the mini- 
mum. Optimal parameters are determined by the maximum power output 
from the device under the constraints provided. 


After conducting experiments, the simulation is readjusted to match the heat 
transfer coefficient of convection (hr) and convective heat transfer coefficient 
in experiments as these values are assumed in the beginning. The TM film 
properties determine the optimum frequency of operation, stroke of the 
device, and operating temperature range of the device. Although frequency 
can be modified by changing the stiffness, the optimum frequency would still 
be determined by the TM film. The key challenge is to design TMGs with 
optimum heat transfer such that heat intake occurs in a limited contact time 
with the magnet and heat dissipation occurs across the balance of the cycle to 
provide regular and periodic AT. If AT in the TM film is not regular, the 
device will overheat, resulting in an extra oscillation or a full halt. An exam- 
ple of additional oscillation and operation halt due to heat buildup is given in 
Figure 53 a and b respectively. It can be observed that the additional oscilla- 
tion is caused by the mismatch between the thermal cycle and mechanical 
oscillation. Ideally, the device should perform with maximum oscillation and 
velocity. 
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Figure 52 Flowchart showing the optimization procedure to obtain the parameters like the 
operating temperature range, dimensions of cantilever and bonding layer ther- 
mal resistance. 
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4.5 Thermal cycle in TMG under resonant self-actuation 


Since the velocity depends on the AT of the film which affects the magnetic 
force (stiffness is constant and hence will not contribute to a change in 
velocity) higher change in magnetization will result in higher velocity 
(matching the tip velocity at the natural frequency of the beam). When the 
ideal conditions change, the device still exhibits resonant self-actuation by 
adjusting its frequency to match the thermal cycle. The thermal cycle is 
affected by the AT in the material, which also determines the magnetic force, 
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Figure 53 Resonant self-actuation (a) with additional oscillation due to uneven heating and 
cooling cycles (b) heat buildup causing the device to stop operating. 


the velocity of the tip and the stroke. However, this window of operation 
depends a lot on the properties of TM film, Tamb, Tsource, and device design. If 
the TM film exhibits a gradual and continuous drop in magnetization, then 
the device shows much more tolerance to change in ideal conditions, but if 
the material is exhibiting an abrupt drop in magnetization over a few Kelvin, 
then the TMG performance would be heavily affected by the change in ideal 
conditions. The frequency of operation also changes during operation in the 
case of non-resonant self-actuation as the thermal cycle is not in balance with 
the mechanical cycle. Figure 54a shows 
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Figure 54 FFT of time domain current measurement, (a) non-resonant self-actuation with 
multiple peaks, (b) resonant self-actuation with a single large peak. 


the Fast Fourier Transform (FFT) for non-resonant self-actuation and it is 
clear that the frequency within oscillations is spread over a span. In the case 
of resonant self-actuation, the frequency of operation is much more stable 
and does not change during the operation provided the source and ambient 
conditions remain the same. Figure 54b shows the FFT for resonant self- 
actuation with a large single peak. 


4.6 Discussion 


The coupled thermo-magneto-mechanical performance of the TMG is de- 
scribed and each physical component of the device is examined using the 
LEM approach. Here, we concentrate on the TM film's AT, the magnetic 
field's dependency on the film's position, and the associated film magnetiza- 
tion which results in changes to the magnetic actuation force needed to 
maintain resonant self-actuation. At contact, the time-dependent temperature 
profile shows a quick temperature rise of 10 °C in just I ms, followed by a 
temperature fall in just 8 ms. This performance is caused by the rapid heat 
absorption during mechanical contact between the magnet and the Ni-Mn-Ga 
film and the slower heat dissipation through convection and conduction via 
the bonding layer and cantilever. The non-conductive bonding layer between 
the cantilever and the Ni-Mn-Ga film serves as a thermal boundary layer and 
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helps to retain the majority of the input heat, which is why the temperature 
fall occurs considerably more slowly. It is required to assume the value of hr 
between the magnet and the film as well as the convective heat transfer 
coefficient because these quantities are challenging to measure and challeng- 
ing to theoretically deduce due to a large number of unknown variables. 


Table 1 Modeling parameters for LEM simulation of the TMG. Part of data pre-published in 
Joseph et al.7°118 


Mechanical parameter Unit Value References 
Pick-up coil mass Kg 6 x 107° This work 
Cantilever density Kg/m? 8500 142 
TM film density Kg/m? 8020 u 
Bonding layer density Kg/m? 1250 129 
Structural and air damping Ns/m 2.8 x 1075 This work 
Impact damping Ns/m 0.1 This work 
Contact stiffness N/m 10000 This work 
Young’s modulus Pa 1 x 1011 me 
Length of the Ni-Mn-Ga film mm 2 This work 
Width of the Ni-Mn-Ga film mm 2 This work 
Thickness of the Ni-Mn-Ga um 10 This work 
film 

Length of the Cantilever mm 5.7 This work 
beam 
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Width of the Cantilever mm 2 This work 

beam 

Thickness of the Cantilever um 20 This work 

beam 

Electrical & Magnetic Unit Value References 

Parameters 

Remanent magnetic field T 1.07 This work 

Number of turns - 400 This work 

Area of pick-up coil m? 1.96 This work 
x 1076 

Pick-up coil internal re- Q 250 This work 

sistance 

Load resistance (RL) Q 400 This work 

Thermal Parameters Unit Value References 

Cantilever thermal conduc- W/mK 109 a 

tivity 

TM film thermal conductivi- W/mK 23.2 a 

ty 

Bonding layer thermal W/mK 0.33 = 

conductivity 

Cantilever heat capacity J/kgK 400 A3 
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TM film heat capacity J/kgK 490 us 
Bonding layer heat capacity J/kgK 2100 12 


Contact heat transfer coeffi- W/m’K 7700-8400 This work 
cient 


Convective heat transfer W/m?K 40-140 This work 
coefficient 
Operation temperature °C 130-170 This work 


Table 1 shows the parameters that have been utilized in the LEM to repro- 
duce the results of the experiment. By matching the simulation with the 
experiments, values of contact heat transfer coefficient and convective heat 
transfer coefficient have been obtained which shows good agreement with 
values suggested from the literature." The actuation force which is deter- 
mined by the magnetic field gradient, magnetization, and the impact force at 
the point of contact with the magnet are both factors that affect the heat 
intake. The maximum attraction force is 5.48mN due to the maximum 
magnetization of 4.7x10° A/m at the point of contact with the magnet. While 
the temperature rises adiabatically from about 91 °C by about 10 °C during 
contact, the magnetization rapidly decreases by about 10° A/m. The magneti- 
zation value depicted in Figure 7 for Ni-Mn-Ga film has a high value of 
AM/AT for unit volume, which is approximately 10* A/(mK) for the volume 
of the TM film. It turns out that the best performance takes place when the 
ferromagnetic transition first starts, while the magnetization is still high. In 
order to complete the TM cycle, sufficiently strong magnetic forces are 
therefore preserved. The cantilever's stiffness is adjusted to balance out the 
forces of inertia and actuation, resulting in a relatively high operation fre- 
quency of 114 Hz. 


The Rı has severe effects on the performance as it causes a variation in the 
current generated from the pick-up coil resulting in varied electromagnetic 
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damping. The transfer efficiency is maximum when the Rı is matched to the 
internal resistance of the pick-up coil, which would be 200 Ohm in the 
presented case. But it can be seen that the power output increases as a result 
of the increase in both frequency and stroke as the Rı increased from 200 to 
300 Ohm. This is due to the decreased current in the circuit resulting in 
reduced electromagnetic damping. Thus, the power increases with an increase 
in Rx even though the transfer efficiency is still reduced when the difference 
between Rı and internal resistance increases. Thus, after increasing the Rı 
beyond 400 Ohm the frequency and stroke still increase but the power output 
is reduced as the increase in frequency and stroke is not enough to compen- 
sate for the subsequent loss in transfer efficiency. This causes the TMG to 
have a peak performance at 400 Ohm in the presented configuration. 


Table2 Presents the experimental and simulation findings for the TMG under investigation 
utilizing a Ni-Mn-Ga film with a 10 um bonding layer. According to the impact 
force at the magnet's surface, as previously demonstrated, the contact hr values are 
shifting somewhat.” Data pre-published in Joseph et al.?® 


Tsource Power Frequency Stroke Heat | Heat dissi- | AT hf 
intake pation 
(°C) [W/ 
CC) | (uW) (Hz) (mm) (mW) (mW) (m’K)] 
Exp Sim Exp Sim | Exp Sim Cond Conv 
170 | 3.08 3.15) 114.1 114.3 1.77 1.78 228 192.3 35.7 | 10.2 8400 
160 | 2.71, 2.72 110.6 110.9 1.67 1.68 196.8 1674 29.4 | 8.79 8350 
150 | 182, 182 991 99.5 | 1.62 1.59 145.8 1277 181 7.56 8000 


140 | 1.32) 137 91 91.2 1.56 1.55 120.5 107.5 13 6.93 7850 


130 | 1.05 1.07 80 82.7 1.52 149 1061| 953 | 10.8 6.79 7700 
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Parts of this section has been published in Joseph et al.’! There are several 
parameters like magnetic field, field gradient, stiffness, pick-up coil turns, 
contact heat transfer coefficient (hr), and bonding layer thermal resistance 
(Ro) that affect the performance of the device. These parameters also deter- 
mine if the device is operating in resonant self-actuation mode, as the balance 


between all of these parameters is required for achieving resonant self- 
actuation. 


5.1 Magnet 
5.1.1 Effect of distance from the magnet surface 
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Figure 55 Magnetic properties of Samarium Cobalt (Sm2Co17) magnet with cross-sectional 
area 3 mm x 3 mm and length of 8 mm at a distance of 1.94 mm from the mag- 
net, (a) magnetic field measured from the center (0 mm) and moving towards 


the edge (1.5 mm), (b) magnetic field gradient derived from the magnetic field 
measurement. 
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Both the magnetic field and magnetic field gradient can be adjusted by 
changing the dimensions of the magnet and by changing the composition of 
the magnetic material which determines the remanent magnetic field. The 
higher the remanent magnetic field, the higher will be the magnetic field and 
field gradient around the magnet. This is clear from the equations Eq ( 18 ) 
and Eq ( 19 ). The material for the magnet is chosen as SmCo with a maxi- 
mum operating temperature of 350 °C and a remanent magnetic field of 1.07 
T. The drop in magnetic flux density is only 0.05 % per 1 °C, showing high 
stability in high-temperature operation. Although it is possible to use Neo- 
dymium magnets (NdFeB) for low-temperature applications below 80 °C, the 
same SmCo magnets are used to keep measurement consistency. Also, the 
magnetic flux density for a NdFeB magnet changes 0.11 % for every 1°C 
change in its temperature. The lateral dimensions of the magnet used in the 
presented work are 3 mm x 3 mm, and a length of 8 mm. Figure 55 shows 
the magnetic field measurements using a precision Gauss meter and resulting 
magnetic field gradient values. 
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Figure 56 Direct comparison of magnetic field measurements and magnetic field simulation. (a) 
Magnetic field along the z direction at the center determined by measurement 
and simulation, (b) magnetic field gradient derived from magnetic field meas- 
urement and simulation along the z direction at the center. 
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The simulated magnetic field and magnetic field gradients are shown in 
Figure 56, and are compared to the measured values for the distance between 
1.9mm and 5.7 mm for which magnetic field measurements are available, 
The measurements are taken from 1.94 mm away from the magnet, this is 
due to the limitation to get the probe closer in the precision Gauss meter. 
Magnetic field simulation from the LEM model could be used to analyze the 
magnetic field along the z direction by considering it as 1-dimensional. 
Figure 56 shows the computed magnetic field and magnetic field gradient 
using LEM simulation. The magnetic field measurements and the simulation 
values are almost identical, thus, validating the magnetic section of the LEM 
simulation. 


5.1.2 Effect of dimensions ofthe magnet 


Using the LEM simulation, it is possible to study how the magnetic field and 
magnetic field gradient are affected by the change in the cross-sectional area 
and the length of the magnet along easy axis. A permanent magnet's magnetic 
easy axis is the direction in which magnetization requires the least energy. 
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Figure 57 LEM simulation of the effect of dimensions of magnet on (a) magnetic field, and (b) 
magnetic field gradient, at the surface of the magnet. 
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In a ferromagnetic material, such as a magnet, with a preferred orientation for 
the magnetic domains within it, it refers to the preferred direction of magneti- 
zation. The material's asymmetry and crystal structure both affect the easy 
axis. As it affects the direction of the magnetic field the magnet produces, 
this easy axis is crucial for understanding how permanent magnets behave in 
applications. Figure 57 shows the effect of dimensions on the magnetic field 
and magnetic field gradient at the surface of the magnet. The magnetic field 
increases when the length of the magnet along the easy axis is increased. The 
number of magnetic domains aligned along the easy axis rises as the length 
along the easy axis' increases, resulting in a higher degree of magnetization. 
A larger magnetic field results from this in turn. As a result, if the magnetiza- 
tion stays constant, the magnetic field of a permanent magnet will grow as 
the length along easy axis does as well. But it decreases as the area is in- 
creased. The same effect is seen on the magnetic field gradient. The depend- 
ence of the magnetic field and the magnetic field gradient on the dimensions 
is crucial in designing a TMG for an application. The increase in the magnet- 
ic field causes an increase in magnetization of the TM film and the increase 
in the magnetic field gradient results in an increase in the magnetic force. 
This shows that having a small cross-sectional area and long magnet along 
easy axis is better for having high surface field density and magnetic field 
gradient. Figure 58 shows the effect of the cross-sectional area and length of 
magnet. Contrary to the magnetic field and field gradient trend at the surface 
of the magnet, it can be seen that at a distance of 2 mm (Figure 58a, Figure 
58b) and 6 mm (Figure 58c, Figure 58d) from the magnet surface, the effect 
of cross-sectional area is different from that observed at the surface in Figure 
57. The peak magnetic field and magnetic field gradient are shifting as we 
move away from the magnet surface. This is because the increase in cross 
sectional area results in steep reduction of magnetic field gradient at the 
surface, thus the reduction in magnetic field in the z direction is lower. This 
results in higher values of magnetic field and field gradients at a distance. 
Thus, as the distance from the magnet surface increases the peak is observed 
towards large cross sectional area magnets. The initial self-actuation cycle is 
not accompanied by the inertia force which drives the subsequent oscilla- 
tions. For achieving the first cycle, the magnetic force needs to be large 
enough at the initial position of the cantilever tip. This means the magnetic 
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field for magnetization of film and magnetic field gradient for magnetic force 
need to be large at a distance (at the initial position of the cantilever tip), 
favoring a large cross-sectional area of the magnet. But this would result in a 
lower surface magnetic field and magnetic field gradient, causing a reduction 
in the generated power output as it is only dependent on magnetic field 
gradient and velocity. The effect of cross-sectional area and length of the 
magnet along easy axis is a fundamental aspect in designing a TMG, as the 
cross-sectional area is determined by the size of the film which needs a large 
aspect ratio to benefit from the large surface to volume ratio. Thus, compared 
to the thickness of the TM film, the area of the film must be large, which 
means the magnet also needs to have a large enough cross-sectional area to 
achieve optimum contact. 
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Figure 58 LEM simulation of the effect of dimensions of magnet, (a) magnetic field at a 
distance (D) of 2 mm from magnet surface, (b) magnetic field gradient at D = 
2 mm, (c) magnetic field at D = 6 mm, (d) magnetic field gradient at D = 
6mm. 
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The length of the magnet along easy axis needs to be comparable to the 
cross-sectional dimensions of the magnet else it would result in the loss of 
benefits of miniature dimensions of the TMG. Considering these factors, it is 
possible to find an optimum dimension of the magnet that is guided by the 
TM film properties, operating temperature, and desired dimensions of the 
final TMG in most cases. 


5.2 Effect of design parameters 


To understand the effects of various parameters, the LEM simulation is 
utilized and as a reference TMG sample operating at 170 °C Tsource is chosen. 
The TMG has 10 um Ni-Mn-Ga film as active material. The design of the 
device is taken as the same as presented in section 4.4. By default, Tamb is 
22 °C, Rı is 400 Ohms, number of turns in coil is 400, thermal resistance of 
the bonding layer is 10K/W and contact heat transfer coefficient is 
8.4kW/m°K unless it is stated otherwise. In each subsection from 5.2.1 to 
5.2.5 one of the parameter is varied while keeping all other parameters 
constant to rule out effects of other parameters. The goal is to describe the 
effects of change in each parameter and to check if there is an optimum value 
of parameter where the device performance is optimum in terms of power 
output. The results would serve as a guideline for optimizing the design of 
the TMG as critical parameters that affect the performance can be tuned for 
maximum power output. The trends of the results would be similar for other 
TMG’s with similar design even though the absolute values would be differ- 
ent. 


5.2.1 Effect of heat intake 


Figure 59 shows the effect of contact heat transfer coefficient (hr) on the 
performance of the TMG like stroke, frequency and power. Stroke and 
frequency increases with an increase in the hr which in turn causes the in- 
crease in the power output. 


“Section 5.2.1 has been published in similar form by the author in Joseph et al.”! 
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Thus, there is no optimum value of contact heat transfer coefficient that could 
be derived from this result. But when the hr is below the critical value of 
8 kW/m’K, the device shows a large reduction in performance driven by a 
large reduction in heat intake. 
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Figure 59 LEM simulation of the effect of Contact heat transfer coefficient (hr) on the perfor- 
mance of the TMG with 10 um Ni-Mn-Ga film, (a) stroke of the cantilever de- 
flection, (b) frequency of oscillation of cantilever, (c) power output from the 
pick-up coil. Legend shown is Tsource. Data pre-published in Joseph et al.”! 


As the hfis reduced below 8 kW/m’K, the stroke decreases, and finally, the 
device enters a secondary mode of operation with very low stroke and high 
frequency which is not suitable for energy harvesting. The performance is 
heavily affected by the reduced value of hr at higher Tsource of 160 °C and 
above compared to lower Tsource of 150 °C and below. Figure 60 shows that 
both the heat intake and AT are reducing with the decrease in hr. When the hr 
falls below the critical value of 8kW/m?K, the decrease in heat intake and 
temperature change ( AT) is accelerated. But the AT is still above the critical 
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limit of 6°C, thus, the device sustains resonant self-actuation, more infor- 
mation about the critical limit is given in section 5.2.2. At lower Tsource Of 
150 °C and below, the effect of hr is not so pronounced as the decrease in 
heat intake and temperature is more gradual. The convection losses tend to 
increase when the frequency 
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Figure 60 LEM simulation of the TMG with 10 um Ni-Mn-Ga film, (a) Heat intake into film, 
and (b) temperature change (AT) within the film. Legend shown is Tsource. Data 
pre-published in Joseph et al.”! 


increases as the convective heat transfer is dependent on the velocity of the 
cantilever tip. But when the TMG enters the secondary mode of operation 
with a small stroke and high frequency, the convective losses are still reduc- 
ing as the stroke is not enough to cause an increase in the losses. The cantile- 
ver tip stays close to the magnet. When the Tsource is increasing, heat intake, 
AT, and heat losses tend to increase. This is due to the higher heat flow as the 
result of the higher temperature difference between the TM film and the 
thermal heat source (in this case magnet). Heat intake at contact has major 
effects on the performance of the TMG as it determines the AT of the TM 
film. There is a minimum AT required in the TMG which is determined by 
parameters like material properties and device configuration. Achieving this 
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critical AT within a very short interval to keep up with the mechanical oscil- 
lation frequency is crucial to sustaining resonant self-actuation. 
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Figure 61 LEM simulation of the TMG with 10 um Ni-Mn-Ga film, (a) conduction losses, (b) 
convection losses. Legend shown is Tsource. Part of data pre-published in Joseph 
et al.” 


5.2.2 Effect of heat dissipation 


The thermal resistance of the bonding layer (R») used to bond the TM film 
onto the tip of the cantilever beam has a major effect on thermal dissipation 
by conduction from the TM film to the cantilever. An optimum R» would 
allow a thermal boundary letting the TM film achieve high AT when in 
contact and then would allow it to transfer all the heat into the cantilever 
beam over the cooling cycle. Figure 62 shows the effect of Rb on the perfor- 
mance of the TMG. The performance of the TMG increases when the R» of 
the bonding layer is increased. Both the stroke and frequency are increasing 
for the TMG and consequently, resulting in a higher power output. 


‘Section 5.2.2 has been published in similar form by the author in Joseph et al.” 
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Figure 62 LEM simulation of Rs on the performance of a TMG with 10 um Ni-Mn-Ga film, (a) 
stroke, (b) frequency, (c) power output. Legend shown is Tsource. Data pre- 
published Joseph et al. *! 


When the value of Rp is increased from 10 K/W to 45 K/W, because the film 
retains more heat without cooling much. If the AT in the TM film is below 
the critical value, a change in magnetic force is not enough to sustain contin- 
uous resonant self-actuation. In the present case, this limit follows from 
Figure 63 which shows the break in operation of the TMG at Tsource of 160 °C 
and above to be around 6 °C. The AT in TM film is directly dependent on 
heat intake and heat dissipation. When the Rp is below 5 K/W, the TM film 
and cantilever become thermally coupled, and at higher Tsource values, the 
temperature of the cantilever increases and causes the device to stop opera- 
tion due to the drop in the TM film’s heat dissipation capability. This causes 
the overall AT in the TM film to drop below the critical temperature limit and 
the TMG drops out of resonant self-actuation Figure 63 and Figure 64 show 
the heat intake, AT, and losses in the TMG respectively. 
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Figure 63 LEM simulation of the effect of Rẹ on the TMG with 10 um Ni-Mn-Ga film, (a) heat 
intake at contact, and (b) temperature change within the TM film. Legend 
shown is Tsource. Data pre-published in Joseph et al.”! 


a) 200 


150 


Conduction loss (mW) 


1004 


b) 40 
S a ae 
E 
| EEE 
© N zen 
= 
Ss P 
3 204 = 
3 en 
> 
= 
[e] 
S ai DE 
A pn 


0 20 40 
Thermal resistance - R, (K/W) 


60 


0 20 40 60 
Thermal resistance - R, (K/W) 


Figure 64 LEM simulation of the effect of Rp on heat dissipation in the TMG with 10 um Ni- 
Mn-Ga film, (a) conduction losses, (b) convection losses. Legend shown is 
Tsource. Part of data pre-published in Joseph et al.?! 
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The heat intake and the AT does drop when the R» reaches critical value of 
10 K/W thermal resistance of bonding layer regardless of the Tsource. The 
heat intake and heat dissipation by conduction are almost identical in trend. A 
large value of R» above the threshold would result in the TM film not being 
able to cool within the cooling cycle as can be observed with the reduction in 
heat intake and heat dissipation beyond 20 K/W. Too low R» would thermal- 
ly couple TM film and the cantilever beam resulting in very low AT within 
the TM film. Such a drop in AT below the critical limit of AT would result in 
the failure of the TMG to operate in resonant self-actuation. 


5.2.3 Effect of cantilever stiffness 


Stiffness determines the magnetic force required for the first cycle of self- 
actuation when the inertia force is missing like in the dynamic operation. The 
stiffness also determines the natural frequency of oscillation of the beam 
which needs to be in balance with the thermal cycle to achieve resonant self- 
actuation. The stiffness of the cantilever beam can be modified by either 
increasing the thickness or width of the cantilever beam. The stiffness in- 
creases when the length of the beam is reduced. 


The moment of inertia (J) of a cantilever beam is given by the equation Eq ( 
23 ). 


Eq (23 ) 


Where b is the width of the cantilever beam and h is the thickness of the 
cantilever beam. 


The stiffness of the cantilever beam (k) is given by the equation Eq ( 24 ). 
k= 3E1/72 


Eq (24) 
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Where | is the length of the cantilever beam. Thus, changing the thickness of 
the beam by a few micrometers would result in a very large increase in 
stiffness. Changing the length of the cantilever beam also increases the 
stiffness with a large effect. Thus, it is possible to choose the thickness and 
width of the cantilever beam in the range of the optimum stiffness of the 
cantilever beam and use the length of the cantilever beam to fine-tune the 
stiffness for optimum performance. This is done to keep the stroke of the 
cantilever beam in the range of the magnetic field and field gradient of the 
chosen magnet. Figure 65 shows the effect of cantilever beam stiffness on 
the performance of the device. 
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Figure 65 LEM simulation of the effect of cantilever stiffness on the performance of the TMG 
with 10 um Ni-Mn-Ga film, (a) stroke of deflection, (b) frequency of opera- 
tion, (c) power output. Legend shown is Tsource. 


Increasing stiffness by decreasing the length until the stroke and bending 
angle are in a reasonable range (in which the bending angle does not cause 
fatigue failure in cantilever) in relation to the dimension of the cantilever. 
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Once the stroke is almost half of the cantilever, the stiffness of the cantilever 
beam is increased so that the length of the cantilever beam can be increased 
to keep the ratio of stroke and length, as this affects the bending angle. The 
performance of the TMG increases when the stiffness of the cantilever beam 
is increased, this continues until the optimum value of stiffness is reached. If 
the stiffness is increased beyond this limit, the device drops out of resonant 
self-actuation at higher Tsource values. Like the effects of hr and Rp, the effect 
of stiffness is more pronounced at higher Tsource values. At lower Tsource 
values, the increase in performance is also not pronounced, and the change in 
performance is observed only at higher stiffness closer to optimum and the 
device is not producing any reasonable power output when the stiffness is too 
low. Figure 66 shows the heat intake and AT in the TM film in the TMG. 
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Figure 66 LEM simulation of the effect of cantilever stiffness on (a) heat intake and (b) AT 
within the TM film of the TMG with 10 um Ni-Mn-Ga film. Legend shown is 


Tsource. 


Unlike the previous cases, the increased heat intake is not the reason for 
increased performance in the case of changing stiffness. The increase in 
stiffness results in a decrease in heat intake as lesser AT is required for the 
device to leave the surface of the magnet due to the higher reset force. The 
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AT is only increasing when the stiffness is increased beyond the critical limit 
after which the resonant self-actuation is lost. Figure 67 illustrates the effect 
of stiffness on the losses from the TM film in the TMG. The conduction 
losses follow an identical trend as the heat intake with slightly reduced values 
due to convection losses. The convection losses are also decreasing even 
though the device is operating at a higher frequency when the stiffness is 
increased. The increased frequency causes an increase in the convective heat 
transfer coefficient, however, the overall heat intake and thermal losses 
required for sustaining resonant self-actuation are reduced due to the increase 
in stiffness causing a decrease in convection. 
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Figure 67 LEM simulation of the effect of stiffness, (a) conduction losses, (b) convection 
losses. Legend shown is Tsource. 


5.2.4 Effect of pick-up coil 


The lateral size of the pick-up coil is fixed by the width of the beam and the 
cross-sectional area of the magnet. Having the pick-up coil larger than the 
cantilever beam is not ideal as it would cause additional undesirable air- 
damping. Having a larger pick-up coil size than the cross-sectional area of the 
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magnet would result in having a pick-up coil at low concentration of magnet- 
ic field lines and low magnetic field gradient. This would still generate 
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Figure 68 LEM simulation of effect of pick-up coil on performance of the TMG with 10 um 
Ni-Mn-Ga film, (a) stroke, (b) frequency, and (c) power output. Legend shown 
is Tsource. 


energy, but the increase in mass of the pick-up coil would cause a larger 
decrease in power output than the increase due to the additional size. For all 
devices in this thesis work, the size of the pick-up coil is kept at 2 mm x 
2 mm x 0.3 mm. Using LEM simulations, it is possible to analyze the effect 
of increasing pick-up coil turns (also the added mass) in the TMG. Figure 68 
shows the effect of a number of turns in pick-up coil on the performance of 
the TMG. Increasing the number of pick-up coil turns results in an increase in 
performance until 400 turns and then the performance of the device decreas- 
es. This effect is due to the additional electromagnetic damping and added 
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mass of the pick-up coil. Changing the stiffness, R», etc. can result in a 
different optimum number of pick-up coil turns. The Tsource range also af 
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Figure 69 LEM simulation of the effect of pick-up coil on (a) heat intake and (b) temperature 
change within the material of the TMG with 10 um Ni-Mn-Ga film. Legend 
shown is Tsource- 


fects the optimum pick-up coil turns in the TMG. Thus, the number of pick- 
up coil turns in a TMG should be optimized based on the operating Tsource 
range, AT, TM film properties, operating load and Tamb conditions. The 
frequency of the TMG is also dependent on the pick-up coil turns, but the 
effect is also affected largely by the Tsource. Figure 69 shows the effect of the 
number of pick-up coil turns on the heat intake and AT within the TM film. 
Although, the temperature of the pick-up coil does increase during operation 
of the TMG, due to its large thermal mass, the temperature does not change 
much after saturation within a few cycles of operation. Thus, thermally the 
addition of the pick-up coil turns does not affect the TMG, but the heat 
transfer is affected indirectly through the change in natural frequency of the 
TMG due to added mass and the added electromagnetic damping. The AT is 
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Figure 70 LEM simulation of the effect of pick-up coil turns on heat dissipation of the TMG 
with 10 um Ni-Mn-Ga film, (a) conduction losses, (b) convection losses. Leg- 
end shown is Tsource. 


increasing with Tsource, but also with the increase in the number of pick-up 
coil turns. The large peak of AT at 300 turns for higher Tsource is a result of the 
TMG operating in non-resonant self-actuation mode for Tsource at 150 °C and 
above. At these points, the device is having a low frequency, stroke, and heat 
intake and losses. The stroke of the device and heat intake of the device 
slightly increase with the increase in pick-up coil turns, but the frequency of 
the device reduces due to the mass which has a higher effect on the perfor- 
mance. Figure 70 demonstrates the dissipation losses in the TMG with 
respect to the changing number of pick-up coil turns. The convection heat 
losses seem to increase with the increased pick-up coil turns even when the 
frequency of the device is low because of the increase in heat intake. Both 
losses and heat intake increase when the Tsource is increased, which is unaf- 
fected by the change in pick-up coil turns. 
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5.2.5 Effect of Curie temperature 


If the magnetization curve of the material is kept the same and only the 
transformation temperature is changed, the minimum operation Tsource chang- 
es accordingly. But consequently, the performance of the device at the 
minimum Tsource is also affected. By keeping the ambient conditions the same, 
LEM simulations are conducted to understand the effect of Tc on the mini- 
mum Tsource. Figure 71b illustrates how the minimum Tsource is reducing with 
the reduction in Tc. To keep consistency the design of the device is kept 
constant. The critical Tsource is seen at 38 °C as the device goes out of reso- 
nant self-actuation. The simulation is done from 38 °C to 98 °C of T. with a 
10 °C step. Figure 71a shows the shifted magnetization curves used, and 
only the values of magnetization values for a magnetic field of 0.05 T are 
shown. 
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Figure 71 (a) LEM simulated values of magnetization vs temperature characteristics shifted for 
decreasing Curie temperature (Te), (b) LEM simulation of minimum Tsource for 
different values of Te. Part of data pre-published in Joseph et al.!?® 


Only the Tc is shifted by 10°C, keeping all the other trends the same. The 
minimum Tsource is also shifting along with the shift in Tc. The change in 
minimum Tsource reduces linearly until Te = 68 °C and then the slope of drop 
changes slightly. One of the advantages of Heusler alloy is the ability to tune 
the Te of the alloy by changing the composition of the alloy and its heat 


®Figure 71b has been published in similar form by the author in Joseph et al. "° 
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treatment properties. The minimum Tsource can be further adjusted by chang- 
ing the magnetization values and design ofthe TMG. 


Having a high magnetization value is key to achieving high performance at 
lower temperature differences as it allows for the use of high stiffness of the 
cantilever which would result in higher frequency operation and thus more 
power output. The curve of the minimum Tsource appears to be saturating as 
the Tamb is kept constant at 22 °C. 
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Figure 72 LEM simulation of the effect of Tc on the performance of the TMG with 10 um Ni- 
Mn-Ga film at minimum source temperature, (a) stroke, (b) frequency, (c) 
power output. 


The device goes out of resonant self-actuation for the same reason, as the Tc 
is getting closer to the Tamb. The device would only be able to work at lower 
Curie temperatures if the design is modified to accommodate these changes. 
To finally get the device working at Tsource values of 38 °C and below and to 
increase performance at lower operating temperatures close to Tamb, several 
solutions could be possible, like; 
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e Change the design of the device to optimize it for higher thermal per- 
formance at lower values of Tsource. 


e To change the magnetization properties of the material 
e To lower the Tamb. 


Such change would result in higher performance requiring very low AT 
within TM film. Figure 72 shows the LEM simulations for the performance 
of the TMG at the minimum Tsource for a given Te (Tsource is not same for all 
points but rather the Tsouree in each case is the minimum temperature possi- 
ble). The device power output falls close to zero when the TMG goes out of 
resonant operation. The power output of the device also reduces as the Te is 
reduced without modifying the design of the device. This also illustrates the 
importance of the TMG design and configuration and its dependence on the 
Te with respect to its operation range and ambient conditions. The heat intake 
and AT seem to follow the trend of power output with continuous reduction 
until the device loses resonant self-actuation. 
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Figure 73 LEM simulation of the effect of curie temperature (Te) on the TMG with 10 um Ni- 
Mn-Ga film (a) heat intake and (b) temperature change in the TM film. 
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The heat intake and AT at T. of 48 °C are only 35 mW and 2.5 °C respective- 
ly as shown in Figure 73. The thermal losses from the TM film also decrease 
with the decrease in Te and Tsource values. It can be seen in Figure 74a that the 
conduction losses saturate as Te reaches 48 °C and slightly increase as the Te 
is reduced further due to failure to achieve resonant self-actuation. The 
minimum Tsource at which the TMG is able to sustain resonant self-actuation 
for Te of 48 °C is 70 °C as shown in Figure 71b as the AT is at its critical 
limit of 2.5 °C. However, the convection losses do not saturate even at a 
Tsouree Of 65 °C and reach a minimum value of 2.3 mW suggesting the low 
heat dissipation capability, due to the lower frequency of operation and the 
lower heat intake. 
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Figure 74 LEM simulation of effect on loss mechanisms of the TMG with 10 um Ni-Mn-Ga 
film, (a) conduction losses, (b) convection losses. 


5.3 Discussion 
The requirements for resonant self-actuation, which enables effective conver- 


sion of thermal energy into electrical energy, are investigated using LEM 
simulations. The effect of distance and dimensions on the magnetic field and 
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field gradient of a permanent magnet is discussed. The material composition 
and remanent magnetic field are kept constant. It is shown that the magnetic 
field is the maximum at the surface of the magnet, and it reduces exponential- 
ly as the distance from the surface of the magnet is increased. The effect of 
dimensions on the magnetic field and field was presented earlier in the 
section. The magnetic field strength determines the magnetization of the TM 
film in the TMG, while the magnetic field gradient causes the magnetic force 
when a magnetized material is placed in it. Thus, to achieve maximum 
magnetic force, the field and field gradient needs to be maximum. At 2 mm 
both the magnetic field and field gradient are already very low and reach their 
minimum. Thus, the TMG is always designed to operate with a stroke close 
to 2 mm so that the device still achieves a large stroke but at the same time is 
not moving out of the magnet’s peak magnetic field and magnetic field 
gradient. 


The dynamics of heat absorption upon mechanical contact with the heat 
source and heat dissipation via heat conduction and convection via the 
oscillating cantilever are discussed. The heat transfer coefficient (hr) turns out 
to be an important design factor when it comes to heat intake. It is crucial to 
build TM films with a large enough surface-to-volume ratio to support high 
heat flows and quick AT in the TM film. For the best thermal contact, a heat 
source with a smooth surface finish is also necessary. The temperature 
difference between the heat source (permanent magnet) and the TM film 
affects how much heat is transferred. Generally speaking, the power output 
will increase as hr increases. According to the results, the critical value of hris 
around 8 kW/(m°K) in the case of the TMG with 10 um Ni-Mn-Ga film. 
Above this point, consistent operation in resonant self-actuation with large 
stroke and frequency is observed yielding high power output. At large Tsource 
values of 160 °C and above, when hr was reduced below the critical limit, the 
device performance degrades as the AT decreases. This is resulting from 
imbalance between heat intake and dissipation. The system can no longer 
maintain resonant self-actuation below a threshold AT of several degrees, 
depending on the thickness of the TM film. The essential AT in the consid- 
ered case is around 6 °C. 
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The thermal resistance (Rv) of the bonding layer is another important design 
factor. It controls how heat is transferred from the TM film to the cantilever, 
which is the heat sink of the TMG. The thickness of the non-conductive 
bonding layer can be modified to vary this parameter, which produces a 
thermal boundary between the TM film and the cantilever. This results in a 
high enough AT at low heat intake. LEM simulation results show that for all 
Tsource values, there is a critical value of R» of around 10 K/W above which 
stable performance occurs and is associated with high values of oscillation 
stroke, frequency, and power. A strong reduction in AT below the critical 
threshold results from an imbalance between heat intake and dissipation, 
which deteriorates the TMG performance. The device performance will be 
greatly reduced if the R» is lower, and it will finally stop operating in the 
resonant self-actuation mode if AT falls below the critical threshold of 6 °C. 
At Tsource values of 160 °C and above, where resonant self-actuation ceases 
operation below the critical limit, the bonding layer's impact on the device's 
performance is more obvious. 


The stiffness of the cantilever also affects the performance of the cantilever 
and choosing the appropriate stiffness is key in optimizing the performance 
of the TMG and also to achieve resonant self-actuation. If the stiffness is too 
low or too large, the device loses resonant self-actuation mode and results in 
a large reduction in power. If the stiffness is too small, the device stays in 
contact for a longer time with the magnet resulting in excess heat flow into 
the TM film, which in turn heats the cantilever causing a very large AT, but at 
the expense of frequency. The stiffness is adjusted in terms of thickness and 
then fine-tuned by adjusting the length of the cantilever so that the device is 
capable of achieving resonant self-actuation mode with a stroke of around 
2 mm. The number of turns in the pick-up coil affects the performance in two 
ways, one with the change in electromagnetic damping and the other due to 
the change in mass of the pick-up coil. In either case, the number of turns in 
the pick-up coil needs to be optimized to achieve maximum output power. If 
all the other parameters are constant it can be seen that the power output is 
maximized at different pick-up coil turns at different Tsource values. This is 
due to the proper balance between the TMG’s electromagnetic damping, 
pick-up coil mass, and the thermal cycling of the device at various Tsource 
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values. This shows the importance of optimizing the design for operating 
conditions like the Tsource and Tamb to achieve maximum power output. 


The Tc of the device has major effects on the TMG’s minimum Tsource re- 
quirement above which it would be capable of achieving resonant self- 
actuation. If the material properties like magnetization and AM/AT at the 
transition are kept the same, it can be seen that the power output reduces as 
the Tc is reduced. This is due to the inability of the device to maintain the 
same conditions of heat intake and heat dissipation for lower Tsource values 
while the Tamb is remaining constant. Increasing the magnetization allowing 
the TMG to have higher stiffness for the cantilever is one of the solutions to 
overcome this loss in power. Increasing the value of AM/AT is also another 
method by which higher power outputs could be achieved, but this also 
results in a lower range of operation temperature as the device would be 
highly sensitive to Tsouree and would only exhibit resonant self-actuation 
within a narrower temperature range. The findings show that device perfor- 
mance can be further enhanced in terms of heat intake and dissipation, 
leading to an increase in electrical power output while operating in resonant 
self-actuation mode. 


Table3 Modeling parameters used in simulation by LEM of the TMG. Part of data pre- 
published in Joseph et al.?!!78 


Mechanical parameter Unit Value References 
Pick-up coil mass Kg 6x10°° | This work 
Cantilever density Kg/m? 8500 ve 
TM film density Kg/m? 8020 wt 
Bonding layer density Kg/m? 920 a 
Structural and air damping Ns/m 2.81075 This work 
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Impact damping Ns/m 0.1 This work 

Contact stiffness N/m 10000 This work 

Young’s modulus Pa 1x 1014 148 

Electrical & Magnetic Unit Value References 

Parameters 

Remanent magnetic field T 1.07 This work 

Area of pick-up coil m? 2.81 This work 
x 1076 

Load resistance (RL) Q 400 This work 

Thermal Parameters Unit Value References 

Cantilever thermal conductiv- W/mK 109 on 

ity 

TM film thermal conductivity W/mK 23.2 bie 

Bonding layer thermal W/mK 0.225 ” 

conductivity 

Cantilever heat capacity J/kgK 400 se 

TM film heat capacity J/kgK 490 8 

Bonding layer heat capacity J/kgK 2100 8 

Convective heat transfer W/m?K 40-140 This work 

coefficient 
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6 Upscaling of TMG 


Parts of this section has been published in Joseph et al.?®!5° This section 
provides the results of the investigation of upscaling and the resulting effects 
on performance metrics, and coupling effects in parallel operating TMGs. 
The investigation utilized TMGs with Ni-Mn-Ga as active TM material. The 
first device with resonant self-actuation is presented in Gueltig et al.” with a 
high power density of 118 mW/cm?. This is normalized to the active materi- 
al volume, which is consistent with theoretical predictions of up to 300 mW/ 
cm? for thin plates.'® No external actuation or influence is required for 
operation, as this is a self-actuating and self-adjusting device, with the only 
requirement of operation being the thermal energy source. But the absolute 
power output from the device is only 2.4 uW, corresponding to a power per 
footprint of 14.2 uW/cm? and in order to make the TMG attractive for most 
applications it is required to increase the amount of raw power output. There 
are two different ways to increase the power output or to upscale the TMG, 
one is to increase the power output from a single device (power output per 
area), second is to add more devices operating in parallel to combine the 
overall output. Section 6.1 shows the effects of upscaling the material thick- 
ness to achieve higher power output from consumed area. Section 6.2 de- 
scribes the effects of devices operating in parallel close to each other. 


6.1 Upscaling TMG by volume of active material 


In order to upscale the power per footprint of the TMG, the amount of active 
material in the device needs to be increased. The scaling of the TM film is 
done by increasing the thickness, so the lateral area of the TMG is kept 
constant for comparison. 


"Section 6.1 has been published in similar form by the author in Joseph et al.”* 
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6 Upscaling of TMG 


The main goal is to scale the TMG without losing benefits of high efficiency 
and rapid heat transfer found in film-based devices while meeting power 
levels for applications at different scales. Cantilever dimensions are adjusted 
accordingly to maintain resonant self-actuation which is required for main- 
taining the large power output while operating at high frequencies. The TMG 
performance is then analyzed using LEM to describe the thermal processes in 
the TMG. The modeling parameters are summarized in Table 4 (p. 140) 


6.1.1 Design ofthe device 


The Ni-Mn-Ga film thickness of three demonstrators is chosen to be 10, 20, 
and 40 um, respectively. Due to limitations in sputtering duration and associ- 
ated film homogeneity, film thicknesses are yet restricted to a maximum of 
around 10 um. Therefore, as shown in Figure 75, the thickness is sequential- 
ly raised by stacking identical Ni-Mn-Ga films of 10 um thickness above one 
another. These films have the same lateral dimensions and magnetic proper- 
ties. The film stack is fixed utilizing conductive bonding layers between the 
films and a non-conductive bonding layer in between the film stack and 
cantilever. The performance of the TM generator is significantly impacted by 
the increase in heat and inertia mass during upscaling. 


b) 


a) 
heated permanent 
magnet 


MSMA film stack 
O cantilever 
] pick-up coil \ 
o substrate 


Figure 75 Design of film stack to increase the thickness (a) schematic of the TMG with 
oscillating cantilever and Ni-Mn-Ga film stack, (b) section of the cantilever 
front showing the film stack with multiple layers of Ni-Mn-Ga films. Adapted 
from Joseph et al.” 
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6.1 Upscaling TMG by volume of active material 


To sustain the conditions of rapid heat transfer and resonant self-actuation, a 
variety of challenges must be resolved. In particular, heat intake during 
contact between the film stack and the heat source (permanent magnet) as 
well as heat dissipation, both of which depend on the dynamics of the canti- 
lever motion, are important. This also includes fabrication technology, 
optimization of cantilever geometry affecting footprint, electrical load re- 
sistance (Rx), and other factors. 


6.1.2 Fabrication of TMG 


Having thermally conductive bonding layers within the TM film layer in- 
creases the average thermal mass and the thermal resistance of the layer. 
Thickness of the TM layers is taken as 10 um, 20 um, and 40 um. However, 
the actual thickness of the whole layer is 10 um, 33 um, and 82 um respec- 
tively due to the added size of the bonding layer between. The bonding layer 
thickness is hard to control due to the size of the conductive particles (silver 
in this case) within the bonding layer. To achieve the minimum thickness of 
the bonding layer between the films and keep the thermal resistance of the 
bonding layer to a minimum, specific steps are taken. The fabrication of the 
TMG is then completed in the same manner as for the TMG with 10 um Ni- 
Mn-Ga film. 


6.1.3 Effects of upscaling Ni-Mn-Ga thickness 


The performance of the TMG increases when the thickness of the TM film is 
increased from 10 um to 40 um. This is due to the additional volume which 
results in the excess magnetization of the TM film which corresponds to 
higher magnetic forces. This allows for modifying the design of the TMG 
with higher stiffness of the cantilever to compensate for the additional mag- 
netic force. The stroke and frequency also increase with the increase in 
thickness of the TM film. 
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Figure 76 Effect of film thickness on the performance of the device, (a) stroke, (b) frequency, 
(c) power output. Data pre-published in Joseph et al.” 


Figure 76 shows the general trend of performance with increasing Tsource 
when the Ni-Mn-Ga film thickness is scaled. It can be observed that the 
power output increases as the thickness of the TM film is increased. For the 
TMG with 10 um TM film, the maximum generated power output is around 
3.1 pW. This value increases to 4 uW for the TMG with 20 um TM film and 
finally reaches a maximum value of 5.4 uW for the TMG with 40 um TM 
film. This is the result of increased stroke and frequency. The thickness of the 
TM film cannot be increased any further as the overall thickness of the stack 
has already reached 82 um for a TM film thickness of 40 um. 


LEM simulations showed that increasing the TM stack thickness by adding 
any additional film would result in too much thermal resistance within the 
stack and effectively result in reduced performance. But having a homoge- 
nous material without any thermal resistance of the bonding layer would 
enable going further than 40 um. The effect of load resistance and source 
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6.1 Upscaling TMG by volume of active material 


temperature on the TMG performance for varying Ni-Mn-Ga film thickness 
is given in Figure A 2 to Figure A 4. Performance results are summarized in 
Table 5 (p.144). From Table 5 it is clear that having a homogenous material 
layer instead of TM stack without thermal resistive layer results in a minor 
performance increase due to added benefit of higher thermal performance. 
The heat intake into the TM film stack on contact is shown in Figure 77a. 


| —— 401m film 
2801 —"— 20um film 
—=— 10um film 


Heat input (mW) 
N 
s 
Temperature change (°C) 


j 
1204 | ee Te 
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120 130 140 150 160 170 180 120 130 140 150 160 170 180 
Temperature (°C) Temperature (°C) 


Figure 77 LEM simulation of the effect of film thickness on (a) heat intake and (b) the average 
temperature change within the TM film. Data pre-published in Joseph et al.”® 


The increase in heat intake when the thickness is increased is due to the 
increase in heat transfer coefficient (hr) as a result of the higher contact force 
from additional magnetization (more material) and impact force owing to 
inertia. This results in a higher frequency of operation as the AT required in 
the material to achieve the magnetization change is now much less as the 
volume of the material is increased. Thus, the AT within the TM film stack 
reduces during operation when the thickness of the film is increased. The 
losses by conduction and convection are also increasing in response to the 
increase in the thickness of the TM film stack. But the increase in conduction 
losses can be seen stagnating as the thickness is increased to 40 um. This is 
due to the thermal resistance of bonding layers within the TM film stack. 


125 


6 Upscaling of TMG 


a) 240 b) 
—— 10um film 
220)“ 20um film 604 
—— 40um film 
=> 2004 > 504 
3 = 
E£ E£ 
u 1804 n 
E 8 40- 
wn N 
2 160 2 
o oO 
3 2 304 
9 1404 3 
no} > 
= 
je] je] 
© 1204 O 204 
4 
1004 
104 " 
80 


120 130 140 150 160 170 180 420 130 140 150 160 170 ` 180 
Temperature (°C) Temperature (°C) 


Figure 78 LEM simulation of the effect of thickness on losses, (a) conduction losses, (b) 
convection losses. Data pre-published in Joseph et al.” 


Since TMGs with 20 um and 40 um, TM film stacks have multiple layers of 
TM film, the temperature of the individual layers is different as they are 
separated by the conductive bonding layer which still has much higher 
thermal resistance compared to the TM film. The temperature of the top layer 
is always higher and the bottom layer is the lowest. But this gradient of 
temperature within the stack is slightly increasing with the increase in Tsource. 
The AT within the individual TM film layers can be observed in Figure 79a. 
Table A 1 provides the maximum temperature (Tmax) and minimum tempera- 
ture (Tmin) that provides the value of AT in the TM film. For the TMG with a 
40 um TM film stack, the AT distribution is given in Figure 79b. The maxi- 
mum AT occurs in the top film at a Tsource of 170 °C and the AT within the 
layers as we move from second to fourth shows only slight reduction. Maxi- 
mum AT at Tsource of 170 °C is 2.9 °C in the top layer, while the bottom layer 
shows only a AT of 2 °C. Figure 79 shows that the top layer experiences the 
maximum change in magnetization during operation and contributes more to 
the change in magnetic force change compared to the layers below. The 
bottom film at 130 °C only experiences 1.6 °C of change and still sustains the 
resonant self-actuation. 
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Figure 79 LEM simulation of (a) temperature change within different film layers in a TM film 
stack (b) temperature profile within the TM film stack in a TMG with 4 films 
(40 um effective thickness). Part of data pre-published in Joseph et al.” 


Even in the case of homogenous TM film, there would be a temperature 
gradient from top to bottom, but this would still be much lower than when the 
TMG is made with multiple layers bonded together. One of the advantages of 
upscaling the TM film layer is that the footprint occupied by the TMG 
reduces to as result of device design. 


This is due to the reduction in the total length of the cantilever to increase the 
stiffness of the TMG to match the increased magnetic force and modified 
thermal cycle. The lateral area occupied by the TMG is reducing with an 
increase in thickness of the TM layer. The efficiency of the TMG in terms of 
magnetic energy converted from the thermal energy is staying almost con- 
stant even when the TMG is scaled. Figure 80 shows that the footprint is 
reduced to 10.6 mm? while the power increased to 5.3 uW. Here the im- 
portant result is upscaling the power output per footprint by about 3.4 times 
to a high value of 50 uW/cm?, when compared to the reference device of the 
TMG with 5 um Ni-Mn-Ga film for the same Tsource of 170 °C. 
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Figure 80 Effect of scaling on power, footprint, and efficiency. Adapted from Joseph et al.” 


6.2 Upscaling TMG by array of parallel devices 


Upscaling the power output of the TMG by adding more devices operating in 
parallel is yet another method by which higher power outputs could be 
achieved from the TMG. The operation of the TMG could either be in syn- 
chronization or without synchronization with each other. However, being 
oscillating devices, it can be assumed that TMGs would have some form of 
cross-coupling between them while operating in close proximity to one 
another, as they have pick-up coils and TM films that change their magneti- 
zation. Both of these might cause an electromagnetic coupling. Moving air 
between the devices might result in mechanical or thermal coupling of the 
device. Understanding such coupling effects between the TMG devices 
operating at close distances is crucial to realize upscaling of the TMG as an 
array of multiple parallel devices. The power output from a single TMG 
device would be limited and in order to increase the power output further 
while still keeping the benefits of resonant self-actuation, rapid heat transfer, 
and thin TM film architecture, it is inevitable to have multiple TMGs operat- 
ing under the same heat source. 


128 


6.2 Upscaling TMG by array of parallel devices 


6.2.1 Layout oftwo device array 


Two devices are operated under the same magnet with enlarged dimensions 
3 mm x 30 mm x 8 mm. The distance between the devices (D) is varied from 
6 mm to 0.3 mm. The lower limit of separation between the devices is due to 
the physical limitation to accommodate the devices without mechanical 
contact between the two devices (due to the wires coming out of the pickup- 
coil. The design of the device is the same as in section 4.4. TMG’s with a 
brass cantilever of 7.2 mm x 2 mm lateral dimensions and thickness of 20 um 
is chosen. TMGs are prepared by keeping close tolerances as the devices 
need to be as identical as possible. 


heated permanent 
magnet 


SMA film A 


cantilever A 
pick-up coil A 


SMA film B 
cantilever B 
pick-up coil B 


substrate 


Figure 81 Schematic of two TMG’s operating in parallel under a long magnet. Adapted from 


Joseph et al.!°° 


One of the TMGs is placed under the magnet and moved across the surface 
while measuring power output, frequency, and stroke to make sure the effects 
of magnetic misalignment and precision position stage misalignments are 
alleviated. Once the magnet and position stages are all aligned, the second 
TMG is placed at a distance of 5 mm. The separation between the TMGs is 
varied in steps. Figure 81 shows the schematic of the TMG. 


‘Section 6.2.1 has been published in similar form by the author in Joseph et al.'*° 
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At each position, ten measurements of a duration of ten seconds are taken 
after giving a few minutes for the device to stabilize after moving the device 
closer. As explained in detail in section 5.1, the change in the dimension of 
the permanent magnet affects the magnetic field and magnetic field gradient 
of the magnet. Therefore, absolute numbers can hardly be compared between 
the results in previous sections 4.4, 5 and 6.1. Figure 82 illustrates the 
variation of magnetic field and magnetic field gradient with respect to dis- 
tance from magnet surface. Both magnets are made of SmCo to be operation- 
al at higher temperatures. The 3 mm long magnet has a remanent magnetic 
field of 1.07 T and the 30 mm long magnet has a remanent magnetic field of 
1.08 T due to a different composition of the magnetic material. 
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Figure 82 Comparison of (a) magnetic field and (b) magnetic field gradient of 3 mm long and 
30 mm long magnet based on LEM simulations. 


The difference in magnetic field due to the slightly higher remanent magnetic 
field is really low as the change is only 0.1 % of the value. The 30 mm long 
magnet has higher magnetic field and slightly higher magnetic field gradient 
at a longer distance from the surface of the magnet compared to the magnet 
with 3mm length used in section 4.4. At the surface of the 3 mm long mag- 
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net, the magnetic field intensity is around 525 mT, while the 30 mm long 
magnet can only generate a magnetic field intensity of 480 mT at the surface 
of the magnet. 


Closer to the magnet, the magnetic field and field gradients are low on the 
longer magnet. This results in lower magnetization of TM film and conse- 
quently lower magnetic forces at the tip of the cantilever beam. Thus, the 
operation frequency will be lower, and the generated power will be lower as 
well due to the lower magnetic field gradient experienced by the pick-up coil. 
The increase in the remanent magnetic field by changing the composition is 
not enough to counter the reduction in the magnetic field and magnetic field 
gradient of the 30 mm long magnet compared to the 3 mm long magnet. The 
TMG with 10 um thick Ni-Mn-Ga film is generating around 1.3 uW under 
the large magnet while it is placed as a single operating device under the 
magnet. At 5 mm away from each other, the devices are assumed to be at the 
same performance level as the distance between them is too large to cause 
any significant coupling effects. 


6.2.2 Performance results of two device array 


When two identical devices with very tight tolerances are operating near each 
other, the performance of the device is affected by the device nearby. Figure 
83 shows the effect of separation between the TMG devices for a Tsource 
temperature of 170 °C and 140 °C. For high source temperature of 170 °C, 
the device shows reduction in performance as the separation is reduced below 
1 mm. The device shows breaks in operation when the separation is further 
reduced to 300 um. On the other hand, at lower source temperature of 
140 °C, the TMG devices seems not to show any coupling effects even at 
close separation of 300 um. It is also visible that the power output even 
increases slightly as the devices are brought closer to each other. Figure 84 
illustrates the behavior of the TMGs when the separation between them is 
varied at different Tsource values from 170 °C to 110 °C. At lower Tsource when 
the frequency stroke and power output are already low, the separation does 
not seem to have any effect. When the temperature of the heat 


Section 6.2.2 has been published in similar form by the author in Joseph et al.!5° 
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Figure 83 Effect of separation on the TMG at Tsource of 170 °C and 140 °C. (a) Stroke of the 
deflection of the TMG, (b) Frequency of operation of the TMG, (c) power out- 
put of the TMG. Part of data pre-published in Joseph et al.'*° 


source is increased, the device starts to show some coupling effects when the 
separation is small. When the Tsource is 150 °C and above, the power seems to 
drop continuously as the separation is reduced below 1 mm. However, the 
devices show coupling effects when separation below 1.5mm. At 140 °C, 
TMG was operating with low power output even at 6 mm separation suggest- 
ing the device was operating closer to the magnet than its optimal operating 
distance where it would generate maximum output. Thus, this shows that if 
the device is operated at suboptimal conditions where the power output is 
reduced, the TMG could be operated at very close distances. 


This also explains the slight increase in power output that was observed as 
the separation was reduced suggesting there is still potential to optimize the 
power output. It is clear that at higher source temperatures, the devices 
become highly unstable and continuously stopoperation resulting in large 
fluctuation in power output. Further examination of Figure 84 reveals that the 
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average stroke is showing reduced values below 1.5 mm, while the frequency 
shows an increase at high source temperatures above 150 °C. This slight peak 
in frequency of operation results in an increase in output of the device before 
it stops operation suggests a thermal build-up in the Ni-Mn-Ga films when 
the devices have a small separation between them. 
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Figure 84 Effect of separation on performance of TMGs. (a) Stroke of deflection, (b) frequen- 
cy of operation, (c) power output. 


Further analysis is only possible by checking individual performance as the 
average values are only showing the device trend. To understand the coupling 
effects that cause break in operation, it is necessary to examine the time 
correlated performance of the TMG devices. Figure 85 shows the output 
generated by the TMG, measured by an Oscilloscope (Keysight 
DSOX1204A). The oscilloscope is connected parallel to the load resistor so 
as to measure the voltage across the resistor. The measurement is taken for a 
long duration of 180 seconds. The device stops operation due to the increased 
temperature of the Ni-Mn-Ga film and restarts itself as soon as the film is 
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cooled enough to resume operation. This is clear from the increase in voltage 
before damped oscillations appear. But the interval of this intermittent break 
is not regular showing that the mutual coupling effects of the TMGs are 
random. The mutual effects could be visualized by plotting the output of the 
TMGs together as shown in Figure 86. The output of the two devices TMG 
A and TMG B denoting two of the TMGs operating in parallel is mutually 
coupled at small separations. 
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Figure 85 Output performance of the TMG (a) at a separation of 1200 um between devices and 
(b) 400 um between devices. Data pre-published in Joseph et al.!5° 


The output of the device TMG B increases when TMG A stops operation and 
shows adaptation to a lower power output operation point as soon as device 
A is back in resonant operation. However, this mutual correlation is not 
present in all the data, which suggests that such a correlation can occur but is 
not required when the devices are operating in parallel. The mutual correla- 
tion is occurring due to the ability of TMG A to cool itself better as the 
device nearby stops operation. This shows the ability of the TMG to adapt its 
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performance according to ambient conditions. In cases where the correlation 
is absent, it can be assumed that the devices are not coupled to each other and 
thus a momentary halt in the operation of the nearby device is not enough to 
cause an increase in the thermal performance of the other device. From 
Figure 85, it can be assumed that there is no mechanical effect due to the air 
movement between the devices other than the indirect effect on the TMG 
thermal dissipation. This is because in case of mechanical coupling between 
devices due to air movement or electromagnetic coupling between coils, the 
TMG would have shown reduced frequency and voltage output before going 
out of resonance. But in the present case the TMG shows an increase in 
frequency and voltage before going out of resonance. This increase in fre- 
quency and voltage is the result of heat buildup as the device is not cooling 
enough, hence the TMG requires less time to heat up, but eventually this 
causes a momentary halt in the operation. 


Voltage (mV) 


Voltage (mV) 


T x T T 
60 80 100 120 140 160 
Time (s) 
Figure 86 Simultaneous performance comparison of the TMG devices. The Tsource is maintained 
at 170 °C and separation is 1000 um, (a) voltage output across load resistance 


for TMG sample A, (b) voltage across load resistance for TMG sample B. Data 
pre-published in Joseph et al.!°° 
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The TMG continues operation once the TM film is cooled enough to be 
attracted to the magnet so that the cycle could be restarted. Thus, it is con- 
cluded that only thermal coupling effects are present without any mechanical 
or electromagnetic coupling effects between the TMG devices operating in 
parallel under a common magnet. 


6.2.3 Layout of three device array 


Having two devices operating in parallel is symmetric with both devices 
experiencing the same ambient conditions. Adding more devices in parallel 
would result in TMG’s experiencing different ambient conditions. The TMGs 
at the ends would have one neighboring device and all in between TMGs 
would have TMGs operating on either side of them. To see the effect of 
having two devices operating on either side of the device, a total of three 
identical devices are prepared. Thereby keeping very tight tolerances are 
placed under the 30 mm long magnet used in the case of two device array 
explained in section 6.2.1. Figure 87 shows the schematic of three TMG 
device array. 


Permanent magnet 
Ni-Mn-Ga film 

Brass (Cu Zn) cantilever 
Substrate 

Pick-up coil 


Figure 87 Parallelization with three devices, with the substrate for the center device rotated by 
90° to allow low separation between devices. 
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The substrate for the center device is rotated by 90° to allow very low separa- 
tion between devices. The overall power output of the center device is re- 
duced to around 0.7 uW. compared to 1.3 uW in the case of devices at the 
end. Like in the case of two device array operating in parallel, the coupling 
effects are present from a Tsource of 150 °C and gets worse as the separation 
between the device is below 1 mm. Figure 88 shows the effect of separation 
on the performance of the center TMG for various Tsource values ranging 
between 140 °C to 160°C. As mentioned earlier the device is operating by 
generating lower power compared to the devices at the side, suggesting even 
at 5 mm the device is having an effect. Reducing the separation from 5 mm to 
1 mm does not affect the performance of the device with power output and 
frequency staying constant. Reducing separation below 1 mm accelerates the 
decrease of power output and causes an increase in the frequency of the 
device due to thermal build-up in the Ni-Mn-Ga film. The average frequency 
increase is different from the case of two device array as in that case the 
average frequency is constant and only an increase in frequency is observed 
locally near the intermittent stop in the TMG operation. 
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Figure 88 Performance of the center device (the TMG with 10 um Ni-Mn-Ga film) in an array 
of three devices as the separation is gradually reduced. The two neighboring 
devices are moved and kept at the same separation. (a) power output of the 
TMG (b) frequency of the TMG. 
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The present case with three device array also exhibits a sudden increase in 
frequency before breaking resonant self-actuation. But on top of it, there is an 
increase in the average frequency of the operation. The power first increases 
and then decreases as the separation is reduced below | mm at higher tem- 
peratures. At 140 °C, the power is increasing slightly as the devices are closer 
with separation below 1 mm. 


6.2.4 Effect of upscaling on coupling effects 


Upscaling the device from having a Ni-Mn-Ga film of 10 um thickness to a 
device with 20 um thick Ni-Mn-Ga provides insight into effects of scaling on 
coupling effects while in parallel operation. It is shown earlier in section 6.1 
that upscaling of the TMG results in increased power and decreased footprint. 
A single TMG device with 20 um Ni-Mn-Ga film under the 30 mm long 
magnet generates 1.8 uW of power compared to the 1.3 uW. But when 
placed in the array, the center device only generates 0.8 uW. When compared 
to the power output of 0.7 uW for the TMG with 10 um Ni-Mn-Ga film, it is 
only an increase of 14.3 %. 
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Figure 89 Performance of center device in the TMG array of three devices with 20 um Ni-Mn- 
Ga film, (a) power output, (b) frequency of operation. 
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At higher Tsouree values the power output of the TMG is decreasing as the 
separation is below 1 mm. But the frequency increases as the device gets 
closer as can be seen from Figure 89. At lower Tsouree values the power 
increases as the separation between the devices is reduced below 1 mm. This 
also suggests thermal imbalance with the TMG requiring additional optimiza- 
tion to operate at full capacity when placed in an array. The highest power 
output is observed at the closest distance of 100 um separation. This is also 
observed in Figure 89 where at 140 °C, the power output increases as the 
separation is reduced. This shows that the design of a TMG with multiple 
devices in parallel would require additional design optimization depending 
on whether the device would be operating in between or at the end of the 
TMG array. 


6.3 Discussion 


The large increase in electrical power for increasing film thickness up to 
40 um while simultaneously decreasing the footprint is a notable result of this 
work. As a result, when increasing film thickness from 10 to 20 um, the 
magnetic power per footprint calculated from the magnetic cycle increases 
from 0.8 mW/cm? to a maximum of 1.1 mW/cm?. A further increase in layer 
thickness to 40 um causes a marginal drop in power density to 1.05 mW/cm?. 
The magnetic power is taken as the enclosed area in the magnetization vs 
magnetic field curve. This curve shows also shows the thermomagnetic cycle 
followed by the TMG, an example of thermomagnetic cycle depicted in 
Figure A 1 for the case of the TMG with 40 um film. This is due to the 
added thermal resistance of the conductive bonding layer between the TM 
film stack (Ni-Mn-Ga) As the heat intake rises further with increasing film 
thickness, this decrease is correlated with a little decline in efficiency. The 
electrical power produced by electromagnetic induction in the pick-up coil, 
however, increases continuously with film thickness and peaks at 5.3 uW ata 
thickness of 40 um for the Ni-Mn-Ga film. 50 uW/cm? is the proportional 
electrical power per footprint. In this instance, only a 3 °C periodic tempera- 
ture change in the film stack is needed. This represents an increase of more 
than 3.4 times when compared to the reference device with a 5 um film 
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thickness from the previous work.” When it comes to length scales in the cm 


range and small temperature differences, this performance rivals that of 


cutting-edge thermoelectric generators. 1! 


Table 4 Model parameters for thermal, electromagnetic, and mechanical components being 
used in LEM simulation of the TMG upscaling. Part of the data pre-published in 
Joseph et al.7°118 


Mechanical Parameters Value Units Reference 
Pick-up coil mass 6 x 107° kg This work 
Cu-Zn density 8500 kg/m? > 

Ni-Mn-Ga density 8020 kg/m? 1# 

Young’s modulus (Cu-Zn) 1x 1011 N/m? 146 

Contact stiffness beam-magnet 10? N/m This Work 
Structural and air damping 1.12 x 1075 Ns/m This Work 
Impact damping 0.1 Ns/m This work 
Length of the Ni-Mn-Ga stack 2 mm This work 
Width of the Ni-Mn-Ga stack 2 mm This work 
Thickness of the Ni-Mn-Ga stack 10-40 um This work 
Length ofthe Cu-Zn cantilever 4.3-5.7 mm This work 
Width of the Cu-Zn cantilever 2 mm This work 
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Thickness ofthe Cu-Zn cantilever 


Electromagnetic Parameters 
Remanent magnetic field strength 
Number of turns of pick-up coil 
Area of pick-up coil 

Electrical resistance (pick-up coil) 


Electrical resistance (load) 


Thermal Parameters 


Range of contact heat transfer 
coefficient 


Range of convective heat transfer 
coefficient 


Thermal conductivity Cu-Zn 
Thermal conductivity Ni-Mn-Ga 


Thermal conductivity of conduc- 
tive bonding layer (BL) 


20 


1.07 


400 


1.96 x 107° 


250 


400 


8400-13300 


150-220 


109 


2.5 


um 


W/m?K 


W/m?K 


W/mK 


W/mK 


W/mK 


This work 


This work 


This work 


This work 


This work 


This work 


This work 


This work 


147 


148 


145 
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Thermal conductivity of non- 0.33 W/mK o 
conductive BL 

Heat capacity Cu-Zn 400 J/KgK = 
Heat capacity of Ni-Mn-Ga 490 J/KgK 14 
Heat capacity of conductive BL 904 J/KgK 145 
Heat capacity of non-conductive 2100 J/KgK 145 

BL 

Area of thermal contact 4 mm? This work 


In order to maintain the resonant self-actuation mode, the strategy to boost 
power depends on tuning the magnetic attraction force, mechanical reset 
force, and heat transfer dynamics. The realized increase in power results from 
the rise in the oscillatory motion’s mechanical energy, which rises with 
increasing mass and frequency as long as resonance conditions are upheld 
and, in turn, sustains magnetic cycles with rising variations in magnetic 
energy. The presented results show that the upscaling strategy is ultimately 
constrained by the heat transfer kinetics for a given ferromagnetic material 
and magnetic field gradient. The scaling dependency of heat input, heat 
dissipation, and the ensuing temperature fluctuations in the Ni-Mn-Ga layers 
are analyzed using LEM for describing the thermal processes in the TMG. 
The dynamic actuation of the cantilever is limited by the thermal processes 
that are result of parameters like TM film thicknesses, load resistance, coil 
parameters, ambient temperature and Tsource values The LEM is also used to 
determine the heat transfer coefficient (hr) values by matching simulation 
results with the experimental results. Table 4 provides the modeling parame- 
ters ofthe LEM simulation model. 


According to LEM simulation results, heat input rises as film thickness 
increases, with conductive heat losses serving as the main mode of heat 
dissipation. As a result, with increasing film thickness from 10 um to 40 um, 
respectively, the corresponding fraction of average retained input heating 
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power increases from roughly 80% to 94%. These fractions only differ 
slightly and exhibit the same thickness dependency at lower Tsource values. 
This finding demonstrates that a significant portion of heat is retained in the 
SMA film stack and contributes to the AT, a necessary condition for achiev- 
ing large power output. But as the film thickness increases, the AT of the film 
stack during oscillation decreases, which has two significant effects. First, the 
effectiveness of increasing film thickness is reduced as AT decreases. Sec- 
ond, the diminishing AT places a severe restriction on the further scaling up 
of film thickness. 


The results of the current study clearly show that a minimum AT is necessary 
for resonant self-actuation, which is dependent on the specifics of material 
properties, operating temperatures, and device design. This limit is between 
two and three degrees Celsius for the TMG with 40 um Ni-Mn-Ga film. 
Table 5 summarizes the key findings of the demonstrator’s performance. The 
attainable AT is no longer sufficient to induce enough change in magnetiza- 
tion AM to support resonant self-actuation when the film thickness is too 
high. The TM film's bonding layers minimize heat transmission, resulting in 
a low temperature gradient (<1 °C) from top to bottom that minimizes the 
contribution of the bottom layers to resonant self-actuation. LEM simulation 
allows to evaluate important temperatures of the TMG demonstrator, which 
are provided in Table A 1. LEM simulation is carried out to examine this 
impact, replacing the stack of four Ni-Mn-Ga films and thermally conductive 
connecting layers with a single homogenous Ni-Mn-Ga layer that is 40 um 
thick, leaving only one non-conductive bonding layer that is 15 um -20 um 
thick. This resulted in a 10% increase in power output which could be 
increased further as the thickness of the actual device was 82 um for the 
TMG with 40 um NI-Mn-Ga film designed for experiments. If the whole 
thickness could be active material that would result in large power output as 
the magnetization of the active material is now larger and, hence, a higher 
stiffness of the cantilever and a higher frequency of operation could be 
attained as the thermal mass would be comparable. Each TMG generates a 
significant amount of power density up to 30 mW/cm? with respect to its TM 
film when operated in a parallel configuration. A parallel architecture of 
numerous TMGs is necessary to raise the power level for particular purposes, 
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because the absolute power of a single TMG is not enough for direct applica- 
tion. However, placing numerous TMGs close to one another with a little gap 
between them raises the problem of cross-coupling, which could impair 
resonant self-actuation. When two devices with close tolerances are operated 
in parallel under the same magnet, the separation does not appear to have any 
impact at lower source temperatures when the frequency, stroke, and power 
output are already low. The TMGs don’t appear to have any coupling effects 
even at an extremely close spacing of 300 um. When the Tsourse increases, the 
TMG begins to exhibit coupling effects at close distances. The power appears 
to continuously decrease as the spacing (D) is decreased below 1 mm when 
the Tsource is 150 °C or above. The frequency of operation slightly peaks and 
causes an increase in the device's output voltage before it ceases operation at 
small separation. Which suggests a heat buildup in the TM film. 


Table 5 Performance results when the TMG is upscaled from having 10 um thick Ni-Mn-Ga 
film to 40 um thick Ni-Mn-Ga. Data pre-published in Joseph et al.”* 


Film Footprint Frequency Stroke AT Magnetic n Electrical 


Thickness power per Power 
footprint 

(um) (mm?) (Hz) (mm) (K) (mW/em’) (%) (uW) 

10 13.4 114 1.8 10.1 0.8 0.047 3.1 

20 12 136 1.9 54 1.1 0.046 4 

40 10.6 145 2.1 2.9 1.05 0.041 5.3 

40? 10.6 146 2.1 3.2 | 1.23 0.045 5.7 


Even though not all of the data exhibits a mutual correlation, it can be ob- 
served that it is possible for mutual correlation to exist when the devices 
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operate in parallel where the performance of one TMG is dependent on the 
neighbor. The ability of the TMG to cool itself more effectively as the neigh- 
boring device ceases functioning is the cause of the mutual correlation. On 
the other hand, this also demonstrates the TMG’s capacity to modify its 
performance in accordance with the surrounding ambient condition. When 
the correlation is missing, it can be inferred that a pause in the operation of 
the nearby device is not enough to produce an increase in the device’s ther- 
mal performance in some cases. According to the observations, the airflow 
between the devices has no further mechanical effects besides an indirect 
impact on the TMG thermal dissipation. This could potentially be evident as 
additional oscillations before the device ceases operation. However, since 
only the frequency spike and amplitude increase are discernible in the raw 
data, it can be said that only thermal coupling effects exist between the 
TMGs operating in parallel under a single magnet, without any mechanical or 
electromagnetic coupling effects. 


Table 6 lists the key parameters of the TMG demonstrator considered for this 
study. When three devices are working in parallel, the coupling effects are 
present for the center device starting at a Tsouree of 150 °C and deteriorate 
when the distance between the devices is reduced to less than 1 mm. In the 
case of two parallel devices, the average frequency is constant, and the only 
rise in frequency that is noticed locally is close to the intermittent pause in 
the TMG operation. This is where the average frequency increase differs 
from the case of two parallel devices. When three devices are operating in 
parallel, the same traits are observed, but there is also a rise in the frequency 
of operation. At higher temperatures, as the spacing decreases below 1 mm, 
the power output decreases. 


As the devices become closer with a spacing of under 1 mm at 140 °C, the 
power output increases slightly. This shows that somehow the TMG in the 
center has found a different stable oscillation frequency and stroke when 
compared to the TMG devices in the end. The thermal equilibrium that is 
reached when resonant self-actuation is achieved results in the high frequen- 
cy operation and large stroke. In this case the thermal equilibrium that is 
reached by the TMG at center to achieve resonant self-actuation is different 
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from the TMG at either ends of the array. This is also visible when the TMG 
are spaced at 5 mm from each other. Two device arrays showed no effect at 
5 mm distance, while the center device in three device array showed low 
output power and frequency even with 5 mm spacing. Hence, it shows that 
the TMG is operating at suboptimal conditions which is made worse by 
bringing the TMG devices closer. This causes an increase in output as the 
heat dissipation is affected causing the device to slowly buildup heat. This 
results in the increased frequency as the separation is reduced. 


Figure A 5 and Figure A 6 show the effect of resistance on the performance 
of the TMG when separation is varied. When operating the devices individu- 
ally, the highest power is still recorded at 400 Ohm. When the resistance is 
lower than 300 Ohm or lower, the power appears to grow as the separation 
decreases, with the highest power being recorded for 200 Ohm. Similar 
behavior is exhibited by the frequency, which has peaked at 200 Ohm and 
600 Ohm and a relatively low frequency at 400 Ohm. Frequency reduces as 
the separation is increased. This performance difference under different load 
resistances implies that the new thermal equilibrium at which the resonance 
is achieved is no longer optimum and has the potential to be improved further 
by device design. 


Table6 TMG parameters for the devices considered for experiments. Part of data is pre- 


published in Joseph et al.”*® 


Parallel operation with two devices 


Parameter Unit Value 
Ni-Mn-Ga width mm 2 
Ni-Mn-Ga length mm 2 
Ni-Mn-Ga thickness um 10 
Cu-Zn cantilever width mm 2 
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Cu-Zn cantilever length mm 
Cu-Zn cantilever thickness um 
Number of turns in pick-up coil - 


Parallel operation with three devices 


Parameter Unit 
Ni-Mn-Ga width mm 
Ni-Mn-Ga length mm 
Ni-Mn-Ga thickness um 

Cu-Zn cantilever width mm 
Cu-Zn cantilever length mm 
Cu-Zn cantilever thickness um 


5.4-6.4 
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Parts of this section has been published in Joseph et al.'°? Thermal energy 
below 200 °C is considered low-grade waste heat as this energy is usually 
taken as lost to the ambient. By adjusting the TMG design parameters men- 
tioned in section 5 and choosing materials with lower Te, it is possible to 
achieve TMGs operating in resonant self-actuation mode that are able to 
convert thermal energy from heat sources from below 100 °C. This enable an 
enormous variety of applications as the heat source temperature would be 
below the boiling temperature of water. As the temperature of the heat source 
is reduced further, TMGs operating with a heat source proximal to ambient 
temperature could be achieved. But when the heat source gets closer to the 
ambient conditions, the thermal performance of the TMG would be heavily 
affected. In this section, TMGs with TM film material of Gadolinium (Gd) 
and hydrogenated Lanthanum Iron Silicon alloy (La-Fe-Si-H) are investigat- 
ed. 


7.1 TMG with Gadolinium 


The TMG with Ni-Mn-Ga as the active material have shown that the concept 
of resonant self-actuation enables miniature TMG devices to be competitive 
with thermoelectric devices.”*°? However, the range of suitable heat source 
temperatures was restricted to 100-170 °C due to the ferromagnetic transition 
temperature Te of the Ni-Mn-Ga films being 98 °C. To utilize the notion of 
resonant self-actuation for the recovery of low-grade thermal energy, further 
suitable thermomagnetic film materials with Te around room temperature 
must be developed.'”® Gd films with a Te around room temperature, a large 
magnetization (M) under a moderate magnetic field, and a large temperature- 
dependent change in magnetization (AM /AT) are promising prospects. 


‘Section 7.1 has been published in similar form by the author in Joseph et al.!” 
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The ease with which they may be removed from the substrate on which they 
are placed helps their incorporation into devices.!5 In this study, a miniature- 
scale Gd-film TMG device is presented that operates in resonant self- 
actuation mode and associated physical performance qualities are analyzed 
near ambient temperature for lowering heat source temperatures from 75 °C 
to 40 °C. 


7.1.1 Material properties 


The active material in this study is a Gd film with a thickness of 40 um. 
Figure 90a depicts typical temperature-dependent magnetization properties 
for applied magnetic field levels ranging from 0.1 to 0.5 T. At the Curie 
temperature Tc, a substantial rapid change in magnetization AM occurs 
within a small temperature change AT induced by a ferromagnetic to para- 
magnetic phase transition. The transition is hysteresis-free due to its second 
order character. An Arrott plot analysis of isothermal M(T) readings yields a 
Te value of 20 °C, which is consistent with the stated value for bulk Gd. 
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Figure 90 (a)Magnetization measurement taken by varying temperature at different magnetic 


fields (b) Magnetization gradient, derived from the magnetization measure- 


ment. Data is pre-published in Joseph et al.!? 


Figure 90b depicts the temperature-dependent change in magnetization M/T. 
An extraction magnetometer was used to characterize the magnetic properties 
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of the Gd film. The maximal AM/AT point shifts somewhat field- 
dependently from 17 °C at 0.1 T to 19 °C at 0.5 T. The substantial variation 
in AM/AT near ambient temperature of up to 2.3 Am*kg!K"! is the primary 
impetus for developing a Gd film-based generator in this study. The maxi- 
mum value of AM/AT is spread into wider temperatures at higher magnetic 
fields. 


7.1.2 Design 


Upscaling of TMGs by array showed that operating TMGs close to each 
other is possible by proper design and parameter adjustment. One of the 
major hurdles to operating devices so close is the free-hanging copper wire 
from the pick-up coil bonded to the connection pin. Since this is through the 
side, it takes up space and sometimes gets stuck to nearby devices. 


permanent magnet 

a) Gd film 
bonding layer 
Cu-Zn cantilever 
pick-up coil 


substrate 


connector pin 


Figure 91 (a) 3D schematic of the double beam TMG with Gadolinium film at the tip, (b) side 
view of the TMG showing the component names. Adapted image from Joseph 
et al.!”? 


To avoid this, fixed electrodes are desirable. One solution to circumvent this 
issue is to split the cantilever beam into two independent beams electrically 
separated and connect the pick-up coil to the cantilever itself at the tip of the 
cantilever. Since the cantilever is made from brass with a lateral width of 
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0.9 mm and 25 um thickness, this approach results in electrical resistance of 
only 0.0165 Ohm. 


a) b) 


Figure 92 (a) Device with overhanging copper wires from pick-up coil 
attached to the connector (b) TMG with double beam brass cantilever as 
electrodes with no overhang. Adapted image from Joseph et al.” 


The resistance of the copper wire (15 um diameter) for the same length is 
larger with a value of 0.68 Ohm. This is due to the cross-sectional area being 
much smaller even though copper has only one third value of electrical 
resistivity compared to brass. Figure 91 shows the new design for the TMG 
with a double beam cantilever with Gadolinium as the TM active material 
and Figure 92 shows the visual comparison between the previous design and 
the new design. 


7.1.3 Fabrication 


The Gd film was created by triode sputtering a 3N purity Gd target onto a 
spinning Si substrate with a diameter of 100mm at a temperature of 
400 °C.'5* Ta buffer and capping layers of 100 nm thickness were also depos- 
ited at 400 °C to preserve the Gd layer. The Ta/Gd/Ta trilayer laminate may 
be easily peeled off the Si substrate because of the higher deposition tempera- 
ture. 154.155 The hcp crystallographic structure of the Gd layer was verified 
by X-ray diffraction. The film is then removed using a tweezer from the 
substrate by pulling. The cantilever is prepared by using an ultra-short, 
pulsed femtosecond laser with a double cantilever shape by keeping the 
connection at one end. Figure 93 shows the steps involved in the fabrication 
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ofthe TMG with Gd film with a double beam cantilever. Step 1 involves the 
preparation of an alumina substrate using a diamond saw. The cantilever is 
bonded to a ceramic substrate (Alumina) using a two-part non-conductive 
bonding adhesive in step II. The pick-up coil is fabricated using the semi- 
automatic coil winder presented in the section 2.7.8. The ends of the copper 
wires are attached to the bottom side of each of the cantilever beams using a 
conductive bond containing silver particles in step II. In step IV the pick-up 
coil is flipped and bonded to the bottom of the cantilever beam using a non- 
conductive bonding layer on top of the electrical bond made with the copper 
wires. This is done to ensure that the electrical bond is provided with addi- 
tional strength to withstand the high impact forces during operation. Sand- 
wiching the electrical connection between the pick-up coil and the cantilever 
ensures that the electrical connection is not affected and does not develop 
additional contact electrical resistance during operation. 


Figure 93 Fabrication steps involved in the preparation of a TMG with Gadolinium as the TM 
active material. Adapted image from Joseph et al.!? 


The Gd film of 40 um thickness is mechanically cut to 2 mm x 2 mm using a 
mechanical cutter and is fixed on the top side of the cantilever using a non- 
conductive bonding layer in step V. Step VI comprises severing the excess 
connected end of the beam after checking its dimensions. A connector bush 
pin is bonded to the substrate and then electrically bonded to the cantilever 
end using conductive adhesive in step VII completing the assembly. Electri- 
cal resistance is checked to ensure the internal resistance is in the range of 
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200 Ohm — 220 Ohm, to ensure no additional contact electrical resistance is 
developed during fabrication. 


7.1.4 Performance 


The cantilever tip shows time-resolved deflections with large strokes and 
high frequencies when resonant conditions are fulfilled, as seen in Figure 94a 
and b. It is important to highlight that resonant self-actuation is bound within 
a narrow range of source temperatures, Tsource. In the case shown in Figure 
94, at an ambient temperature of 11 °C, resonant self-actuation begins around 
40 °C and suddenly ends at 71°C. The stroke maintains above 1.5 mm 
throughout this temperature range and reaches a maximum of 2 mm at 60 °C. 
The frequency, on the other hand, begins at 65 Hz and progressively climbs 
to 106 Hz at 65 °C before abruptly decreasing at the critical maximum source 
temperature. Figure 93c shows how the stroke and frequency performance 
connects with the power output. The power output is already 1.3 uW at the 
lower limit of the heat source temperature of 40 °C, equivalent to a power per 
footprint of 10.3 uW/cm? (8.1 mW/cm? power per Gd volume), which is 
among the highest values recorded for Gd-based TMGs when the temperature 
difference between heat source and ambient (Tsource-Tamb) is only 29 °C. 
When the temperature of the source is raised, the power output increases 
gradually at first, then rapidly, peaking at 3 uW at 65 °C, equating to a power 
per footprint of 23.8 uW/cm? (18.8 mW/cm?). The load resistance of the 
electrical circuit was adjusted to its ideal value throughout these experiments. 


In these investigations the electrical circuit's load resistance was adjusted to 
its ideal value of roughly 500 Ohm, which is larger than the internal re- 
sistance of the pick-up coil, which is around 220 Ohm. The power output 
may be enhanced further by raising the load resistance above 220 Ohm, as 
shown in Figure A 7 in Appendix, since the reduction in electromagnetic 
damping leads to an increase in stroke and frequency. The capability of the 
Gd-based thin film TMG to operate under resonant self-actuation conditions 
is critical for practical applications, particularly when operating near room 
temperature. As a result, the ambient temperature range, Tamb, within which 
resonant self-actuation occurs for the specified Gd film must be investigated. 
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Figure 94 Mechanical and electrical performance of the Gd-film TMG as a function of heat 
source temperature (Tsource) at 11 °C ambient temperature: a) Cantilever-front 
stroke during resonant self-actuation, b) cantilever oscillation frequency, c) 
electrical power output at a load resistance of 500 Ohm. Data pre-published in 
Joseph et al.'°? 


Figure 95 depicts the actuation strokes, oscillation frequencies, and electrical 
power output of the TMG for various Tamb values ranging from 9 to 19 °C, 
while the source temperature, Tsource, varies between 35 and 75 °C. The 
results show that resonant self-actuation occurs within a narrow range of 
ambient temperatures, and that the TMG can maintain it up to a Tamb of 
19 °C, beyond which it becomes unsustainable. At 19°C, the TMG still 
generates roughly 1 uW at 50 °C and 1.4 uW at 65 °C Tsource. Furthermore, 
when Tamb increases to 19°C, the allowed range of source temperatures 
moves towards higher values. For example, the lower and higher limits of 
allowed source temperatures shift from 40 °C <Tsource <70 °C at Tamb of 9 to 
50 °C <Tsouree <75 °C when Tamb is increased by 10 °C. In terms of stroke and 
frequency, a maximum stroke of roughly 2 mm occurs at Tamb of 11 °C, while 
the stroke gradually decreases with rising Tamb, reaching a minimum value of 


155 


7 TMG near ambient temperature 


800 um at Tamb of 19 °C. Furthermore, changes in Tamb inside the resonance 
area have no effect on frequency. However, at low source temperatures, the 
frequencies are low and tend to decline more as Tamb increases until resonant 
self-actuation ceases to occur. 
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Figure 95 Effect of ambient temperature on the performance of the device, (a) stroke of cantile- 
ver deflection, (b) frequency of cantilever oscillation, (c) power output. Re- 
printed from Joseph et al.” 


Furthermore, the TMG's power output has a clear maximum at Tamb of 11 °C, 
which corresponds to the maximum stroke and frequency at source tempera- 
tures ranging from 55 to 70 °C. Figure A 8 to Figure A 12 provide further 
charts summarizing the mechanical and electrical performance of the TMG 
device as a function of load resistance and heat source temperature at various 
ambient temperatures. Heat intake and heat dissipation owing to conduction 
and convection, as well as temperature fluctuations in the Gd film, were 
simulated using LEM. Figure 96 depicts an overview of the heat intake and 
temperature change in Gd film. This demonstrates that the heat intake is 
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proportional to the heat source temperature Tsource. The largest heat absorption 
occurs at low ambient temperatures, ranging from 70 mW at 40 °C to 97 mW 
at 70 °C. Heat intake reduces with rising ambient temperature, with minimum 
values of 40 mW recorded at an ambient temperature of 19 °C. The heat 
transfer coefficient (hr) at contact between the Gd film and the heat source is 
the most important parameter determining heat intake. According to LEM 
simulations, hr should be as large as feasible. A heat transfer coefficient more 
than a critical value is necessary, at which steady resonant oscillation with a 
large stroke and frequency occurs. The needed minimal heat transfer coeffi- 
cient in the present the case of the TMG with Gd as active material is deter- 
mined to be 4000 Wm’K°', which was consistent with observations in section 
5.2.1 and published data.’8?!!” 


The results shown in Figure 97 show that the simulated heat dissipation from 
the Gd film to the cantilever by heat conduction through the bonding layer 
has the same Tamb and Tsource dependency as the heat intake. The largest heat 
dissipation by heat conduction was recorded at lower ambient temperatures, 
ranging from 60 mW at Tsource of 40 °C to 79 mW at Tsource of 70 °C. 
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Figure 96 LEM simulation of the effect of ambient temperature on thermal performance (a) 
average heat intake into the Gadolinium film on contact, (b) average tempera- 
ture change of the Gadolinium film. Reprinted from Joseph et al. "° 
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It is worth mentioning that heat conduction accounts for the major proportion 
of heat dissipation when compared to heat convection. As a result, the canti- 
lever is functioning as a critical primary heat sink. Heat convection from the 
Gd film to the ambient air, on the other hand, correlates with the stroke and 
frequency of oscillation. The consistence between heat dissipation by con- 
duction and convection and simulated heat intake validate the overall con- 
sistency of the data. The bonding layer's thermal resistance (R») between the 
Gd film and the cantilever is identified as the essential parameter influencing 
conductive heat dissipation. As a result, its value must be tailored during the 
TMG device's design to match the heat intake during contact with the heat 
source for effective heat transfer. 
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Figure 97 LEM simulation of the effect of ambient temperature on the loss mechanism (a) 
conduction losses from the Gadolinium film into the cantilever, (b) convection 


losses into the ambient air from the Gadolinium film. Reprinted from Joseph et 
al.'52 


Temperature changes (AT) in Figure 98b are relatively small due to the high 
frequency of oscillation and short heat transfer durations, ranging from 
1.2 °C at low ambient temperatures of 9 °C to around 0.45 °C at high ambient 
temperatures of 19 °C. Despite minor AT values, the resonant self-actuation 
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conditions are still sustained and the TMG can operate in a steady manner. 
Aside from the AT seen in Figure 98b, the cantilever temperature Tcant is 
another important thermal parameter to examine. Figure 98a shows the 
simulated cantilever temperature Tcant at various ambient and source tempera- 
tures. Tcant rises with rising Tamb and source temperatures Tsource, with mini- 
mum values of 25 °C and highest values of 34 °C. However, increasing Tcant 
causes a decrease in heat dissipation from the Gd layer, limiting heat input. 
As illustrated in Figure 96, a decrease in heat dissipation and an increase in 
Tamb are both connected. As a result, the rise in Tcant, heat dissipation, and 
heat intake for rising ambient temperature are all interconnected. These 
results are interesting because a small temperature change (AT), as shown in 
Figure 96b, is sufficient for stable operation in resonant self-actuation mode. 
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Tcant is observed to be 28.4 °C at the optimal power point, which occurs at an 
ambient temperature of 11 °C and a source temperature of 65 °C. The tem- 
perature of the Gd film (Trim) oscillates roughly between 28.5 and 29.4 °C 
throughout this time, resulting in a AT of 0.9 °C, as shown in Figure 98b. 
Figure 98c depicts the equivalent heat intake in small pulses lasting roughly 
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0.5 ms, each of which coincides with the contact between the heated magnet 
and the Gd layer. Following cooling takes roughly 9 ms. In contrast, with an 
ambient temperature of 19 °C and a source temperature of 50 °C, the lowest 
limit of operation is found. The cantilever temperature Tcant is roughly the 
same in this situation, at 28.6 °C, while the temperature of the Gd film Tfilm 
only shows a AT of only 0.45 °C. According to our observations, this value 
of AT acts as a lower limit for resonant self-actuation. Below this threshold, 
the magnetization change AM falls below the critical limit, and the magnetic 
attraction force can no longer compensate for damping losses. 


7.2 TMG with La-Fe-Si-H 


To achieve high power output and to develop TMGs that are operating at 
ambient temperatures above and near room temperature, TM films with Te 
near room temperature are required. La-Fe-Si (Lanthanum Iron Silicon) alloy 
is a very promising alternative with tunable T. by hydrogenation, and Man- 
ganese or Cobalt addition around the region of 27 °C to 50 °C. The advantage 
of La-Fe-Si alloy is the very large magnetization values and sharp drop at Tc. 
These properties are very important for the presented TMGs to operate near 
room temperatures with high power output and frequency. 


7.2.1 Material properties 


The magnetization of the La-Fe-Si-H film reaches a maximum value of 
171 Am?/kg when the material is ferromagnetic. After the transition at 48 °C, 
the material still retains a magnetization value of 15 Am’/kg. The high value 
of 14.6 Am?/kgK for AM/AT is unprecedented and results in a very sharp 
drop in magnetization within a few degrees. It can be seen in Figure 99 that 
the whole transformation is confined to a few degrees around Tc. Figure 99 
shows that the Curie temperature (Tc) is slightly shifting to higher tempera- 
tures at higher magnetic fields. The average value of Te is 48.3 °C. The La- 
Fe-Si-H film is prepared by melt spinning. By arc melting in an environment 
of Argon (Ar) gas, LaıFe11.6S11.4 alloys are prepared. 
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Figure 99 (a)Magnetization measurement for the La-Fe-Si-H sample, (b) AM/AT for La-Fe-Si-H 
sample with respect to changing temperature. 


Then ingot fragments are placed through a nozzle into a quartz tube. High- 
purity Ar is then injected into the chamber after the air had been evacuated to 
create a vacuum. By employing a pressure difference, the samples are induc- 
tion melted and expelled through the nozzle. Cu wheel is spun, and the 
molten alloy solidifies into a ribbon as it gets in contact with the Cu wheel as 
mentioned in section 2.7.2. The as-spun ribbons are subsequently annealed. 
Further information regarding the preparation of La-Fe-Si samples using melt 
spinning can be obtained from Gutfleisch et al.!°° Unlike DC RF magnetron 
sputtering, the melt spinning technique does not yield a very smooth surface. 
The surface is rough and further polishing is limited due to the material being 
very brittle and would break with the slightest pressure on the films. But the 
advantage is that it is possible to achieve a large thickness in the range of 
50 um which is large compared to DC RF magnetron sputtering of Ni-Mn-Ga 
films which is limited to 10 um due to film homogeneity. The samples are 
then hydrogenated to bring the T.e above room temperature to finish the 
fabrication of La-Fe- Si-H ribbons. 
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7.2.2 Fabrication 


The device fabrication and design of the TMG are similar to the TMG with 
Ni-Mn-Ga films. The fabrication process involved is described in section 
4.4.2. 50 um La-Fe-Si samples are cut using an ultra short pulse femtosecond 
laser to a lateral size of 2 mm x 2 mm. A single cantilever beam is made from 
a brass sheet with a thickness of 30 um, which is also structured using the 
laser. The pick-up coil is fabricated out of 15 um enamel-coated copper wire 
with a total of 400 turns around a core of 1 mm? area. The cantilever is 
attached to the substrate using a thermally conductive but electrically isolat- 
ing bonding layer. Both the pick-up coil and La-Fe-Si are attached to the tip 
of the cantilever, with the former on the bottom side and later on the top side 
of the cantilever using a non-conductive bonding layer. The ends of the 
copper wires from the pick-up coil are attached to the connector bush pin 
which is connected to the substrate. The electrical connection is made using 
an electrically conductive bonding layer. The surface of the La-Fe-Si-H 
sample is polished by hand to improve the surface roughness. This step is 
done last as the material is very brittle, causing it to break if attempted to 
polish in a free-standing state. 


7.2.3 Performance 


Since the AM/AT value for La-Fe-Si-H film is very large and confined to a 
few degrees, the resulting device would exhibit very high performance in a 
small temperature window and would experience performance loss as the 
conditions are changed. This reduction in tolerance is a consequence of the 
very small temperature range in which the material transforms, unlike Ni- 
Mn-Ga films or Gd films both of which exhibit a more gradual change in 
magnetization over a wide range of temperatures. Figure 100 shows the 
performance measurement of the TMG with La-Fe-Si-H as TM film. The 
large power output is confined to a single Rt of 400 Ohm. At 90 °C source 
temperature (Tsource) the device is capable of generating 4.2 uW. The TMG 
can sustain the resonant self-actuation from a minimum Tsource of 70 °C which 
is limited by the T. of the material. At a Tsource of 70 °C, the maximum power 
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output of the TMG is around 2.2 uW, which increases as the Tsource is in- 
creased. But after increasing the Tsource beyond 90 °C, the performance of the 
device reduces. The device loses resonant self-actuation when the TMG is 
operated with a Tsource of more than 
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Figure 100 performance of the TMG with La-Fe-Si-H as TM film, (a) stroke of the deflection of 
the cantilever tip during resonant self-actuation, (b) frequency of oscillation of 
the cantilever, (c) power output from the pick-up coil 


100 °C. LEM simulations are used to determine the heat intake in the La-Fe- 
Si-H film on contact. Figure 101a shows the heat intake values correspond- 
ing to the value of the Tsource. It can be seen that the average value of heat 
intake is much lower than that of the TMG with Gd as TM film and much 
lower than that of the TMG with Ni-Mn-Ga as Tm film. Figure 101b illus- 
trates that the TMG is capable of sustaining resonant self-actuation even 
though the AT within the material is below 1 °C. The AT is 0.22 °C at a Tsource 
of 70 °C and increases to 0.5 °C at a Tsource of 100 °C. 
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Figure 101 LEM simulation of the thermal processes, (a) heat intake into the La-Fe-Si film, (b) 
temperature change within the La-Fe-Si film 
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Figure 102 LEM simulation of thermal losses, (a) conduction losses from the La-Fe-Si film to 


the cantilever through the R», (b) convection losses from the film directly to the 
ambient. 
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The TMG is able to operate in resonant self-actuation with such low AT in 
the material due to the very large AM/AT. Figure 102 shows the LEM 
simulation values of heat dissipation from the La-Fe-Si-H film. This shows 
that conduction is the main form of dissipation and the trend of the change in 
conduction losses is similar to the change in heat intake. The change in 
convection loss on the other hand is slightly different from the conduction as 
this value is also dependent on the Tamb, frequency of operation, and the 
convective heat transfer coefficient. The TMG is operating at a maximum 
frequency of 137 Hz even with the low contact heat transfer coefficient (hr), 
suggesting that the performance of the TMG could be increased many folds 
and the minimum Tsource could also be reduced further by increasing the hr 
between the magnet and the TM film, which could be done by minimizing 
the surface roughness of the TM film. 


7.3 Discussion 


Because of its capacity to function in the resonant self-actuation mode, which 
allows for high frequency and stroke, film-based thermomagnetic generators 
(TMGs) are capable of efficient conversion of thermal energy into electrical 
energy. However, the useful range of heat source temperatures for these 
devices is restricted by the thermomagnetic material's ferromagnetic to 
paramagnetic transition temperature (Tc). For example, a Ni-Mn-Ga film 
with a Te of 98 °C at 0.5 T can only function between 100 and 170 °C. To 
recover low-grade waste heat, appropriate thermomagnetic films with Te 
close to room temperature must be designed and applied. Gd, which has a Te 
of 20 °C, is one interesting candidate material for this purpose. 


The thermo-magneto-mechanical performance of a Gd-film TMG at room 
temperature is investigated in this work as the heat source temperature is 
gradually reduced from 75 °C to 35 °C to assess its potential for low-grade 
waste heat recovery. The LEM simulation model of the TMG is reported in 
section 4.2 and presented also by Joseph et al.!”® to describe the TMG 
performance. This entails regulating the resonance frequency, improving the 
heat input, matching heat dissipation, and determining the best performance 
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parameters for optimum power and efficiency during resonant self-actuation 
at the highest feasible ambient temperature. While film-based TMG devices 
exhibit self-adapting performance, research on a Gd-film TMG discloses a 
severely limited range of ambient and source temperatures necessary for 
resonant self-actuation. 


Table 7 Modeling parameters for the LEM simulation ofthe TMG with Gadolinium. Data pre- 
published in Joseph et al." 


Mechanical parameter Unit Value References 
Pick-up coil mass Kg 6x 10° | This work 
Cantilever density Kg/m? 8500 ue 
Gadolinium density Kg/m? 7900 14 
Bonding layer density Kg/m? 920 1 
Structural and air damping Ns/m 2.8x 107" This work 
Impact damping Ns/m 0.1 This work 
Contact stiffness N/m 10000 This work 
Young’s modulus Pa 1x 1011 146, 
Cantilever length m 0.0045 This work 
Cantilever width m 0.0009 (each) This work 
Cantilever slit (gap between m 0.0002 This work 
beams) 
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Cantilever thickness 


Electrical & Magnetic Parame- 
ters 


Remanent magnetic field 


Number of turns 


Area of pick-up coil 


Pick-up coil internal resistance 


Load resistance (RL) 


Thermal Parameters 


Cantilever thermal conductivity 


Gadolinium thermal conductivity 


Bonding layer thermal conduc- 
tivity 


Cantilever heat capacity 


Gadolinium heat capacity 


Bonding layer heat capacity 


Unit 


0.000025 | This work 


Value References 


1.07 This work 


400 This work 


2.81 x 1076 This work 


250 This work 


500 This work 


Value References 
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Contact heat transfer coefficient W/m?K 4600-6600 | This work 


Convective heat transfer coeffi- W/m?K 40-140 This work 
cient 


The highest limit of ambient temperature (Tamb) corresponds to the tempera- 
ture at which AM/AT is maximum (19 °C at 0.5 T), beyond which resonant 
self-actuation is no longer sustainable, regardless of heat source temperature 
(Tsource). TMG devices, at Tamb = 19 °C, can still generate 1 yW at minimum 
and maximum heat source temperatures of 50 and 75 °C, respectively. At 
lower ambient temperatures, such as Tamb = 11 °C, resonant self-actuation 
quickly kicks in, yielding a power output of 1.3 uW at lower heat source 
temperatures, with Tsouree = 40°C. At just 29°C temperature difference 
between heat source and ambient, the associated relative power per footprint 
of 10.3 u W/cm? (8.1 mW/cm? power per Gd volume) is among the highest 
values found for Gd-based TMGs. The resonant self-actuation of the TMG is 
lost abruptly above 71 °C. The efficiency of the TMG, given by the ratio of 
electrical output power and thermal power intake, is 0.0021 % at Tsource = 
40 °C, while the corresponding relative efficiency, given by the ratio of 
absolute and Carnot efficiency, is estimated to be 0.045 %. Thereby, the 
cantilever temperature is assumed to be the heat sink temperature, as it 
accounts for 80 % of heat dissipation. The efficiency of the related Gd film 
material is 0.02 %, as evaluated by the ratio of output magnetic energy to 
thermal input energy, where output magnetic energy is obtained from the 
thermomagnetic cycle of the Gd film illustrated in Figure A 13 (Appendix). 
The notion of resonant self-actuation does not rely on utilizing the whole 
temperature range of the magnetic transition for energy conversion, but rather 
on generating the highest electrical power output with the least amount of 
temperature change of the film (Trim), rather than optimum. 


The La-Fe-Si-H boasts a much higher value for the magnetization in its 
ferromagnetic state and also shows a very sharp drop in magnetization with a 
very high value of 14.6 Am?/kgK for the AM/AT. This sharp drop comes with 
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advantages like high frequency of operation due to the very low AT required 
in the material for achieving the required magnetization to release the canti- 
lever tip from the magnet surface. The fabrication process of melt spinning 
does present a set of challenges for achieving resonant self-actuation with the 
prepared La-Fe-Si-H samples. The surface roughness results in a very low 
thermal contact between the magnet and the La-Fe-Si film. The minimum 
Tsource is not only limited by the Te of the material, but also by the very low hr 
between the TM film and the magnet. Even though the magnet is polished to 
a mirror finish, the TM film had a very rough surface due to the fabrication 
process (melt spinning) The roughness is large enough to cause a major 
reduction in hr. This is confirmed by LEM simulations by matching the 
performance of the TMG suggesting a hr value between 2800 and 3200 
W/m?K.Thus, even though the device does possess the potential to operate at 
higher frequencies (over 200 Hz, estimated from the natural frequency of the 
assembly), the device operation is limited to a maximum frequency of 
137 Hz due to the heat transfer times required to intake required amount of 
heat at contact. The material also has a very low thermal conductivity of 
8 W/mK which further restricts the heat dissipation from the material. 


Table 8 Modeling parameters for the TMG with La-Fe-Si-H as the active TM film 


Mechanical parameter Unit Value | References 

Pick-up coil mass Kg 6x 1076, This work 

Cantilever density Kg/m? 8500 13 

La-Fe-Si-H density Kg/m? 7245 152 

Bonding layer density Kg/m? 920 14 

Structural and air damping Ns/m 2.8 x This work 
x 10° 
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Impact damping 
Contact stiffness 
Young’s modulus 


Cantilever length 


Cantilever width 


Cantilever slit (gap between 
beams) 


Cantilever thickness 


Electrical & Magnetic 
Parameters 


Remanent magnetic field 
Number of turns 


Area of pick-up coil 


Pick-up coil internal resistance 
Load resistance (RL) 
Thermal Parameters 


Cantilever thermal conductivi- 


ty 


170 


Ns/m 


Unit 


Q 
Q 
Unit 


W/mK 


0.1 
10000 
1x 1011 


0.0053 


0.0009 
(each) 


0.0002 


0.000025 


Value 


1.07 
400 


2.81 
x 1076 


250 
500 
Value 


109 


This work 


This work 


146 


This work 


This work 


This work 


This work 


References 


This work 


This work 


This work 


This work 


This work 


References 
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La-Fe-Si thermal conductivity W/mK 8 9 
Bonding layer thermal con- W/mK 0.225 = 
ductivity 

Cantilever heat capacity J/kgK 400 = 
La-Fe-Si heat capacity J/kgK 600 a 
Bonding layer heat capacity J/kgK 2100 $% 
Contact heat transfer coeffi- W/m?K 2800-3200 This work 
cient 

Convective heat transfer W/m?K 40-140 | This work 
coefficient 


After several design optimizations including the Re, cantilever stiffness, and 
polishing of the surface of the La-Fe-Si-H film, resonant self-actuation is 
achieved. In resonant self-actuation, the TMG is able to generate a maximum 
power output of 4.2 uW at a Tsource of 90 °C. Since the temperature span in 
which the material transformation occurs is curbed to a few degrees, the 
performance of the device also shows similar trends. The device exhibits a 
lack of tolerance to the ambient conditions, unlike the TMG with Gd film. 
Even while changing the Rı, it has a major effect on the electromagnetic 
damping at the cantilever tip causing the performance to have a clear peak at 
400 Ohm. When the resistance is changed from this value, the performance 
of the TMG reduces drastically. This is also compounded by the lower field 
gradient experienced by the pick-up coil. 


In most TMGs the surface of the TMG active material is placed on the 
surface of the magnet surface in such a way that the pick-up coil experiences 
the higher field gradient at the edges of the magnets. But in this case, due to 
the lower contact hr, the TMG has to be placed close to the center of the 
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magnet causing a large reduction in the magnetic field gradient experienced 
by the pick-up coil during actuation even though the effect on the actuation 
itself is not much. This effect can be observed in Figure 100, where the 
frequency and stroke are reduced at 400 Ohm even though it is the point at 
which the power output is also the maximum. This is caused due to the 
additional current through the pick-up coil at this point. At other Rıs the 
conditions are different, and the device fails to achieve the same level of 
magnetic field gradient throughout the actuation cycle. In other TMGs, the 
tolerance of the device to adjust itself to the changing ambient conditions 
help the device to avoid such a performance peak and to form a more gradual 
increase or decrease in performance. 


It is also unprecedented that resonant self-actuation is achieved with just a AT 
of 0.22 °C. This is observed at a Tsource of 70 °C. Since currently, the surface 
roughness is still large to have a major impact on the hs, it can be assumed 
that the Tsource values including maximum, minimum, and the Tsource at which 
maximum output occurs would shift towards lower temperatures as the hr is 
improved. This can be achieved either by polishing the surface further or to 
coat the surface with a material with high thermal conductivity to even out 
the surface imperfections. The heat intake into the La-Fe-Si-H film at the 
cantilever tip is low when compared to the TMG with Ni-Mn-Ga film. 


Adding manganese to La-Fe-Si-H has shown great potential in reducing the 
Te of the material which is also another way in which the minimum operating 
temperature could be reduced further and brought closer to the room tempera- 
ture.' Figure 18 shows an example of a La-Fe-Si-Mn-H sample with its 
magnetization dropping abruptly at around 25 °C. This shows further poten- 
tial to bring the TMG operation temperature closer to room temperature and 
crossing the point of 37 °C for the Tsource. This would result in the develop- 
ment of miniature TMGs operating with body heat, which is of special 
interest for miniature scale medical applications. 


172 


7.3 Discussion 


Table 9 TMG parameters for devices considered for experiments. Part of data pre-published in 


Joseph et al.'°? 


TMG based on Gadolinium 


Parameter 

Gadolinium width 
Gadolinium length 
Gadolinium thickness 

Cu-Zn cantilever width (each) 
Cu-Zn cantilever length 
Cu-Zn cantilever thickness 
Number of turns in pick-up coil 
TMG based on La-Fe-Si 
Parameter 

La-Fe-Si-H width 

La-Fe-Si-H length 
La-Fe-Si-H thickness 

Cu-Zn cantilever width (each) 
Cu-Zn cantilever length 
Cu-Zn cantilever thickness 


Number of turns in pick-up coil 


Unit 


um 


Unit 


um 


Value 


400 


400 
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This section presents a detailed evaluation of the TMG performances pre- 
sented in Sections 4-7 and the influences of all design parameters. The TMG 
designs will be compared with respect to the target objectives, i.e. maximum 
power output, optimum efficiency and adjustment of the range of source 
temperature. 


8.1 Discussion of results 


8.1.1 Lumped element model 


The presented multi-physical lumped element model (LEM) of the cantilever- 
based TMGs has been validated with experimental results. Thus, it is used as 
a powerful tool to design and optimize the TMG performances. Identifying 
the different physical domains of the TMG is the first step in preparing a 
LEM. The presented TMG is identified to have four physical domains in- 
volved, namely (1) magnetic, (2) mechanical, (3) thermal and electrical 
domains. Each domain is modeled by taking into consideration its effect on 
other domains. Most importantly the LEM simulation is the key in ensuring 
that the design of the TMG is compatible with resonant self-actuation before 
the device is fabricated. 


(1) Magnetic domain: Modeling of magnetic properties of the perma- 
nent magnet used in the TMG is presented. The analytical model is 
capable of calculating generated magnetic field and magnetic field 
gradient of the permanent magnet, based on the dimensions of the 
magnet and the remanent magnetic field (dependent on the material 
composition of the permanent magnet). The magnetization of the 
material is temperature dependent, and the values are directly taken 
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(2) 


(3) 


(4) 


as a lookup table prepared from the experimental measurement of 
the magnetization. 


Mechanical domain: The mechanical oscillations of the cantilever 
beam with tip mass are modeled by solving the partial differential 
equations (PDEs) for oscillations in z direction corresponding to the 
1-Dimensional assumption of the spring damper system. The mag- 
netization of the active TM film, inertia, stiffness of the cantilever, 
and magnetic field gradient determines the effective force at the can- 
tilever tip, which would either cause the tip to bend towards the 
magnet or away from it depending on its direction. This effective 
force is modeled by the LEM by solving the thermal, magnetic, elec- 
trical, and mechanical sections simultaneously. Since each section of 
the model is dependent on the outputs of other sections of the model, 
the initial point is taken as the position of the cantilever at zero actu- 
ation with no reset force present. 


Thermal domain: The thermal section of the model shows the most 
sensitivity toward performance optimization and it is observed that a 
minute change in parameters in the thermal section would result in a 
large deviation in performance. The values like heat transfer coeffi- 
cient (hr) and convective heat transfer coefficients are assumed ini- 
tially to design the basic TMG. Once the design is completed, the 
prototype is fabricated and experimentally tested. The assumed val- 
ues are readjusted to match the experimental values to represent the 
TMG properly. 


Electrical domain: The electrical section of the model uses Faraday 
law of induction to determine the electrical power conversion in the 
pick-up coil. The mechanical oscillations of the cantilever tip corre- 
spondingly results in a varying magnetic field gradient in the pick- 
up coil generating a voltage. The generated voltage is converted to 
current accordingly based on the internal resistance of the TMG and 
the load resistance connected. The electromagnetic damping force 
resulting from the conversion of oscillations to voltage is also mod- 
elled in this section. 


8.1 Discussion of results 


LEM simulation allows to interpret the results of the experiments in a mean- 
ingful way and to understand the thermal and magnetic processes within the 
TMG, both of which are difficult to measure. The AT of the active TM film 
and the temperature of the cantilever are crucial in achieving resonant self- 
actuation and LEM simulation gives a detailed conception of these processes. 
In the presented example of the capability of LEM, a TMG based on 10 um 
Ni-Mn-Ga film as active TM film is analyzed in detail. The experimental 
analysis of the device is also completed and presented to validate the model. 


8.1.2 Impact of parameters on the performance 


Many parameters influence the design and performance of the TMG. These 
parameters could be broadly classified into (1) TM film parameters, (2) 
cantilever material parameters, and (3) geometrical and operation parameters. 
The TM film parameters depend on the composition of the material and heat 
treatment conditions. 


(1) TM film parameters 


Following are the material parameters that affect the design and performance 
of the TMG: 


e Magnetization (M) 

e Change in magnetization with temperature change (AM/AT) 
e Curie temperature (Te) 

e Surface roughness 

e Thermal conductivity 

e Thermal capacity 


Magnetization: The value of magnetization determines the maximum possi- 
ble stiffness of the cantilever as in the cold start case when the cantilever is at 
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its initial position, the magnetic force must be strong enough to bend the 
cantilever tip to the magnet surface. A high value of the magnetization would 
enable the design of the TMG with high frequency provided the thermal 
cycle is also capable of high frequency operation. 


AM/AT: The value of AM /AT affects the range of Tsource values in which the 

device is able to sustain resonant self-actuation and also the AT required for 
achieving self-actuation. A high value of AM/AT would result in very high 
frequency operation as the heat intake and temperature required for resonant 
self-actuation is low. A second-order transformation without hysteresis, high 
value of magnetization, and high value of AM/AT would be the ideal charac- 
teristics of the material for the TMG as this combination would result in a 
very high operation frequency and high power output. 


Curie temperature: Te in this work is taken as the temperature at which the 
material exhibits the maximum AM/AT. This value determines the operation 
temperature of the TMG in which the device would be able to sustain reso- 
nant self-actuation. TM film material with higher Te would require higher 
Tsource to achieve self-actuation and consequently even higher temperature to 
achieve resonant self-actuation. Material with closer to room temperature 
(slightly above) values of Te would enable the development of the TMG with 
near room temperature operation. 


Surface roughness: Section 2.5.6 shows how the surface roughness of the 
contact surfaces affects value of hr. Surfaces with higher surface roughness 
would result in a lower value of hr. This would mean a longer heat transfer 
time to achieve the same amount of heat flow. For the presented design of the 
TMGs, it is beneficial to have a smooth surface with a very low surface 
roughness to achieve high values of hr. 


Thermal conductivity and thermal capacity: Thermal conductivity and ther- 
mal capacity of the TM film would affect the heat intake and the heat dissipa- 
tion capability of the material. Material with higher thermal conductivity 
tends to have a lower value of thermal capacity and vice versa. The higher 
thermal conductivity and lower thermal capacity would result in a large AT 
for a small amount of heat flow. The heat flow required to achieve the same 
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temperature would increase with an increase in the thermal capacity of the 
material. 


(2) Cantilever material parameters 


The following material parameters are the ones affecting the design and 
performance of the TMG: 


e Young’s modulus 

e Thermal conductivity 
e Thermal capacity 

e Corrosion resistance 


Young’s modulus: The young’s modulus affects the stiffness of the material. 
The higher the Young’s modulus, the higher would be the stiffness for the 
same dimensions of the cantilever beam. 


Thermal conductivity & thermal capacity: The thermal conductivity of the 
cantilever material affects the thermal dissipation of the film. The convective 
heat transfer from the cantilever to the ambient air and the thermal conduc- 
tion from the film to the cantilever is both affected by the thermal conductivi- 
ty of the cantilever material. The higher thermal capacity of the material is 
also desirable, and this would result in a stable temperature in the cantilever 
during operation without much thermal cycling. 


Corrosion resistance: The material chosen for the cantilever should have high 
corrosion resistance. This is important for the durability of the TMGs and 
also to maintain the performance over time. 


(3) Geometrical and operation parameters 


Following are the important design parameters that have major effects on the 
performance of the TMG and its ability to achieve resonant self-actuation: 
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e Magnetic field 


e Magnetic field gradient 


Ambient temperature 
e Heat source temperature 


Stiffness of the cantilever 


Contact heat transfer coefficient 


Bonding layer thermal resistance 
e Pick-up coil turns 


Load resistance 


Magnetic field: The magnetic field is affected by the dimension of the mag- 
net and the material of the permanent magnet. The magnetic field at the 
surface increase as the aspect ratio (ratio between the height along the magnet 
along easy axis and lateral cross-sectional area) increases. But this would also 
result in a lower magnetic field at a distance from the magnet surface. A 
higher magnetic field is desirable both at the surface (during operation) and at 
a distance (for the cold start of the device) for a TMG with cantilever based 
design. 


Magnetic field gradient: Magnetic field gradient at the surface increases with 
the increased aspect ratio of the magnet dimensions. The magnetic field 
gradient is slightly lower at a distance for a higher aspect ratio. A higher 
magnetic field gradient is desirable to have a higher magnetic force and 
higher power output from the pick-up coil. 


Ambient temperature (Tamb): Ambient air is the final sink into which all the 
heat is dissipated. For a TMG, the heat sink of the device is the cantilever. 
Which finally dissipates it to the ambient air. Thus, the temperature of the 
ambient air is an important factor that determines if the device is capable of 
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sustaining resonant self-actuation. Which requires high rates of heat intake 
and heat dissipation to achieve the minimum AT in the material. 


Heat source temperature (Tsource): This forms a major part of the TMG design, 
as the TMG design could be adjusted to operate at different Tsource values. But 
the minimum value of Tsource is limited by the T. of the TM film. The TMG 
does tolerate a range of heat source temperatures at which it will exhibit 
resonant self-actuation. The heat source temperature affects the performance 
of the TMG and usually has an optimum Tsource at which the power output is 
maximum which is unique for the TMG depending on its other parameters. 


Stiffness of the cantilever: Maximum stiffness of the cantilever is limited by 
the magnetic force at the initial state of the cantilever position. The stiffness 
is given by both the material property Young’s modulus and the geometric 
properties like length, width and thickness of the cantilever. But the designed 
stiffness would be tailored to match the natural frequency of the beam with 
the thermal cycling frequency of operation. Thus, the final frequency is also 
affected by the thermal processes within the TMG. 


Contact heat transfer coefficient (hr): The heat flow at contact is the driving 
force of the TMG operation. This heat flow is governed by the hr at contact, 
which needs to be maximized to achieve optimum performance. Achieving a 
smooth surface at the contact is crucial to have large values for the hr A 
critical minimal value of hr of 8 kW/m?K is observed for a TMG with 10 um 
Ni-Mn-Ga film below which resonant actuation fails. 


Bonding layer thermal resistance (Rt): The value R» between the TM film 
and the cantilever beam is important for all TMGs. A critical minimal value 
of Rb of 5 K/W is observed for a TMG with 10 um Ni-Mn-Ga film. Below 
this value, the cantilever and the TM film becomes thermally coupled result- 
ing in failure to achieve resonant self-actuation. Increasing the thermal 
resistance results in a very sharp increase in the power output. But power 
output saturates when Rb is increased over 10 K/W. Showing only marginal 
increase when the Rois increased further. 
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Pick-up coil turns: The number of turns in the pick-up coil has a direct impact 
on the amount of energy converted from mechanical oscillation to electrical 
energy. The electromagnetic damping is caused by energy conversion from 
mechanical energy into electrical energy. Thus, it is desirable to have maxi- 
mum electromagnetic damping to maximize energy conversion. But too 
much electromagnetic damping would result in the device falling out of 
resonant self-actuation. Hence the number of turns of the pick-up coil needs 
to be optimized depending on the configuration of the TMG. A TMG with a 
material having a large magnetization value would possess large stiffness for 
the cantilever which could tolerate a higher number of pick-up coil turns as 
the device possesses higher magnetic force and reset force. 


Load resistance (RL): The Rı connected to the TMG would limit the current 
through the pick-up coil. Thus, higher Rı at the ends of the pick-up coil 
would result in lower electromagnetic damping due to the lower current 
flowing through the pick-up coil. Increasing the frequency and thus resulting 
in higher power output. The transfer efficiency is maximum when the Rı is 
matched to the internal resistance of the pick-up coil. Thus, the transfer 
efficiency is reduced as the Rı value is moving further away from the value 
of the internal resistance. But as the Rı is increased further, the transfer 
efficiency keeps on decreasing and after a point, it would result in a reduction 
in power output even when the frequency and stroke are increasing. This 
causes the system to have an optimum Rı at which the TMG would generate 
the maximum output. This value also depends on the design and configura- 
tion of the TMG. 


The validated LEM model is utilized to perform the parameter study to 
understand the effect of different design parameters on the performance of 
the TMG including the thermal processes within the TMG. It is observed that 
the critical minimal value of AT of the active TM layer determines the limit 
below which the TMG would lose its resonant self-actuation. This value is 
affected by both heat intake and heat. The heat intake is controlled by the 
temperature difference between the magnet and the TM film and the value of 
hr. But the hr is dependent on the contact area which is governed by the 
surface roughness of the contact surfaces. The magnet in each case of the 
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TMG is polished to a mirror finish to achieve maximum possible contact. 
The sputtered films have low surface roughness. Thus, the combination of 
these films and the polished magnet surface would result in a large portion of 
the film being in contact and hence a large value of hr. The melt spun film on 
the other hand would have a rough surface, these films need to be polished to 
reduce the surface roughness so that at least reasonable values of hr can be 
achieved between the magnet surface and the TM film on contact. 


8.1.3 Upscaling of TMG 


The method to increase power focuses on the tuning of the magnetic attrac- 
tion force, mechanical reset force, and heat transfer dynamics in order to 
sustain resonant self-actuation. As long as resonance conditions are main- 
tained, the mechanical energy of oscillatory motion increases with increasing 
mass and frequency, sustaining magnetic cycles with rising variations in 
magnetic energy. This increase in mechanical energy is what leads to the 
realized increase in power. The results demonstrate that the heat transfer 
kinetics for a specific ferromagnetic material and magnetic field gradient 
eventually limit the upscaling technique. 


According to LEM simulations, conductive heat losses serve as the primary 
source of heat dissipation, and heat input grows as film thickness increases. 
As a result, the corresponding fraction of average retained input heating 
power increases from about 80% to 94% with increasing film thickness 
from 10 um to 40 um, respectively. These fractions display the same thick- 
ness dependence at lower Tsource values and only slightly deviate from one 
another. This study shows that a sizeable fraction of heat is maintained in the 
SMA film stack and contributes to the AT. The AT of the film stack during 
oscillation, however, diminishes as the film thickness rises, which has two 
significant effects. First, as AT reduces, the effectiveness of increasing film 
thickness decreases. Second, the continued scaling up of film thickness is 
severely limited by the declining AT due to the increased thermal mass of the 
stack. 
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The absolute power of a single TMG is insufficient for most applications. 
Therefore, a parallel architecture of many TMGs is required to boost the 
power level for specific uses. Cross-coupling, which could hinder resonant 
self-actuation, is a concern when many TMGs are placed next to one another 
with a small separation between them. At lower temperatures, when the 
frequency, stroke, and power output are already low, the separation does not 
appear to have any effect when two devices with near tolerances are run in 
parallel under the same magnet. Even at an extremely close spacing of 
300 um, the TMGs don't seem to have any coupling effects. At higher Tsource 
values, thermal cross-coupling effects are observed between devices when 
the separation between the devices is reduced below 1 mm. Correlation in 
performance between devices is observed for some of the instances suggest- 
ing it is possible for devices to develop a synchronization in their operation 
where its performance is dependent on the neighboring device due to thermal 
coupling. 


8.1.4 Expanding the range of operation temperature 


TMGs with Ni-Mn-Ga has shown the potential of the concept of resonant 
self-actuation but this is limited to temperatures well above 100 °C due to the 
Te of the material. Thus, in order to expand the range of heat source tempera- 
ture under which resonant self-actuation could be achieved, Gadolinium film 
and La-Fe-Si based alloys are used. Dzekan et al.!6°16* gives an overview of 
novel thermomagnetic materials for power generation. Section 5.2.5 has 
presented the effect of Tc on the minimum heat source temperature. 


Using the Gd film under optimum conditions, a large power output of 3 uW 
resulting in a power per footprint of 23.8 1 W/cm? was achieved at a Tsource of 
65 °C and Tamb of 11 °C. This opens up a wide variety of applications as the 
temperature required for achieving resonant self-actuation is well below the 
boiling point of water. But this investigation revealed strict boundaries of 
temperatures for both the ambient (Tamb) and heat source (Tsource) where the 
TMG could maintain resonant self-actuation. Since the Te of Gd film was 
20 °C, maximum Tamb was limited to 19°C above which the TMG stops 
operation. Even though, the TMG is able to sustain resonant self-actuation at 
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a Tsource of 40 °C, the minimum value of Tsource increases to 50 °C as the Tamb 
is increased to 19 °C.But even at the maximum Tamb of 19 °C and minimum 
Tsource Of 50 °C, the TMG is capable of generating 1 uW resulting in a power 
per footprint of 8.7 uW/cm?. While under optimum ambient conditions the 
value reaches 10 uW/cm? even at the minimum Tsource of 40 °C. 


Due to the large magnetization and low Te of La-Fe-Si based alloys, these 
materials are suitable for developing TMGs with high performance near room 
temperature. For the investigation presented in this work, La-Fe-Si ribbon 
fabricated by melt spinning and later hydrogenated are used. This TMG 
generated a very large power output of 4.2 uW while operating at a Tsource of 
90 °C despite having very low value of hs. This low value of hr is the result of 
the fabrication process which makes the surface very rough and thus mini- 
mizing the contact area between TM film and magnet during operation. But 
in order to reduce the Tsource under which resonant self-actuation could be 
achieved for TMGs with La-Fe-Si and to be able to generate maximum 
output, further optimizations are required. By using Gd, La-Fe-Si-H and Ni- 
Mn-Ga, a whole range of TMGs are presented. These devices make use of 
Tsource starting from 40 °C all the way up to 170 °C. The TMG based on Gd 
operated from 40 °C to 70 °C, while the TMG based on La-Fe-Si-H operated 
between 70°C and 100 °C. The Ni-Mn-Ga operates between 110°C and 
170 °C. 


8.2 Benchmarking 


All demonstrators in this work are fabricated with similar dimensions with an 
overall lateral footprint of the final device in the range from 10 mm? to 
14 mm?. All TMGs presented are compared with each other and then a few of 
them have been compared to other demonstrators from literature. 


8.2.1 TMGs within this work 


The peak performance of the device at higher Tsource values and the perfor- 
mance at the minimum operating temperatures are shown in Figure 103. The 
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maximum power per footprint of 50 uW/cm? is observed for a TMG with 
40 um thick Ni-Mn-Ga film at Tsource of 170 °C and then drops to a value of 
28 uW/cm? when the Tsource is reduced to 130 °C. At lower Tsource values of 
90 °C, the TMG with La-Fe-Si-H as the active TM film generates a power 
per footprint of 38 uW/cm?. This marks a milestone in TM energy generation 
as the waste heat from water boilers is now accessible to this technology as 
the Tsouree is below 100°C. The power per footprint value reduces to 
20 uW/cm? when the Tsource is reduced to 70 °C. Gadolinium based TMG is 
able to generate a power per footprint of 24 uW/cm? from an even lower 
Tsource of 65 °C. At a heat source of 40 °C which is the minimum for all the 
TMGs discussed in this work, the TMG with Gd film as TM film is still capa- 
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Figure 103 Comparison of thermomagnetic generators presented as part of this work (a) power 
per footprint of devices, (b) temperature change within active material of the 
device, (c) frequency of operation of the devices. 
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ble of gen-erating 10 uW/cm?. This is also another milestone in the develop- 
ment of TMGs near room temperature as we are getting closer to human body 
temperature. All the demonstrators presented in this work are operating with 
high frequencies. The highest frequency of 146 Hz is observed for the TMG 
with 40 um Ni-Mn-Ga film. La-Fe-Si-H based TMG is operating with a 
frequency of 137 Hz, but this is restricted due to the low value of hr. Once 
this is improved, it does have the potential to operate at close to 200 Hz 
(assumption based on the natural frequency of cantilever assembly) at a low 
Tsouree of 90 °C. The TMG with Gd film shows that it is possible to sustain 
resonant self-actuation with a frequency of 107 Hz and 54 Hz at 70 °C and 
40 °C Tsource values respectively. 


Figure 103 also shows the AT in the TM films for the respective TMG 
demonstrators. It can be observed that AT is higher in thinner TM films and is 
lower as the thickness of the sample is increasing resulting in the large 
thermal mass of the thick films. The Ni-Mn-Ga films show much larger AT 
even at higher thicknesses compared to Gd and La-Fe-Si-H films. Ni-Mn-Ga 
film shows a AT between 3 °C to 10°C for the TMG with TM film thick- 
nesses between 40 um to 10 um. Gd film shows a AT of 0.45°C at a Tsource 
value of 40 °C and the La-Fe-Si-H shows the least AT value of 0.22 °C at a 
Tsource of 70 °C. Even at the higher Tsource of 90 °C, the AT is only 0.43 °C. 
Unlike, frequency and power output, the AT required for operation should be 
as low as possible. This is also dependent on the thermal mass (volume of the 
material). But if the material volume is assumed to be the same, then a lower 
AT requirement for resonant self-actuation is desirable as it would also result 
in a lower heat intake requirement. 
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Table 10 Performance of different TMGs are presented. 


Material Thick- | Tsource | AT Size | Pow-| Power Fre- 
ness of er per quency 
film footprint 


(um) CCO | (°C) | (mm?) | uw) nn (Hz) 


Ni-Mn-Ga | 40 130 2 10.6 |3 28.3 104 


40 170 2.9 5.3 | 50 146 


20 130 2.7 12 2.9 | 24.2 98 


20 170 5.4 4 32.8 136 

10 130 5.9 13.4 |1 7.8 79 
10 170 10 3 23 114 

Gadolinium | 40 40 1 12.6 |1.26 |10 64 
65 0.9 3 23.8 107 
La-Fe-Si-H | 50 70 0.2 11 2.2 |20 127 
90 0.4 4.2 |38 138 


8.2.2 Benchmarking 


TMGs with 40 um Ni-Mn-Ga film, Gd film, and La-Fe-Si-H film are selected 
for comparison with other similar TMG demonstrators presented in the 
literature. Several published examples are considered that use the TM energy 
conversion principle and TM films to convert thermal energy into electrical 
energy. Many more potential energy-generating systems ha-ve been pro- 
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posed, but they primarily establish the actuation state, which transforms 
thermal energy into mechanical energy.!!'!® It can be seen that in terms of 
power per footprint the 40 um Ni-Mn-Ga is the peak performing TMG. The 
closest competition is provided by the TMG presented by Waske et al.!!° But 
even though the power output is considered without taking into account the 
energy required for powering the pump used to circulate the heat transfer 
fluid, the power per footprint is 32 uW/cm? while operating at a frequency of 
1 Hz. The TMG with 40 um Ni-Mn-Ga film generates 56 % more power per 
footprint while operating at 146 Hz. 
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Figure 104 Comparison with similar devices from literature, only best performing thermomag- 
netic generators from this work are included, (a) power per footprint of devic- 
es, (b) temperature change within the active material of the device, (c) frequen- 
cy of operation of the devices. Waske et. al.,''" Ahmim et. al.,''* Rodrigues et. 
al.,'"® Chun et. al., Lallart et. al.,''’ Srivastava et. al.,'® Chun et. al.,!!? Guel- 
tig et. al.,”” Carlioz et. al.'”' Adapted image from Joseph et al.'”? 
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It can also be seen that all demonstrators presented in this work are among 
the smallest devices. The frequency of the operation is also among the 
highest with no other device being in even the same order of magnitude. The 
AT in the material under operation is among the lowest. Both TMGs with La- 
Fe-Si-H film and Gd film show a AT that is more than an order of magnitude 
lower than the other devices. 


Table 11 Performance of TMG from literature sources. 


Source Tsource | AT | Size Power | Power per | Frequency, 
footprint 


(°C) | (°C) | (mm?) | (uw) | (HW/em’) | (Hz) 


Waske et al.!!° 50 30 3900 1240 31.8 1 
Ahmin et al.!'* 56 35 3316.6 | 4.2 0.13 0.12 
Rodrigues et al.!!5 | 65 25 |8800 | 481.4 | 5.47 0.011 
Chun et al.!”° 70 80 10000 | 158 1.58 3 
Lallart et al.!!” 115 23 28 0.0058 | 0.021 0.35 
Srivastava et al.!°8 | 250 10 |176.7 |50 28.3 0.05 
Chun et al.!”° 50 60 | 500 80 16 9 


Gueltig et al.” 160 15 17 2.4 14.12 84 
8 


Carlioz et al.!?! 50 10 |6079 \3.2 0.053 3.1 
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9 Conclusions and outlook 


This section concludes the work on power generation by resonant self- 
actuation with a summary of the presented work and an outlook to further 
develop resonant self-actuated thermomagnetic generators. 


9.1 Conclusions 


With the increasing demand for connected IoT (internet of things) devices 
which are forming a grid to sense its environment conditions, remote power 
generation is ever more in demand. The ability to generate useful amounts of 
power on site would enable widespread adaption. Since this would reduce the 
maintenance efforts, this development is also beneficial to the economy in 
long term. Since much of the waste energy is in the form of low-grade 
thermal energy in the environment which would otherwise be lost, a technol- 
ogy to harvest this energy has a widespread interest. Thermal energy harvest- 
ing using thermoelectric devices utilizing the Seebeck effect is popular at 
higher temperatures, but when it comes to low heat source temperature 
regime below 200 °C and miniature size of a few millimeters, this method of 
thermal energy harvesting lacks the power density required to stay competi- 
tive. 


Using thin films and the recently established concept of resonant self- 
actuation, this thesis presents several TMGs operating with high power 
output per footprint. In addition an extensive study on device performances 
and the factors affecting the performance is presented based on LEM simula- 
tions. 


The Key results can be summarized as follows; 


1. LEM model and parameter study: Different domains of the model are 
dedicated to the different physical principles that drive the TMG. The 


191 


9 Conclusions and outlook 


192 


magnetic part of the model is utilized to derive the magnetic field and 
magnetic field gradient of the permanent magnet at various positions. 
The mechanical domain of the TMG is used to model and understand 
the oscillation behavior due to the stiffness, mass, damping, and mag- 
netic force. The magnetic force is dependent on the magnetic field, 
which causes the magnetization of the TM film and the magnetic field 
gradient which causes the magnetic force. The thermal domain of the 
LEM simulation models the heat flow through the different layers using 
thermal resistances and capacitances. The validated LEM gives a valua- 
ble conception of the thermal processes within the TMG and supports 
the development of novel thermomagnetic generators that operate in 
resonant self-actuation mode. This allows the model to derive the tem- 
peratures of different components of the TMG during its transient op- 
eration. Since the TMG never reaches thermal equilibrium and always 
operates in transient mode, understanding the thermal processes in the 
TMG is crucial to achieving resonant self-actuation mode. The electrical 
domain of the LEM simulation is used to model the energy generation 
from the mechanical oscillation that is converted to electrical energy 
with the pick-up coil that moves in the changing magnetic field gradi- 
ent. All domains of the LEM simulation are solved simultaneously so 
that mutual interactions between the different cycles can be modeled 
appropriately. The simulation model is validated by comparing, stroke, 
frequency and power output ofthe TMG from experiments. To design a 
TMG with a TM film, initially, the hr is postulated as this value is very 
hard to measure accurately. After fabricating the designed TMG, the 
device is tested and then the simulation is readjusted to match the per- 
formance of the TMG by matching the hr. 


Different parameters affect the performance of the TMG. The material 
of the magnet and its dimension affect the magnetic field and magnetic 
field gradient generated by the permanent magnet. It is also shown how 
having a higher aspect ratio (between the height of the magnet along the 
easy axis and cross-sectional area) would result in a higher magnetic 
field gradient which would increase the magnetic force, but would also 
result in a lower magnetic field at a distance, as a result, causing a re- 
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duction in magnetization of the TM film. Hence choosing the proper 
dimensions of the permanent magnet is essential in achieving optimum 
performance of the TMG as the maximum stroke, stiffness (hence fre- 
quency), and power output would be dependent on this. The effect of 
different thermal processes on the performance of the TMG is also dis- 
cussed. It is shown there is a critical value for both hr (8 kW/m’K) and 
Rp (5 K/W) between the TM film and cantilever that affect the heat in- 
take into the material and heat dissipation from the material respective- 
ly. Above the critical value increasing the hr increased the output power 
and frequency of the TMG. Similar observations are made for increas- 
ing Rv. But after reaching the peak value for the AT, the AT is slightly 
showing a reduction suggesting that further increase might cause a re- 
duction in power output and there would be a maximum critical value 
for the bonding layer R» above which resonant actuation would fail. But 
the investigated Ru of 45 K/W corresponds to a bonding layer thickness 
of 40 um for material with thermal conductivity of 0.225. This is al- 
ready four times the size of the 10 um Ni-Mn-Ga film used for investi- 
gation suggesting the higher limit is beyond reach in practical applica- 
tions. 


The effect of cantilever stiffness is also studied in detail, and it is found 
that increased stiffness would result in higher power output and fre- 
quency due to the increased natural frequency of cantilever assembly. 
But beyond a critical stiffness, the device stops operating as the device’s 
mechanical oscillation cycle does not match with the TMG’s thermal 
cycling. The heating and cooling cycles of the TMG needs to be in bal- 
ance with the oscillation cycle in order to sustain resonant self- 
actuation. The pick-up coil turns in the device are also affecting the per- 
formance. It is possible to optimize the number of turns in the pick-up 
coil according to the Tsource, magnetic force, and stiffness of the cantile- 
ver to achieve the maximum power output. This is the result of optimiz- 
ing the electromagnetic damping to generate maximum power output. 
The Tc of the TM film used as the active material in the TMG also af- 
fects the performance defines the operating temperature range of the 
TMG. The TMG with high T. of the TM film would require higher 
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Tsouree While the TMG with lower Ts of the TM film would require lower 
Tsource Values. As the Tc is approaching the Tamb, additional effects due to 
lack in dissipation capability are observed that cause a reduction in 
power output. Further an increase in magnetization is required to over- 
come this limitation. The load resistance (RL) connected to the TMG al- 
so affects the performance as this regulates the current through the pick- 
up coil. The current in the pick-up coil causes electromagnetic damping 
at the cantilever tip. Thus, the Rr also needs to be matched to the TMG 
to achieve maximum output. It is also observed that increasing Rı 
would cause an increase in power generated due to the reduced electro- 
magnetic damping, resulting in higher frequency and stroke. But the 
transfer efficiency decreases as the value of the Rx is increased further 
(more than the internal resistance of the pick-up coil) resulting in an op- 
timum value of RL. 


2. Upscaling: The upscaling of the TMG is investigated to achieve high- 
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er power output. 


a). Increase in active material volume: The TMG is successfully up- 
scaled from generating 2.4 uW to 5.3 uW resulting in a gain of 3.4 
times power output per footprint with values reaching 50 uW/cm? by 
increasing the active material volume. The effects of upscaling are 
studied in detail using experiments and LEM simulation to under- 
stand the thermal and mechanical effects of scaling. The upscaling is 
limited by the lack of homogeneous thin films of Ni-Mn-Ga films 
beyond 10 um thickness, demanding the fabrication of thin film 
stacks prepared by stacking Ni-Mn-Ga using a conductive bonding 
layer. The thermal mass of the stack is causing the AT in the materi- 
al to drop and 40 um is the limit beyond which it is not possible to 
increase the power output due to the additional passive thermal load 
of the conductive bonding layer. 


b). Increasing the lateral size: The lateral size of the effective area of 
power generation of the TMG could be upscaled by adding multiple 
TMGs in parallel. Adding multiple TMGs in parallel under a single 
magnet surface may result in coupling effects between the devices. It 
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is observed that only thermal coupling effects affect the performance. 
For the TMG using 10 um Ni-Mn-Ga films, the coupling effects are 
observed only at higher Tsource values above 150 °C. This shows the 
tolerance ofthe TMGs to changes in the ambient condition and the po- 
tential to upscale in terms of the number of devices operating in paral- 
lel close to each other. The thermal coupling effects are caused due to 
the inability of the TMGs to dissipate the heat that is taken in during 
contact effectively because of the neighboring devices. It is possible 
for the TMGs to operate in correlation to each other, with the TMG 
increasing its output as soon as the neighboring TMG stops operation 
and as soon the neighboring TMG resumes operation, the primary de- 
vice reduces its power output to accommodate the thermal coupling. 


3. Expanding the range of operation temperature: Operation near 
room temperature is achieved using Gd and La-Fe-Si-H as the active 
TM films. The TMG using La-Fe-Si-H film with Te of 48 °C is able to 
operate in resonant self-actuation mode even at a Tsource of 70 °C. The 
TMG using a Gd film, on the other hand, with a Te of 20 °C enable op- 
eration in resonant self-actuation at a minimum Tsource of 40 °C. Since 
the change in magnetization of these two materials is different, it is pos- 
sible to understand how both, the gradual change of magnetization and 
the abrupt change of magnetization, could be utilized to develop devic- 
es, which can either sustain resonant self-actuation in a wide range of 
ambient conditions or devices that can perform with large power output 
and very high frequency at a small range of temperature, respectively. 
The TMGs using Gd film is generating 3.1 uW at its peak performance 
state with a Tsource of 65 °C while also generating 1 uW at 40 °C. The 
TMG using a La-Fe-Si-H film generates around 4.2 uW at a Tsource of 
90 °C and is also capable of generating 2.2 uW of power at a Tsource of 
70 °C. The performance of the TMG with La-Fe-Si-H film is limited by 
the surface roughness of the material as a result of the fabrication pro- 
cess. 


With this thesis work, TMGs operating within a range of 40°C to 
170 °C for Tsource values are developed. A LEM simulation tool is de- 
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veloped for design and optimization of film based TMGs for various 
other materials and layouts. These preliminary results point out a route 
for future investigations, which are discussed in more detail in section 
9.2. 


9.2 Outlook 


Developing thermomagnetic generators operating near room temperature 
comes with a lot of challenges. Although many of these challenges have been 
tackled during the course of this thesis, further work needs to be undertaken 
to achieve more TMGs, e.g., that may utilize body heat as the heat source. 
For optimization, the procedures below, in particular, need to be taken into 
account: 


e Extending the LEM simulation model for impact of coupling effects in 
TMGs operating in parallel 


e Optimizing bi-stable TMGs by integrating a pyroelectric element 


e Improving the surface roughness of La-Fe-Si-H samples to increase 
power output 


e In vacuum operation could be achieved if support of convection could 
be completely avoided. 


Developing novel materials with high magnetization (M), large mag- 
netization changes with respect to temperature (AM/AT) and Te 
near room temperature 


Continuous innovation in the field of TM energy generation is required to 
achieve near-room temperature resonant self-actuation in TMGs which are 
capable of generating a large amount of power from a small area. Since the 
absolute power output of the TMG would be small considering its miniature 
dimensions, it is inevitable that multiple devices need to be operated side by 
side to achieve sufficient energy output required to continuously power a 
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sensor node. To optimize the power output and thermal performance of the 
TMGs, the thermal coupling between the devices need to be investigated 
experimentally and by using a LEM simulation model. 


TMG with La-Fe-Si-H film shows a lot of potential, and the device design 
suggests high frequency operation that is attainable with good surface contact 
with the magnet. Improving the surface roughness of the La-Fe-Si-H sample 
would increase the hr between the film and the magnet. This would increase 
the power output many folds and also improve the minimum Tsource at which 
the device is capable of sustaining resonant self-actuation. The surface 
roughness could be improved either by an additional layer of conductive 
material and then polishing the material or by choosing an alternative fabrica- 
tion process for the material. Adding manganese would allow further tuning 
of the material properties of the material to achieve operation at even lower 
Tsource values. In particular, accessing Tsource less than 37 °C by TMGs based 
on La-Fe-Si-Mn-H films would open up a large field of new applications 
particularly in the medical field. 


Currently, the thermal dissipation in the TMG is actively controlled by 
utilizing thermal convection to ambient air. This is achieved as the convec- 
tion is increased many folds as soon as the device is rapidly oscillating and is 
negligible when the device is in contact as the cantilever is not moving. If 
this active component of the TMG is integrated into the cantilever design to 
avoid support of convection, then it is possible to develop TMGs capable of 
operating in a vacuum by relying solely on conductive dissipation of heat 
energy. This would allow encapsulation of the TMG and increase the thermal 
performance many folds by taking the effects of thermal convection away. 
Multiple devices can also be operated at a very close distance as the thermal 
coupling effects could be completely avoided by avoiding convective dissipa- 
tion. 


Further, new novel ferromagnetic materials with high value of magnetization 
in ferromagnetic state and low T. near room temperature need to be devel- 
oped. It would also be possible to integrate materials showing low or zero 
hysteresis first-order transformation into the present design to achieve reso- 
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nant self-actuation which could potentially generate large amounts of power 
owing to its abrupt first-order transformation. 
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Table A 1 Important values of temperatures for the TMG upscaling 


Device film thickness: 40 um 


Source temperature | 130 140 150 
Sink temperature 87 87 87 
(Teant) 

Tiop (max) 91 91.1 91.2 
Tiop (min) 88.9 88.8 88.6 
Toot (max) 90.3 90.4 90.4 
Toot (min) 88.8 88.7 88.6 
Temperature change 2.1 2.3 2.6 
(Trop) 

Temperature change 1.5 1.7 1.8 
(Toot) 


Device film thickness: 20 um 


Source temperature | 130 140 150 
Sink temperature 91 91 91 
Tiop (max) 96.5 96.6 96.7 


Trop (min) 93.8 92.9 91.5 


Appendix 


Toot (max) 96.4 
Toot (min) 93.8 


Temperature change | 2.7 
(Trop) 


Temperature change | 2.6 
(Toot) 


Device film thickness: 10 um 
Source temperature | 130 
Sink temperature 93 

Tiop (max) 100.9 
Tiop (min) 94,9 


Temperature change | 6 
(Trop) 
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Figure A 1 The magnetization cycle in TMG, with values corresponding to operation in resonant 


self-actuation mode of TMG with 40um Ni-Mn-Ga film. Adapted from Joseph 


et al.” 


Rt has a significant impact on TMG performance. When Rı is matched to the 
internal resistance of the pick-up coil, which is roughly 220 Ohm, the transfer 
efficiency is greatest. However, as Rx is raised, the power increases owing to 
reduced electromagnetic damping, leading in a greater frequency and stroke 
of operation. Simultaneously, the transfer efficiency decreases, and a distinct 
maximum occurs at 500 Ohm. Even if the frequency and stroke are still 
growing at higher Rı, the decrease in electromagnetic damping is insufficient 
to compensate for the loss in transfer efficiency. 
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Figure A 2 The mechanical and electrical performance of the TMG device with Ni-Mn-Ga 


10um as a function of load resistance (Rı) and heat source temperature (Tsource). 
The supplied color code can be used to calculate numerical values. (a) Stroke 
during resonant self-actuation, (b) oscillation frequency, (c) electrical power 


output. 
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Figure A 3 The mechanical and electrical performance of the TMG device with Ni-Mn-Ga 


20um as a function of load resistance (Rı) and heat source temperature (Tsource). 
The supplied color code can be used to calculate numerical values. (a) Stroke 
during resonant self-actuation, (b) oscillation frequency, (c) electrical power 


output. 
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Figure A 4 The mechanical and electrical performance of the TMG device with Ni-Mn-Ga 
40um as a function of load resistance (Rı) and heat source temperature (Tsource). 
The supplied color code can be used to calculate numerical values. (a) Stroke 
during resonant self-actuation, (b) oscillation frequency, (c) electrical power 
output. 
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Figure A 5 Effect ofresistance and separation on performance of TMG at center in three device 
array with 10m Ni-Mn-Ga. (a) Power output of TMG, (b) frequency of opera- 
tion of TMG. 


243 


Appendix 


Power (IWW) b) Frequency 


600 

0.6 75.9 

J 74.8 

. 500 73.6 

. 725 

. 713 
400 

E 70.2 

. 69.1 

300 679 

i 66.8 
200 

1 2 3 4 5 1 2 3 4 5 


Separation (mm) Separation (mm) 


E 
c 
Q 
8 
© 400 
© 
© 
7 
® 
X 


Resistance (Ohm) 


Figure A 6 Effect ofresistance and separation on performance of TMG in three device array with 
20um Ni-Mn-Ga film. (a) Power output of TMG, (b) frequency of TMG opera- 


tion. 
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Figure A 7 The mechanical and electrical performance of the Gd TMG device as a function of 
load resistance (Ri) and heat source temperature (Tsource) at 11 °C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output. 
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Figure A 8 The mechanical and electrical performance ofthe Gd TMG device as a function of 
load resistance (Rı) and heat source temperature (Tsource) at 19 °C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output determined 
from the induced current of the pick-up coil. 
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Figure A 9 The mechanical and electrical performance of the Gd TMG device as a function of 
load resistance (R) and heat source temperature (Tsource) at 17°C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output determined 
from the induced current of the pick-up coil. 
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Figure A 10 The mechanical and electrical performance of the Gd TMG device as a function of 


load resistance (Ri) and heat source temperature (Tsouree) at 15 °C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output determined 
from the induced current of the pick-up coil. 
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Figure A 11 The mechanical and electrical performance of the Gd TMG device as a function of 
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load resistance (Ri) and heat source temperature (Tsource) at 13 °C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output determined 
from the induced current of the pick-up coil. 
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Figure A 12 The mechanical and electrical performance of the Gd TMG device as a function of 
load resistance (RL) and heat source temperature (Tsource) at 9 °C ambient tem- 
perature were investigated. The supplied color code can be used to calculate 
numerical values. (a) Cantilever front stroke during resonant self-actuation, (b) 
cantilever front oscillation frequency, (c) electrical power output determined 
from the induced current of the pick-up coil. 
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Figure A 13 The thermomagnetic cycle of the Gd-film TMG was simulated using LEM for Tsource 
of 40 °C and Tamb of 11 °C. In an ideal world, the cycle would consist of four 
phases shown by stages I-IV. But in the actual cycle, stages III, IV and I are 
combined into an continuous path as the device is not maintaining adiabatic 
conditions in stage III and I, nor does it have a heat sink at stage IV. We evalu- 
ate the material's efficiency using the simulated thermomagnetic cycle, which 
is represented by the ratio of output magnetic energy EMag to input heat ener- 


; EMa 
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Table A 2 LEM simulation values of major thermal properties of TMG 
TMG based on Gadolinium film 
Ambient temperature: 9°C 


Source 40°C 45°C 50°c | 55°C | 60°C 65°C | 70°C 
temperature 
(Tsource) 
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